Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,616,835	T→A	16.5%	I358K (ATA→AAA)	yneK →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,835	0	T	A	16.5%	128.0 / 2.4	60	I358K (ATA→AAA)	yneK	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (33/16);  new base A (5/6);  total (38/22)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

AGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAACACTGTTTCCCGCAAAGTGGCGTGGCTACGG  >  NC_000913/1616604‑1617069                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           aGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGGGAAAAAAAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑AATACCCGCa                                                                                                                                                                                                                                     >  2:451951/1‑239 (MQ=255)     cTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCACTTTATGAATCCCAATGACAGAGATTTGAATGACGTGTTGAAGACAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCCAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGCAGACCT‑GTTGAAAACCAG‑GCCTCCCTCCCCTGT‑CAG‑A‑CCACA‑TATAcccaacaca                                                                                                                                                                                                                                  >  2:496194/1‑232 (MQ=255)        gTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTTTAAAAAATGAAAAACCGCCAAAAGTT‑AC‑TCACAAAGACAT‑GGTGAAAAAAAG‑GCCTCCCCCCCCTGT‑CAG‑A‑ACACC‑AAAACCCGAAACAGcaa                                                                                                                                                                                                                              >  1:24045/1‑237 (MQ=255)                ttttCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGATCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCCCCAAAACGcc                                                                                                                                                                                                                       >  2:508813/1‑238 (MQ=255)                                   gccattacacaAATTTAAGAATCCAATTAAAAGTAATTTAAATGACGTTTTGAGGCCAGAATCATTAAGTGATAAAAAAATTAACAAACATTTGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATAAACCCCCGCCAAAAGTT‑CC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCCCGCCCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCAtt                                                                                                                                                                                                   <  1:816607/232‑1 (MQ=255)                                     aCTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCCTTaa                                                                                                                                                                                                 >  2:270157/1‑239 (MQ=255)                                      cTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCCCGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACCCCAGCCCCCCCCGAACGTCATTaa                                                                                                                                                                                                 >  1:23600/1‑238 (MQ=255)                                      cTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTAAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAACGTT‑AC‑GACCGAAGACCT‑GGTGAAAAACAG‑GACTCGCTCCCCTGT‑AAG‑A‑CCACC‑TAAACCAGCAACCGCCAAAACCCCAGCCCCCCCTACACGTCATaaa                                                                                                                                                                                                 >  1:664329/1‑238 (MQ=255)                                        aTCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAACCATATGGCGCGCAGGA‑CCCTCCAAC‑TGTTTGTAAAAAATGAAAAACGCCCACAAGTT‑AT‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCCCGCACCCCGGT‑CTG‑A‑CCACC‑TAAACCCTCCACCGCCAAAAGCCCACAccccccctaacaccctaaacac                                                                                                                                                                                              >  1:1111924/1‑223 (MQ=255)                                         tCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAAGGAAACAACCCAAAAAGTT‑AC‑TCACGAAAACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGG‑CAG‑A‑CCACA‑TATACCCCCAACCGACAAAACCCAAGACCCCCCTGAACCGCACTAAAccc                                                                                                                                                                                             >  1:1074246/1‑238 (MQ=255)                                         tCACCAATTTATGAATCCCAATGAAAGAGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑CCACGAAGACCT‑GTGGAAAAAAAG‑ACCTCGCACCCCCTC‑CAA‑A‑CCACC‑TAAAACCGCAACCACCAAAAAGCCAGACCCCCCTGAACgcactaaaaccg                                                                                                                                                                                             >  1:364002/1‑228 (MQ=255)                                                tttATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCCCCAAACAGCCAGACCTCCCTGAACGTCATTAAACCGTGATGtt                                                                                                                                                                                      >  2:1182930/1‑239 (MQ=255)                                                   aTGAATCCCAATGAAAGTGATTTGAATGCCGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGGAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGCCCT‑GGTGAAAACCAG‑GCCTCGCGCCCCTTT‑CAT‑A‑CCACC‑TATACCCGCAACCGCCAAAACCCCAAAACTCCCTGAACGCCATAAAAcagatatgtgaac                                                                                                                                                                                   >  1:551985/1‑227 (MQ=255)                                                   aTGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGCCCTCCCTGAACGTCATTAAACCGTGATGTTAc                                                                                                                                                                                    >  2:463048/1‑238 (MQ=255)                                                   aTGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAATT‑AC‑TCACGAAGACCT‑GGTGAAAACCAG‑GCCTCGCTCCCCTGC‑CAG‑A‑CACCC‑AAAACCCGCAACCCCCAAAACCCCAAACCTCCTCAAACACCATTAAACCGtaccgttccc                                                                                                                                                                                   >  1:1198277/1‑230 (MQ=255)                                                                  aGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGt                                                                                                                                                                    >  2:904290/1‑239 (MQ=255)                                                                           aaTGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGCCATTAAAACCGAATGTTACAAACCCCCTGACCCGCGaaaaaaac                                                                                                                                                           >  1:561488/1‑237 (MQ=255)                                                                                     tGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCagag                                                                                                                                                 >  2:19869/1‑239 (MQ=255)                                                                                         gCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGGGAAAGAATCAGCGTCAGAGAAAc                                                                                                                                             >  2:790762/1‑239 (MQ=255)                                                                                                          atgaaAAAATAAAAAAACATAGGGC‑CGGCGAA‑CCCCAAAACTTGTTTAAAAAAAAAGAAACCCCCCAAAATATA‑AC‑TCCCAAAGACCG‑GGGAAAAAACAG‑CCCCCGCCCCCCTGT‑‑AG‑A‑CCCCCTTAACCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAg                                                                                                                             <  1:412019/238‑1 (MQ=255)                                                                                                             atAAAATTATCAAACATTGCGGCCGGAGCA‑CCACAACAC‑TGTTTGAAAAACATAAACCCCCGCCAAAATTA‑AC‑ACACAAAGCCCT‑‑GTGAAAAACAGTGCCGCGC‑CCCCTGG‑CAA‑ACCCACC‑TAACCCCACAACCGCCAAAGCTCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATca                                                                                                                         <  1:1182930/239‑1 (MQ=255)                                                                                                                  aTTAACAAACATTGGGCCACCAGGA‑CCCTCAAAC‑TGTTGGTAAAAAATGAAACCCCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAAt                                                                                                                    <  2:188360/239‑1 (MQ=255)                                                                                                                   ttttccaaataatgggcggcAGGACCCCTAAAAG‑TGTTGGAAAAAAATGAACCCCCCCCAAAGGTT‑AC‑‑CACAAAAGCCTGGGTAAAAACCAG‑GCCTCGCTCCCCGGT‑CAG‑A‑CCTCC‑TATCCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAAt                                                                                                                    <  1:112623/221‑1 (MQ=255)                                                                                                                     acaaaaaATAGGGGGCGCAGGA‑CCCTAAAAT‑GTTTTGAAAAAAAGGAACCCCCGCAAAAAGTT‑CC‑CCACGAAGCCCT‑GGTAAAAAACAG‑GCCTCGCTCCCCTGC‑AAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCa                                                                                                                 <  2:561488/236‑1 (MQ=255)                                                                                                                     acaaaaCATATGGCGCGCAGGC‑CCCCCAAAC‑GGTTTGAAAAAAATGAACCCCCGCAAAAAGTA‑CC‑TCACGAAGACCG‑GGTGAAAAACAG‑CCCTCCCTCCCCTGT‑AAG‑A‑CCCCC‑TATACCCACAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCa                                                                                                                 <  2:664329/236‑1 (MQ=255)                                                                                                                       aaaacattgtgggcagcgca‑CCCACAAAA‑TGTTTTTAAAAAAAGAAACCCCCACAAAAATT‑CCACCAAGAAACCCT‑GGTAAAAACCAG‑CCCTCGC‑CCCCTCTGCAA‑C‑CCACC‑TATCCCCCCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCAt                                                                                                                <  1:1322381/220‑1 (MQ=255)                                                                                                                          aCATAGGGCGCGAGGGA‑CCCAAAAAT‑GGTTTGAAAAAAATGAAACCCCGCAAAAAGTT‑AC‑ACACGAAGACCG‑GGTGAAAACCCG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCCCC‑TATCCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTa                                                                                                            <  2:100542/239‑1 (MQ=255)                                                                                                                             tatGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGGAAAAAATGAAACCCCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTAAAAAACGG‑GCCCCGCTCCCCTGT‑CAG‑A‑CCACC‑TATCCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTc                                                                                                         <  2:861310/239‑1 (MQ=255)                                                                                                                             tatGGCGCGCAGGA‑CCCTAAAAC‑TGTTTGAAAAAAATGAAACCCCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTAAAAAACAG‑GCCCCGCTCCCCTGT‑CAG‑A‑CCACC‑TATCCCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTc                                                                                                         <  2:436793/239‑1 (MQ=255)                                                                                                                              atGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAAAGCGATCCAGATCACAAATGCATTGAATTc                                                                                                         >  2:892089/1‑238 (MQ=255)                                                                                                                                ggCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGAATTccca                                                                                                      >  2:344282/1‑239 (MQ=255)                                                                                                                                  cgcgcAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCCCAAATGCATTGCATTCCcatc                                                                                                    >  2:972675/1‑239 (MQ=255)                                                                                                                                        agaccaaaaaaaC‑TTTTTTAAAAAAAAAAAACCCCCCCAAAATTT‑AC‑CCCAGAACACCG‑GGACAAAAACAG‑GCCCCACTCCCCTGC‑CAA‑C‑CCCCC‑TAAACCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTaa                                                                                               <  1:122612/230‑1 (MQ=255)                                                                                                                                           ccccctaaaC‑TGTTTG‑AAAAAATTAAACCCCCCCAAAAATTGCC‑CCCCAAAGTCCG‑GGTAAAAAACAG‑CCCTGGCCCCCCTTT‑AAG‑A‑CCCCC‑TATCCCCACAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTaa                                                                                               <  2:327437/230‑1 (MQ=255)                                                                                                                                               tCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGAGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAAAGGa                                                                                                                                          >  2:754519/1‑189 (MQ=255)                                                                                                                                                    c‑TGTTTTAAAAAAATGAAACCCCGCAAAAAGTT‑CC‑CAACAAAGACCT‑GGTAAAAAACAG‑GCCCCGCTCCCCTGT‑CAG‑A‑CCACC‑TATCCCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGAt                                                                                   <  2:24045/239‑1 (MQ=255)                                                                                                                                                        tttGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCCCAAATGCATTGTATTCCCATCATTAACCGTTTTAAGATAAt                                                                                >  2:587595/1‑239 (MQ=255)                                                                                                                                                        tttGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTAtt                                                                                                          >  2:1106713/1‑213 (MQ=255)                                                                                                                                                        tttGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCAt                                                                                >  2:903395/1‑239 (MQ=255)                                                                                                                                                        tttGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCAAATCATTAACCGTTTTAAATTCAt                                                                                >  2:164242/1‑239 (MQ=255)                                                                                                                                                        tttGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGACTGAAAGAATCAGCGTCAGAGAAACGGCAAACCCGCTCCAGATCACAAAAGCAATTTATTAACATCATTAACCGTTTTAAAACCAt                                                                                >  1:535135/1‑239 (MQ=255)                                                                                                                                                          tGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACCTCATTAACCGTTTTAAGATAAtt                                                                               >  2:986810/1‑238 (MQ=255)                                                                                                                                                            tAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGactgtctct                                                                            >  2:122612/1‑231 (MQ=255)                                                                                                                                                                           cGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACACCGGAAAGAATCACCGTCAGAGAAACGGAAAACGCGATCACAACCACAAATGCATTGTATTCACATAATTAACCGTTTTAAGACAACTTACTCACttttttgaaa                                                              >  1:528062/1‑233 (MQ=255)                                                                                                                                                                            ccaaaatGTT‑AC‑TCACAAACACCT‑GTGAAAAAACAG‑GCCCCCCCCCCCGTC‑CAA‑C‑CCCCC‑TCAACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACt                                                            <  1:163156/236‑1 (MQ=255)                                                                                                                                                                                aaaGTC‑AC‑CCCCGAACACCG‑GGTAAAAAACAG‑GCCCCGCTCCCCTGT‑CAG‑A‑CCCCC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCAcc                                                        <  2:341686/239‑1 (MQ=255)                                                                                                                                                                                        c‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATAAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCCTTTCATCACTTTTTCGCAACTTACCCGATaa                                                  >  2:147658/1‑238 (MQ=255)                                                                                                                                                                                         cccacaaaacccT‑GGCAAAAACCAG‑GCC‑CCCCCCCCTGT‑CCGTC‑CCACC‑TAACCCCGCAACCGCCAAAACGCCAGCCCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGAGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAAt                                                 <  1:270157/229‑1 (MQ=255)                                                                                                                                                                                               aaacccccggaaaaaaaCAG‑CCCCCCCCCCCCGGA‑CCG‑A‑CCCCC‑TATCCCCGCAACCGCCAAAACTCCAGGCCTCCCGGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGtt                                            <  1:904290/227‑1 (MQ=255)                                                                                                                                                                                                   ccT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGaa                                                                                                                                                                                                          >  1:836043/1‑76 (MQ=255)                                                                                                                                                                                                   ccT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGaa                                                                                                                                                                                                          <  2:836043/76‑1 (MQ=255)                                                                                                                                                                                                    ccgggaaaaaaaCGG‑GCCTCCCTCCCCTGT‑CAG‑A‑CCCCC‑TATACCCGAAACCGGCAAAACGCCAGACCCCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGac                                       <  2:1198277/232‑1 (MQ=255)                                                                                                                                                                                                       tggaaaaacCAG‑GCCCCCCCCCCCGCT‑CAA‑C‑CCCCC‑TTTACCCCTAACCGCCCAAACGCCCGACCTGCCGGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGacaac                                    <  1:986810/235‑1 (MQ=255)                                                                                                                                                                                                       ggataaacaCGG‑GCCCCCCTCCCCTGT‑CAG‑G‑CCCCC‑TAACCCCGCAACCGCCAAAACGCCAGGCCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGacaac                                    <  2:1186033/235‑1 (MQ=255)                                                                                                                                                                                                         tGAAAAACAG‑GCCTCGCTCCACTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTCCCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACCCGCTCCAGATCACAAATGCATTGAATTCCAATAATTAACCGTTTTAAGTTCATTTCATCCTTTTTTGCCAACTAACCCAATAATTTATTATAAcaaaaa                                  >  1:1138051/1‑239 (MQ=255)                                                                                                                                                                                                             aaaCAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACaacaaac                                >  2:150292/1‑237 (MQ=255)                                                                                                                                                                                                                  g‑gCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCCCATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAACACTGtt                          >  2:1180305/1‑238 (MQ=255)                                                                                                                                                                                                                  g‑gCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAg                                                                                     >  2:541227/1‑179 (MQ=255)                                                                                                                                                                                                                                  t‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCCCTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAAAACTGTTTCCCCCAAAatggcgt          >  2:501668/1‑239 (MQ=255)                                                                                                                                                                                                                                             cc‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAACACTGTTTCCCGCAAAGTGGCGTGGCTACgg  >  2:359365/1‑239 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           AGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGA‑CCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑AC‑TCACGAAGACCT‑GGTGAAAAACAG‑GCCTCGCTCCCCTGT‑CAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAACACTGTTTCCCGCAAAGTGGCGTGGCTACGG  >  NC_000913/1616604‑1617069

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