Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	16,725	T→A	intergenic (+168/+26)	insL1 → / ← hokC	IS186 transposase/toxic membrane protein, small

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	16,725	0	T	A	100.0%	186.1 / NA	52	intergenic (+168/+26)	insL1/hokC	IS186 transposase/toxic membrane protein, small
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (37/15);  total (37/15)

GACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTT  >  NC_000913/16495‑16954                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       gacATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCACCCCTTAAGTTAGCGCTAAGGGG‑ATc                                                                                                                                                                                                                               <  1:2823/239‑1 (MQ=11)       accccgCCATCGCTGGATTCCCCCCCAAGAAGGCCCGGATCCAAAAAGAAGAACTACCTCGTTGGGGAGAAAAAAAAAAATGGTCATCTGGAGCTTAAGGGTGGCCATTCGTGGGAC‑AGTACCCCGGACACCCAACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCCCCGGGGCGTCGCATTCTTCAGATGGTTCACCCCTTAAGTTAGCGCTAAGGGG‑ACCACcccc                                                                                                                                                                                                                         <  2:569986/234‑1 (MQ=11)             ccATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGt                                                                                                                                                                                                                   >  1:630219/1‑239 (MQ=21)             ccATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGt                                                                                                                                                                                                                   >  2:902938/1‑239 (MQ=21)                        tttCCCCCCAATAAGCCCGGAACCCGAAAGGAAGATATAACTCGTTGTGTAAAAGAAAAAATAGGGTCTTCTGGTGCTTACCGGTTGCCATTCGTGGGAC‑ATTTCCCCTGACAGCCTACAAAACGCAATTGAAGAACGGGGCGCATCTTTTTACCGAGGCCCGGGGGCGTCTCATTCTCCAGATGGTTCACCCCTTAAGTTGGCCCAAATGGG‑ACCACCCCCCGCCGTTGCTCTTACt                                                                                                                                                                                                         <  2:777526/238‑1 (MQ=2)                                         ccGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGGGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTTGCATTCTTCAGATGGTTCAACCCTTAAGTTACCGCTAATGGG‑ATCACTCCCCGCCGTTCCACTTACTCGGATTCGTAAGCAGTGa                                                                                                                                                                                       >  2:236936/1‑239 (MQ=255)                                                         gaaCTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAAcc                                                                                                                                                                            >  2:181488/1‑234 (MQ=255)                                                                           aaaaaaCAAAAAGGGTCATCGGGGGTTTACAGGTGGCCATTCGGGGGAA‑AGTACCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGGGGCATCGTCTTAACGAGGCCCCGGGGCGCCGCATTCTTCAGATGGTTCACCCCTAAAGTTAGCGCTAAGGGG‑ATCCCTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGAtgtg                                                                                                                                                     <  1:902938/235‑1 (MQ=18)                                                                              tAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGGAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGaa                                                                                                                                                   >  1:1320660/1‑238 (MQ=255)                                                                              tAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGaa                                                                                                                                                   >  1:1176900/1‑238 (MQ=255)                                                                                        ggTCATCTGGAGCTTACAGGGGGCCATTCGTGGGAC‑AGTATCCCTGACACCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACGGGGGCGTCCCATTTTCCAGATGTTTCACCCCTAAAGTTAGCGCTAATGGG‑ATCCCCCCCCGCCGTTGCTCTTACTCg                                                                                                                                                                                                       <  1:645506/176‑1 (MQ=14)                                                                                                  aGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAg                                                                                                                                                                                             >  2:645506/1‑176 (MQ=255)                                                                                                     taacgggTGG‑CATTCGGGGGACAAGTACCCCGGACGCCCAAAAAACCGAAATGGAAAACCGGGGGCCATCTTCTAAAGGGGCCACGGAGGCGCAGCTTTCTCCAGATGGTCCACCCTTAAAGTTACCCCTAAGGGG‑ACCCCTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTtctc                                                                                                                               <  2:498537/230‑1 (MQ=255)                                                                                                       aCAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTAcc                                                                                                                         >  1:239385/1‑239 (MQ=255)                                                                                                              gCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGcgcc                                                                                                                   >  1:1151111/1‑238 (MQ=255)                                                                                                                               tATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCgg                                                                                                                                                          >  1:1217557/1‑183 (MQ=255)                                                                                                                               tATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTcaca                                                                                                                                            >  1:1196354/1‑197 (MQ=255)                                                                                                                                 cccttGACATCCTACAAAACCCAATGGAAGAACGCGGGCCATCCTCTTAACGAGGCCCGGAGCCGTCGCTTTCTCCAGATGTTCCACCCCTTAAGTGACCCCTAAGGGG‑ATCACTCCCCGCCGTTGCTCTTACTCgg                                                                                                                                                                                                      <  2:1077540/133‑1 (MQ=255)                                                                                                                                 ccccTGACGGCCTAAAAAACGCAATTGAAGAACGCGAGGCATCGTCTTACCGAGGCCCCGAGGCGCCCCTTTCTTCAGTTGTTTCACCCCTTAATTTGGCGCTAAGGGG‑ATCACCCCCCGCCGTTGCTCTTACTCgg                                                                                                                                                                                                      <  2:938626/136‑1 (MQ=17)                                                                                                                                    cTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGa                                                                                             >  1:1248402/1‑239 (MQ=255)                                                                                                                                      gACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCg                                                                                                                                                                                                       >  1:569264/1‑131 (MQ=255)                                                                                                                                      gACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAAGGGG‑ATCACTCCCCGCCGTTGCTCTTACTCg                                                                                                                                                                                                       <  2:569264/131‑1 (MQ=255)                                                                                                                                           ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGc                                                                                                                                                                                            >  1:938626/1‑137 (MQ=255)                                                                                                                                           ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGc                                                                                                                                                                                            >  1:1077540/1‑137 (MQ=255)                                                                                                                                           ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGa                                                                                       >  1:200074/1‑238 (MQ=255)                                                                                                                                           ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAAAGATGACGAt                                                                                      >  1:317063/1‑239 (MQ=255)                                                                                                                                             taaaaacCGAAATTGAAGACCGCGGGGCATCTTCTTAACGGGGCCCCGAGGCGCCCCATTCTTCAGATGTTTCACCCCTTAAGTTAGCGCTAATGGG‑ATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGtt                                                                                                                                                             <  2:524225/163‑1 (MQ=255)                                                                                                                                                aaaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATAACGATGACGATCAggg                                                                                 >  1:1031934/1‑239 (MQ=255)                                                                                                                                                       aatgaaagaaCGCGAGGCATGCTCTAAAGGGGGCCCCGGGGCGCCGCTTTCTTCAGATGGTTCACCCCTTAATTTAGCGCAAGGGGG‑ATCCCTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAAcc                                                                                                                                                 <  2:122526/163‑1 (MQ=255)                                                                                                                                                       aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGaa                                                                                                                                                   >  1:524225/1‑166 (MQ=255)                                                                                                                                                       aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACaa                                                                           >  1:114736/1‑238 (MQ=255)                                                                                                                                                                 cgcgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGt                                                                                                                                       >  1:122526/1‑168 (MQ=255)                                                                                                                                                                 cgcgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATa                                                                                                 >  1:83038/1‑206 (MQ=255)                                                                                                                                                                      ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGAAGATCAGGGCGACAATCATCGCCCTAt                                                               >  1:328695/1‑235 (MQ=255)                                                                                                                                                                                       gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAAGAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGGGACAATCATCGCCTTATGCTGCTTCATTGCTCtcttc                                           >  1:992609/1‑238 (MQ=255)                                                                                                                                                                                       gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCtt                                                                                                                                    <  1:920234/149‑1 (MQ=255)                                                                                                                                                                                       gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCtt                                                                                                                                    >  2:920234/1‑149 (MQ=255)                                                                                                                                                                                       gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCtcttc                                           >  1:966364/1‑238 (MQ=255)                                                                                                                                                                                       gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCtcttc                                           >  1:208144/1‑238 (MQ=255)                                                                                                                                                                                                 gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGc                                                                    >  1:968283/1‑203 (MQ=255)                                                                                                                                                                                                            cAGTTGGTTCACCCCTTAGGTTAGCCCAAGGGGG‑ATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGAt                                                                                                                                                        <  2:1004941/108‑1 (MQ=255)                                                                                                                                                                                                                    tCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGa                                                                                                                                                                                                     >  1:1058826/1‑55 (MQ=255)                                                                                                                                                                                                                    tCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCCTTGTTTTCTTCTCCTTGa                                    >  2:947270/1‑216 (MQ=255)                                                                                                                                                                                                                    tCAACCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGACACAGATGACGACCAGGGCGAcc                                                                            >  1:446992/1‑175 (MQ=255)                                                                                                                                                                                                                      aaCCCTTAAGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTca                                                                                                                                              >  1:1004941/1‑108 (MQ=255)                                                                                                                                                                                                                        cccTTAATTTACCGCTAAGGGG‑ACCACCCCCCGCCGTTGCTCTTACTCGGa                                                                                                                                                                                                     <  2:63053/51‑1 (MQ=255)                                                                                                                                                                                                                             aaGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGTATGTGAACCTCACAGAGGTCTTTTCTCTTTACCAGCGCCGCCACTACGGCGGTGATACAGATTACGATCATGGCCACAATCATCGCCTTTTGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTTCa        >  1:729309/1‑235 (MQ=255)                                                                                                                                                                                                                             aaGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCCCTTATCCTTGACCTTACCGTCAGTAAGAGGCACCCTACa        >  1:324139/1‑235 (MQ=255)                                                                                                                                                                                                                             aaGTTAGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGCGATACAGATGACGATCAGGGCGGCAATCCTCGCCTTATGCTGCTTCATTGCTCTCCTCTCCTTGTCCTTACGGTCAGTAAGAGGCACCTTACa        >  2:820924/1‑235 (MQ=255)                                                                                                                                                                                                                                ttAGCGCTAAAGGGCATCCCCCCCCGCCGTTGCTCTTACTCGGATGCGTAAGCCGTGAAAACACCAACCTCCGGCGGGCCAGTTCGGAGGGGAACCCCACAGAGGTTTTTGCTCGTTACCAGGGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGtt  <  1:820924/239‑1 (MQ=255)                                                                                                                                                                                                                                  aGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGcc                                                                                                                                                                                           >  1:63053/1‑51 (MQ=255)                                                                                                                                                                                                                                  aGCGCTAATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGAGATACAGATGACGATCAGGGCGACAAACATCGCCTTATGCTGCTTTATTGCTCTCTTCTCCTTGAACTTACGGGACGTAAGAGGGAATCTACAtgt     >  1:291825/1‑233 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       GACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGAC‑AGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑ATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTT  >  NC_000913/16495‑16954

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