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

GATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_pae_pgi/1979134‑1979424                                                                                                                                                  |                                                                                                                                                  aatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                                >  1:1394610‑M2/146‑149 (MQ=255)  atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                               >  2:592137‑M2/145‑149 (MQ=255)     ggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                                                                            <  1:277536‑M2/8‑1 (MQ=255)      gtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                                                                           >  1:827767‑M2/141‑149 (MQ=255)      gtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaataTATAAGTCa                                                                                                                                           >  1:541726‑M2/141‑149 (MQ=255)                ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                                                                 >  2:113344‑M2/131‑149 (MQ=255)                     taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                            >  1:80556‑M2/126‑149 (MQ=255)                         atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                        >  2:694417‑M2/122‑149 (MQ=255)                          tttttccgaaagcgacgcaactggactcgtcactacacgcacatacaccggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                       <  2:898812‑M2/29‑1 (MQ=255)                              tcctaaatcgacgcaaatgtcatcgtcactccacgcaaacacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                                   <  1:530433‑M2/33‑1 (MQ=255)                                ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                 >  1:1604412‑M2/115‑149 (MQ=255)                                  aaatcgacgcaactgtactggtcactagacggacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactctattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                                                               <  2:327172‑M2/37‑1 (MQ=255)                                        acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                         >  2:201851‑M2/107‑149 (MQ=255)                                         cgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                                                                                        >  1:535622‑M2/106‑149 (MQ=255)                                          gcaactgtgcgcggcacgagacggacatacgacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                       <  2:223593‑M2/45‑1 (MQ=255)                                                gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                 <  1:113344‑M2/51‑1 (MQ=255)                                                ggactcgtcactacacgaacagacaacggaggggggctgagatttccgatcaggcgtaatgcagatcaaacaaaacactcgactacataacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                 <  2:748562‑M2/51‑1 (MQ=255)                                                atactcttccctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                                 <  2:1210729‑M2/51‑1 (MQ=255)                                                 tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                                <  2:1026573‑M2/52‑1 (MQ=255)                                                 tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                                <  1:1012975‑M2/52‑1 (MQ=255)                                                     cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                            <  1:562372‑M2/56‑1 (MQ=255)                                                     cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                            <  1:790642‑M2/56‑1 (MQ=255)                                                      gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                           >  1:1264576‑M2/93‑149 (MQ=255)                                                        cactacacgcacatacaacggagggaggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattactgaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                         <  2:1010910‑M2/59‑1 (MQ=255)                                                         actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                                                                        >  1:149264‑M2/90‑149 (MQ=255)                                                          ctacacgcacctagaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAACAAGTATTTcgctacgcta                                                                                       <  2:151256‑M2/61‑1 (MQ=255)                                                          ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTtcgctaagctc                                                                                       >  1:161477‑M2/89‑148 (MQ=255)                                                             cacgcacatacaacggaggggggctgcgattttcaagaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTGTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                                                                    <  1:1138823‑M2/64‑1 (MQ=255)                                                              acgcacataccacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                   <  1:592137‑M2/65‑1 (MQ=255)                                                                 cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                                                                >  2:326269‑M2/82‑149 (MQ=255)                                                                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                             >  1:665776‑M2/79‑149 (MQ=255)                                                                    atacaacggagggaggctgcgattttcaataatgcgtgatgtcagtcatcaaaaaccctccataactgaacagaaatgGATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                             >  2:1151069‑M2/79‑149 (MQ=255)                                                                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                            >  1:49821‑M2/78‑149 (MQ=255)                                                                        aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                         >  2:1426111‑M2/75‑149 (MQ=255)                                                                          cggaggggggcttcaattttaaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGGCATACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                       <  2:515483‑M2/77‑1 (MQ=255)                                                                             aggggggccgcgatttaaaaaaatgcgtgatgcagatcacaaaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTAAAAAAAAAAAAAAAGTTTTTCGCTACGCTAGCAAACAAAAATGAAACAg                                                                    <  2:1394610‑M2/80‑1 (MQ=255)                                                                                   gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                                                              >  1:1110608‑M2/64‑149 (MQ=255)                                                                                     tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                            <  1:694417‑M2/88‑1 (MQ=255)                                                                                         attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGGTACGCTAGCAAGCAAAAATGTAACAAGAATAATCGaaa                                                        >  2:116615‑M2/58‑149 (MQ=255)                                                                                           tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg                                                      >  2:368287‑M2/56‑149 (MQ=255)                                                                                             tcaataatgcgtgatgcagatcacacaaaccactcaactacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTGCAATCTGGTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg                                                    <  1:1305464‑M2/96‑1 (MQ=255)                                                                                               aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                                  >  2:1419847‑M2/52‑149 (MQ=255)                                                                                               aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTTGGGTGTATTTCAATCTGTTAAACAGTTTTTCGCTACGCTAGAAAGAAAAAATGAAACAGTAATAATCGCAATCGGAt                                                  >  2:430276‑M2/52‑149 (MQ=255)                                                                                                 taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                                >  2:224102‑M2/50‑149 (MQ=255)                                                                                                   atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTg                                              >  2:503516‑M2/48‑149 (MQ=255)                                                                                                   atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTg                                              >  2:504487‑M2/48‑149 (MQ=255)                                                                                                       gtgaagcagagcacacaaaacaatgaaatacttaacacaaataTATAAGTCAGACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                          <  1:116615‑M2/106‑1 (MQ=255)                                                                                                            gcagatcacacaaaccactcaattacttaacataaatgTATCAGGCATACTTTCGTATTGGGTGAAGTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACCGTc                                     >  2:1547429‑M2/39‑149 (MQ=255)                                                                                                            gcagatcacacaaaacactcaattacttaacataaattTATAAGTCATACTTTTATTTTGGCTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCCACCACGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAgtg                                     <  2:1224209‑M2/111‑2 (MQ=255)                                                                                                                   acacaaaacactccattacgtaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                              <  1:210133‑M2/118‑1 (MQ=255)                                                                                                                   acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                              <  1:326269‑M2/118‑1 (MQ=255)                                                                                                                   acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                              >  1:1101089‑M2/32‑149 (MQ=255)                                                                                                                   acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                              >  1:567912‑M2/32‑149 (MQ=255)                                                                                                                       aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                             >  1:171621‑M2/28‑146 (MQ=255)                                                                                                                       aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAc                             <  2:171621‑M2/119‑1 (MQ=255)                                                                                                                         aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAcc                        >  2:1130240‑M2/26‑149 (MQ=255)                                                                                                                             ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                    <  2:149264‑M2/128‑1 (MQ=255)                                                                                                                             ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                    >  1:747695‑M2/22‑149 (MQ=255)                                                                                                                              tcaattagttaacataaatgTATAAGTCATACTATTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                   <  2:1262110‑M2/129‑1 (MQ=255)                                                                                                                                aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCTCTACGCTAGCAAGCAAACATGAAACAGGAATAATCTAAATGCGATGTTGCGCACAGTCAAAATAACTCACCATAaata                 >  2:773384‑M2/19‑149 (MQ=255)                                                                                                                                   tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc              <  1:503516‑M2/134‑1 (MQ=255)                                                                                                                                   tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc              <  1:504487‑M2/134‑1 (MQ=255)                                                                                                                                    acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                             <  2:884665‑M2/103‑1 (MQ=255)                                                                                                                                    acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                             >  1:884665‑M2/15‑117 (MQ=255)                                                                                                                                     cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat            <  2:1264576‑M2/136‑1 (MQ=255)                                                                                                                                     cttaacataaatgTATAAGTCATACTTTTGTGTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAATGTCAGTTGAGCAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAAAAatcat            >  1:1578245‑M2/14‑149 (MQ=255)                                                                                                                                      ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAGAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc           >  1:1061619‑M2/13‑149 (MQ=255)                                                                                                                                           ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATTATCTGCta      >  1:901147‑M2/8‑149 (MQ=255)                                                                                                                                           ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta      >  1:1083057‑M2/8‑149 (MQ=255)                                                                                                                                               atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa  >  1:1103647‑M2/4‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_pae_pgi/1979134‑1979424

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: ‑