Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,883	1	.	C	93.6%	143.7 / 2.8	48	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (26/18);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 8.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815648‑3816132                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                            aaaaCAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACCCCTGCTGTCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGTTTTCATGCCTTCGCTCCTCATCTTACTTTCCTACAGACAAAAAAAAGGCGACTCACCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGcct                                                                                                                                                                                                                                               <  2:413354/250‑2 (MQ=255)      aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACGCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCcaaca                                                                                                                                                                                                                                            >  1:54795/1‑249 (MQ=255)          cGGCGTTGAAGAAATAGTGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCACCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTGCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGa                                                                                                                                                                                                                                      <  1:234549/251‑1 (MQ=255)            gCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGt                                                                                                                                                                                                                                    <  2:240898/251‑1 (MQ=255)                      aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCGTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                                                                                                                                                                                                                          <  1:201172/251‑1 (MQ=255)                       aTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                                                                                                                                                                                                         >  1:143284/1‑251 (MQ=255)                       aTAGGCGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTCAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCACCAAGATGAGTAGCTCTTCATg                                                                                                                                                                                                                         >  1:84003/1‑251 (MQ=255)                            ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggggg                                                                                                                                                                                                                     >  1:251563/1‑248 (MQ=255)                             gCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGaa                                                                                                                                                                                                                   <  2:145333/251‑1 (MQ=255)                               ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTATCTCTTCATggggtaacg                                                                                                                                                                                                                 >  1:381668/1‑245 (MQ=255)                                   gtgcCCGGATTTCAGCGTAAACCCGCCAAACTTTAACACCTTCTTGCTAAGCGCAAATTCAATAAACGGGCGCTG‑ATATGGGTTCATGCCTTCGCTCCTCATCTTACTGTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                                             <  2:1130/249‑1 (MQ=255)                                   gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                                             <  1:228665/251‑1 (MQ=255)                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACATGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACACCACAGTCGCCTTAAAAATCAGTTTGACCGCACCGCCTGCTGCGTCGCTACAATGACTTCGAT‑TCC‑CCATCGGGCCAGAGCCAACAAGCTGAAGAGAACTTCATGGAGGTACGGCTCGCCtt                                                                                                                                                                                                        >  2:253815/1‑251 (MQ=255)                                               cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCACGCCTTCTGCCGt                                                                                                                                                                                                 >  1:279869/1‑251 (MQ=255)                                                       aCTCGCCAAACTTTAACTCCTGCTTTCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa                                                                                                                                                                                         <  1:337750/251‑1 (MQ=255)                                                                           tgcttgctAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCt                                                                                                                                                                     >  2:93311/1‑251 (MQ=255)                                                                              ttgctAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACACGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCg                                                                                                                                                                  >  2:464236/1‑251 (MQ=255)                                                                               tgctACGCGCAAAGTCAATAAACTGGCG‑TGTATATTTTTTCATGCATTTGTCCCTACTCTTCTTTTTCTACAGACACAAAAAAGGCGCTTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGGCGCTACAATTGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTGCTGCCGTACCCTGCACTTCAATGATGCGCCAGTCTTc                                                                                                                                                                   <  2:355917/249‑1 (MQ=255)                                                                               tgctAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTc                                                                                                                                                                   <  2:445278/249‑1 (MQ=255)                                                                                     gcgcAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatcc                                                                                                                                                           >  2:363587/1‑250 (MQ=255)                                                                                         aaaTTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCGCGTCTTCGGTCATCACTAc                                                                                                                                                       >  1:165324/1‑251 (MQ=255)                                                                                                    aCTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCgg                                                                                                                                           >  1:11686/1‑251 (MQ=255)                                                                                                             g‑ATAAGGTTTCATGCCTCCGCTCCGCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCCGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTACA‑TGTCGGTCTCTGCg                                                                                                                                   <  2:356532/250‑1 (MQ=255)                                                                                                                  tgggttCATGCCTTCGCCCCTCATCTTCCTTTTCCTCAGACAAAAAAAAGGTGACTCCTCAGTCGACTTAAAATTCAGTTTGCTAGCGCCTCCTTCTGCGTTGCGACAATTGATTCGCT‑TCCACC‑ACTGGCCAGAGCCAACAAGATGAGTCTCTCTTCATGGGTGAACGGCTCGCCGTTTGCCGTCCCCTGCTCGTCAATGATGCGACCTTCTTCTGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCag                                                                                                                               <  2:343305/246‑1 (MQ=255)                                                                                                                  ggttttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCATTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGATCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTTTGCGGCagt                                                                                                                              <  2:339475/248‑2 (MQ=255)                                                                                                                              ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCACCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑Gt                                                                                                                  >  2:171801/1‑250 (MQ=255)                                                                                                                                  ctcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTcc                                                                                                             >  1:181062/1‑251 (MQ=255)                                                                                                                                     ctcATCGTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCACCAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTCCGTTCA‑TGTCGGTCTATGCGGCAGAGTCTT‑AACAGTATTcc                                                                                                             >  2:260831/1‑248 (MQ=255)                                                                                                                                       cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATc                                                                                                        <  2:251563/251‑1 (MQ=255)                                                                                                                                            tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATTGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCa                                                                                                     <  1:306917/249‑1 (MQ=255)                                                                                                                                                tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATGCCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCg                                                                                                <  2:23140/250‑1 (MQ=255)                                                                                                                                                       aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCg                                                                                        >  1:323768/1‑251 (MQ=255)                                                                                                                                                          cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTc                                                                                     <  1:73570/251‑1 (MQ=255)                                                                                                                                                          aaaaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTc                                                                                     <  1:336319/250‑1 (MQ=255)                                                                                                                                                          aaaaaaaaaGGCGAATCAACAGTCGCCTTAAAAATCAGCTTGCCAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTc                                                                                     <  2:44233/250‑1 (MQ=255)                                                                                                                                                               aaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAt                                                                                >  1:199217/1‑251 (MQ=255)                                                                                                                                                                aaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCCAC‑GTATTCCAGATCGCCAACCGCTTCGCCGTTccccatt                                                                               >  2:205197/1‑245 (MQ=255)                                                                                                                                                                     cGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTCCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTCCGTACA‑TGTCGGCCTCTGCGGCAGAGTCTTCAAC‑GAAATCCAGATCGAAAACCGCTTCGCCGttacccattcctcc                                                                          >  2:2551/1‑239 (MQ=255)                                                                                                                                                                        cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                          >  1:346534/1‑248 (MQ=255)                                                                                                                                                                        cTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                                          >  2:73706/1‑248 (MQ=255)                                                                                                                                                                         tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGa                                                                       >  2:72390/1‑250 (MQ=255)                                                                                                                                                                          catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TTTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAAtttcgccagca                                                                      >  2:225179/1‑241 (MQ=255)                                                                                                                                                                            tcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACt                                                                   >  2:378295/1‑251 (MQ=255)                                                                                                                                                                                tCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCCAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACACTTCCGACAGAAACTGCgg                                                               >  2:127833/1‑251 (MQ=255)                                                                                                                                                                                        aaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                                     >  1:83732/1‑127 (MQ=255)                                                                                                                                                                                        aaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc                                                                                                                                                                                     <  2:83732/127‑1 (MQ=255)                                                                                                                                                                                                        cgcctccttctccggcggtaaaagggATTCGATGCCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTgg                                       <  2:317084/227‑1 (MQ=255)                                                                                                                                                                                                                    cGTCGCTACAATGGATTCGAT‑TCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGCAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGc                           >  2:388224/1‑251 (MQ=255)                                                                                                                                                                                                                                     ggagtccc‑cc‑GCGGGCCAGAGCCAGCAAGTTGAGTGGCTTTCCATGGTGGAACGTCTCGCCTGCTTCTGTCCCCTGCATTTCAATGTTGCGCCCGTCTTCGGTCATCGCT‑GGTGGAGTGTCGGGCTCTGCGTCAGAGTCTCCAAC‑GTATCCCAGATCGCAAAACTCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACTATCCCGTTCATCGGAGTGGTTTTCAGCTTGCTGTTTTCCGCCAGCttt          <  2:37428/245‑2 (MQ=255)                                                                                                                                                                                                                                       attccc‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTCCCATCCCTTTCATCGGTTTAGTTTTCAGCTTGCCGTTTTCCACCAGCTTCtt        >  2:111236/5‑250 (MQ=255)                                                                                                                                                                                                                                          ccc‑cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt       >  1:169128/2‑249 (MQ=255)                                                                                                                                                                                                                                           cc‑cc‑TCGGGCCATAGCCAAAAAGTTGAGGAGCCCGGGATGGTTGAACGGCGCGCCTTCTGCCTTCCCCTGCACTTCAATGATGCGCCCTTCTGCGGTCAGCCCTTCGGTCA‑TGTCGGTCTCTGCTGCAGAGTCTTCAAC‑GTATTCCAGATCGCAACCCGCTTCGCCGGTCACAAGTCCGCCAGAAACTGCTGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGCTTTCCACCAGCTTCTGTAgc    <  2:381668/251‑1 (MQ=255)                                                                                                                                                                                                                                              cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcgc  >  1:468631/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                            TAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTG‑ATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCA‑TGTCGGTCTCTGCGGCAGAGTCTTCAAC‑GTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815648‑3816132

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