Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_pae_pgi	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_pae_pgi	1978503	1979278	776	71 [0]	[0] 70	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 1978502	0 (0.000)	67 (0.770)	55/266	0.7	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913_3_pae_pgi	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_pae_pgi/1978354‑1978648                                                                                                                                                      |                                                                                                                                                   tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                                                                    >  1:1157629‑M1/1‑149 (MQ=255)     aaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                                                                >  1:1394610‑M1/1‑145 (MQ=255)      aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               >  2:592137‑M1/1‑144 (MQ=255)         ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                            <  1:277536‑M1/149‑9 (MQ=255)          gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                           >  1:827767‑M1/1‑140 (MQ=255)          gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATAtataagtca                                                                                                                                           >  1:541726‑M1/1‑140 (MQ=255)                    ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                                 >  2:113344‑M1/1‑130 (MQ=255)                         tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                            >  1:80556‑M1/1‑125 (MQ=255)                             aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                        >  2:694417‑M1/1‑121 (MQ=255)                              tttttCCGAAAGCGACGCAACTGGACTCGTCACTACACGCACATACACCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                       <  2:898812‑M1/149‑30 (MQ=255)                                  tCCTAAATCGACGCAAATGTCATCGTCACTCCACGCAAACACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                                   <  1:530433‑M1/149‑34 (MQ=255)                                    cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  1:1604412‑M1/1‑114 (MQ=255)                                      aaaTCGACGCAACTGTACTGGTCACTAGACGGACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCTATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                               <  2:327172‑M1/149‑38 (MQ=255)                                            aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                         >  2:201851‑M1/1‑106 (MQ=255)                                             cGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                        >  1:535622‑M1/1‑105 (MQ=255)                                              gCAACTGTGCGCGGCACGAGACGGACATACGACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                       <  2:223593‑M1/149‑46 (MQ=255)                                                    ggactCGTCACTACACGAACAGACAACGGAGGGGGGCTGAGATTTCCGATCAGGCGTAATG‑CAGATCAAACAAAACACTCGACTACATAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                 <  2:748562‑M1/147‑52 (MQ=255)                                                    gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                 <  1:113344‑M1/149‑52 (MQ=255)                                                    ataCTCTTCCCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                 <  2:1210729‑M1/148‑52 (MQ=255)                                                     tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                <  2:1026573‑M1/149‑53 (MQ=255)                                                     tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                <  1:1012975‑M1/149‑53 (MQ=255)                                                         cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                            <  1:790642‑M1/149‑57 (MQ=255)                                                         cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                            <  1:562372‑M1/149‑57 (MQ=255)                                                          gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                           >  1:1264576‑M1/1‑92 (MQ=255)                                                            cACTACACGCACATACAACGGAGGGAGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                         <  2:1010910‑M1/149‑60 (MQ=255)                                                             aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                                        >  1:149264‑M1/1‑89 (MQ=255)                                                              cTACACGCACCTAGAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaacaagtatttcgctacgcta                                                                                       <  2:151256‑M1/149‑62 (MQ=255)                                                              cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctc                                                                                       >  1:161477‑M1/1‑88 (MQ=255)                                                                 cacGCACATACAACGGAGGGGGGCTGCGATTTTCAAGAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttgtgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                                                                    <  1:1138823‑M1/149‑65 (MQ=255)                                                                  acGCACATACCACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                   <  1:592137‑M1/149‑66 (MQ=255)                                                                     cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                                                                >  2:326269‑M1/1‑81 (MQ=255)                                                                        aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                             >  1:665776‑M1/1‑78 (MQ=255)                                                                        aTACAACGGAGGGAGGCTGCGATTTTCAATAATGCGTGATGTCAG‑TCATCAAAAACCCTCCATAACTGAACAGAAATGgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                             >  2:1151069‑M1/1‑78 (MQ=255)                                                                         tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                            >  1:49821‑M1/1‑77 (MQ=255)                                                                            aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                         >  2:1426111‑M1/1‑74 (MQ=255)                                                                              cGGAGGGGGGCTTCAATTTTAAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataaggcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                       <  2:515483‑M1/149‑78 (MQ=255)                                                                                 aGGGGGGCCGCGATTTAAAAAAATGCGTGATG‑CAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattaaaaaaaaaaaaaaagtttttcgctacgctagcaaacaaaaatgaaacag                                                                    <  2:1394610‑M1/149‑81 (MQ=255)                                                                                       gCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                              >  1:1110608‑M1/1‑63 (MQ=255)                                                                                         tGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                            <  1:694417‑M1/149‑89 (MQ=255)                                                                                             aTTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcggtacgctagcaagcaaaaatgtaacaagaataatcgaaa                                                        >  2:116615‑M1/1‑57 (MQ=255)                                                                                               tttCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                      >  2:368287‑M1/1‑55 (MQ=255)                                                                                                 tCAATAATGCGTGATG‑CAGATCACACAAACCACTCAACTACTTAACATAAATGtataagtcatacttttgttttgggtgtattgcaatctggtaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                    <  1:1305464‑M1/149‑97 (MQ=255)                                                                                                   aataatGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                                  >  2:1419847‑M1/1‑51 (MQ=255)                                                                                                   aataatGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgatttgggtgtatttcaatctgttaaacagtttttcgctacgctagaaagaaaaaatgaaacagtaataatcgcaatcggat                                                  >  2:430276‑M1/1‑51 (MQ=255)                                                                                                     taatGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                                >  2:224102‑M1/1‑49 (MQ=255)                                                                                                       atGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggttgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                              >  2:503516‑M1/1‑47 (MQ=255)                                                                                                       atGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                              >  2:504487‑M1/1‑47 (MQ=255)                                                                                                           gTGAAG‑CAGAGCACACAAAACAATGAAATACTTAACACAAATAtataagtcagacttttgttttgggtgtatttaaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                          <  1:116615‑M1/149‑107 (MQ=255)                                                                                                                g‑CAGATCACACAAACCACTCAATTACTTAACATAAATGtatcaggcatactttcgtattgggtgaagttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaccgtc                                     >  2:1547429‑M1/1‑38 (MQ=255)                                                                                                                g‑CAGATCACACAAAACACTCAATTACTTAACATAAATTtataagtcatacttttattttggctgtatttcaatctgttaaaaagtttttcgctacgccaccacgcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtg                                     <  2:1224209‑M1/149‑112 (MQ=255)                                                                                                                        acacaAAACACTCCATTACGTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              <  1:210133‑M1/149‑119 (MQ=255)                                                                                                                        acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              >  1:1101089‑M1/1‑31 (MQ=255)                                                                                                                        acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              <  1:326269‑M1/149‑119 (MQ=255)                                                                                                                        acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              >  1:567912‑M1/1‑31 (MQ=255)                                                                                                                            aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                             >  1:171621‑M1/1‑27 (MQ=255)                                                                                                                            aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                             <  2:171621‑M1/146‑120 (MQ=255)                                                                                                                              aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                        >  2:1130240‑M1/1‑25 (MQ=255)                                                                                                                                  cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                    >  1:747695‑M1/1‑21 (MQ=255)                                                                                                                                  cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                    <  2:149264‑M1/149‑129 (MQ=255)                                                                                                                                   tCAATTAGTTAACATAAATGtataagtcatactattgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                   <  2:1262110‑M1/149‑130 (MQ=255)                                                                                                                                     aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaacatgaaacaggaataatctaaatgcgatgttgcgcacagtcaaaataactcaccataaata                 >  2:773384‑M1/1‑18 (MQ=255)                                                                                                                                        tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  1:504487‑M1/149‑135 (MQ=255)                                                                                                                                        tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  1:503516‑M1/149‑135 (MQ=255)                                                                                                                                         aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                             >  1:884665‑M1/1‑14 (MQ=255)                                                                                                                                         aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                             <  2:884665‑M1/117‑104 (MQ=255)                                                                                                                                          cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat            <  2:1264576‑M1/149‑137 (MQ=255)                                                                                                                                          cTTAACATAAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaatgtcagttgagcaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaaaaatcat            >  1:1578245‑M1/1‑13 (MQ=255)                                                                                                                                           ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttagaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc           >  1:1061619‑M1/1‑12 (MQ=255)                                                                                                                                                aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataattatctgcta      >  1:901147‑M1/1‑7 (MQ=255)                                                                                                                                                aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta      >  1:1083057‑M1/1‑7 (MQ=255)                                                                                                                                                    aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  >  1:1103647‑M1/1‑3 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATG‑CAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_pae_pgi/1978354‑1978648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑