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%	125.0 / 2.2	55	I358K (ATA→AAA)	yneK	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (32/13);  new base A (4/6);  total (36/19)
Fisher's exact test for biased strand distribution p-value = 7.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAACACTGTTTCCCGCAAAGTGGCGTGGCTACGG  >  NC_000913/1616604‑1617069                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           aGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGGGAAAAAAAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑AATACCCGCa                                                                                                                                                                                                                                     >  1:451951/1‑239 (MQ=255)     cTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCACTTTATGAATCCCAATGACAGAGATTTGAATGACGTGTTGAAGACAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCCAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGCAGACCTGTTGAAAACCAG‑GCCTCCCTCCCCTGTCAG‑A‑CCACA‑TATAcccaacaca                                                                                                                                                                                                                                  >  1:496194/1‑232 (MQ=255)        gTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTTTAAAAAATGAAAAACCGCCAAAAGTT‑ACTCACAAAGACATGGTGAAAAAAAG‑GCCTCCCCCCCCTGTCAG‑A‑ACACC‑AAAACCCGAAACAGcaa                                                                                                                                                                                                                              >  2:24045/1‑237 (MQ=255)                ttttCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGATCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCCCCAAAACGcc                                                                                                                                                                                                                       >  1:508813/1‑238 (MQ=255)                                   gccattacacaAATTTAAGAATCCAATTAAAAGTAATTTAAATGACGTTTTGAGGCCAGAATCATTAAGTGATAAAAAAATTAACAAACATTTGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATAAACCCCCGCCAAAAGTT‑CCTCACGAAGACCTGGTGAAAAACAG‑GCCCCGCCCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCAtt                                                                                                                                                                                                   <  2:816607/232‑1 (MQ=255)                                     aCTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCCTTaa                                                                                                                                                                                                 >  1:270157/1‑239 (MQ=255)                                      cTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCCCGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACCCCAGCCCCCCCCGAACGTCATTaa                                                                                                                                                                                                 >  2:23600/1‑238 (MQ=255)                                      cTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTAAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAACGTT‑ACGACCGAAGACCTGGTGAAAAACAG‑GACTCGCTCCCCTGTAAG‑A‑CCACC‑TAAACCAGCAACCGCCAAAACCCCAGCCCCCCCTACACGTCATaaa                                                                                                                                                                                                 >  2:664329/1‑238 (MQ=255)                                         tCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAAGGAAACAACCCAAAAAGTT‑ACTCACGAAAACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGGCAG‑A‑CCACA‑TATACCCCCAACCGACAAAACCCAAGACCCCCCTGAACCGCACTAAAccc                                                                                                                                                                                             >  2:1074246/1‑238 (MQ=255)                                         tCACCAATTTATGAATCCCAATGAAAGAGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACCCACGAAGACCTGTGGAAAAAAAG‑ACCTCGCACCCCCTCCAA‑A‑CCACC‑TAAAACCGCAACCACCAAAAAGCCAGACCCCCCTGAACgcactaaaaccg                                                                                                                                                                                             >  2:364002/1‑228 (MQ=255)                                                tttATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCCCCAAACAGCCAGACCTCCCTGAACGTCATTAAACCGTGATGtt                                                                                                                                                                                      >  1:1182930/1‑239 (MQ=255)                                                   aTGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGCCCTCCCTGAACGTCATTAAACCGTGATGTTAc                                                                                                                                                                                    >  1:463048/1‑238 (MQ=255)                                                   aTGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAATT‑ACTCACGAAGACCTGGTGAAAACCAG‑GCCTCGCTCCCCTGCCAG‑A‑CACCC‑AAAACCCGCAACCCCCAAAACCCCAAACCTCCTCAAACACCATTAAACCGtaccgttccc                                                                                                                                                                                   >  2:1198277/1‑230 (MQ=255)                                                                  aGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGt                                                                                                                                                                    >  1:904290/1‑239 (MQ=255)                                                                           aaTGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGCCATTAAAACCGAATGTTACAAACCCCCTGACCCGCGaaaaaaac                                                                                                                                                           >  2:561488/1‑237 (MQ=255)                                                                                     tGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCagag                                                                                                                                                 >  1:19869/1‑239 (MQ=255)                                                                                         gCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGGGAAAGAATCAGCGTCAGAGAAAc                                                                                                                                             >  1:790762/1‑239 (MQ=255)                                                                                                          atgaaAAAATAAAAAAACATAGGGC‑CGGCGAACCCCAAAACTTGTTTAAAAAAAAAGAAACCCCCCAAAATATA‑ACTCCCAAAGACCGGGGAAAAAACAG‑CCCCCGCCCCCCTGT‑AG‑A‑CCCCCTTAACCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAg                                                                                                                             <  2:412019/238‑1 (MQ=255)                                                                                                             atAAAATTATCAAACATTGCGGCCGGAGCACCACAACAC‑TGTTTGAAAAACATAAACCCCCGCCAAAATTA‑ACACACAAAGCCCT‑GTGAAAAACAGTGCCGCGC‑CCCCTGGCAA‑ACCCACC‑TAACCCCACAACCGCCAAAGCTCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATca                                                                                                                         <  2:1182930/239‑1 (MQ=255)                                                                                                                  aTTAACAAACATTGGGCCACCAGGACCCTCAAAC‑TGTTGGTAAAAAATGAAACCCCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAAt                                                                                                                    <  1:188360/239‑1 (MQ=255)                                                                                                                     acaaaaaATAGGGGGCGCAGGACCCTAAAAT‑GTTTTGAAAAAAAGGAACCCCCGCAAAAAGTT‑CCCCACGAAGCCCTGGTAAAAAACAG‑GCCTCGCTCCCCTGCAAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCa                                                                                                                 <  1:561488/236‑1 (MQ=255)                                                                                                                     acaaaaCATATGGCGCGCAGGCCCCCCAAAC‑GGTTTGAAAAAAATGAACCCCCGCAAAAAGTA‑CCTCACGAAGACCGGGTGAAAAACAG‑CCCTCCCTCCCCTGTAAG‑A‑CCCCC‑TATACCCACAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCa                                                                                                                 <  1:664329/236‑1 (MQ=255)                                                                                                                          aCATAGGGCGCGAGGGACCCAAAAAT‑GGTTTGAAAAAAATGAAACCCCGCAAAAAGTT‑ACACACGAAGACCGGGTGAAAACCCG‑GCCTCGCTCCCCTGTCAG‑A‑CCCCC‑TATCCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTa                                                                                                            <  1:100542/239‑1 (MQ=255)                                                                                                                             tatGGCGCGCAGGACCCTCAAAC‑TGTTTGGAAAAAATGAAACCCCGCCAAAAGTT‑ACTCACGAAGACCTGGTAAAAAACGG‑GCCCCGCTCCCCTGTCAG‑A‑CCACC‑TATCCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTc                                                                                                         <  1:861310/239‑1 (MQ=255)                                                                                                                             tatGGCGCGCAGGACCCTAAAAC‑TGTTTGAAAAAAATGAAACCCCGCCAAAAGTT‑ACTCACGAAGACCTGGTAAAAAACAG‑GCCCCGCTCCCCTGTCAG‑A‑CCACC‑TATCCCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTc                                                                                                         <  1:436793/239‑1 (MQ=255)                                                                                                                              atGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAAAGCGATCCAGATCACAAATGCATTGAATTc                                                                                                         >  1:892089/1‑238 (MQ=255)                                                                                                                                ggCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGAATTccca                                                                                                      >  1:344282/1‑239 (MQ=255)                                                                                                                                  cgcgcAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCCCAAATGCATTGCATTCCcatc                                                                                                    >  1:972675/1‑239 (MQ=255)                                                                                                                                       agaccaaaaaaaC‑TTTTTTAAAAAAAAAAAACCCCCCCAAAATTT‑ACCCCAGAACACCGGGACAAAAACAG‑GCCCCACTCCCCTGCCAA‑C‑CCCCC‑TAAACCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTaa                                                                                               <  2:122612/230‑1 (MQ=255)                                                                                                                                          ccccctaaaC‑TGTTTG‑AAAAAATTAAACCCCCCCAAAAATTGCCCCCCAAAGTCCGGGTAAAAAACAG‑CCCTGGCCCCCCTTTAAG‑A‑CCCCC‑TATCCCCACAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTaa                                                                                               <  1:327437/230‑1 (MQ=255)                                                                                                                                              tCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGAGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAAAGGa                                                                                                                                          >  1:754519/1‑189 (MQ=255)                                                                                                                                                   c‑TGTTTTAAAAAAATGAAACCCCGCAAAAAGTT‑CCCAACAAAGACCTGGTAAAAAACAG‑GCCCCGCTCCCCTGTCAG‑A‑CCACC‑TATCCCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGAt                                                                                   <  1:24045/239‑1 (MQ=255)                                                                                                                                                       tttGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCCCAAATGCATTGTATTCCCATCATTAACCGTTTTAAGATAAt                                                                                >  1:587595/1‑239 (MQ=255)                                                                                                                                                       tttGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTAtt                                                                                                          >  1:1106713/1‑213 (MQ=255)                                                                                                                                                       tttGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCAt                                                                                >  1:903395/1‑239 (MQ=255)                                                                                                                                                       tttGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCAAATCATTAACCGTTTTAAATTCAt                                                                                >  1:164242/1‑239 (MQ=255)                                                                                                                                                       tttGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGACTGAAAGAATCAGCGTCAGAGAAACGGCAAACCCGCTCCAGATCACAAAAGCAATTTATTAACATCATTAACCGTTTTAAAACCAt                                                                                >  2:535135/1‑239 (MQ=255)                                                                                                                                                         tGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACCTCATTAACCGTTTTAAGATAAtt                                                                               >  1:986810/1‑238 (MQ=255)                                                                                                                                                           tAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGactgtctct                                                                            >  1:122612/1‑231 (MQ=255)                                                                                                                                                                          cGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACACCGGAAAGAATCACCGTCAGAGAAACGGAAAACGCGATCACAACCACAAATGCATTGTATTCACATAATTAACCGTTTTAAGACAACTTACTCACttttttgaaa                                                              >  2:528062/1‑233 (MQ=255)                                                                                                                                                                           ccaaaatGTT‑ACTCACAAACACCTGTGAAAAAACAG‑GCCCCCCCCCCCGTCCAA‑C‑CCCCC‑TCAACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACt                                                            <  2:163156/236‑1 (MQ=255)                                                                                                                                                                               aaaGTC‑ACCCCCGAACACCGGGTAAAAAACAG‑GCCCCGCTCCCCTGTCAG‑A‑CCCCC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCAcc                                                        <  1:341686/239‑1 (MQ=255)                                                                                                                                                                                       cccacaaaacccTGGCAAAAACCAG‑GCC‑CCCCCCCCTGTCCGTC‑CCACC‑TAACCCCGCAACCGCCAAAACGCCAGCCCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGAGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAAt                                                 <  2:270157/229‑1 (MQ=255)                                                                                                                                                                                       cTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATAAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCCTTTCATCACTTTTTCGCAACTTACCCGATaa                                                  >  1:147658/1‑238 (MQ=255)                                                                                                                                                                                                 ccgggaaaaaaaCGG‑GCCTCCCTCCCCTGTCAG‑A‑CCCCC‑TATACCCGAAACCGGCAAAACGCCAGACCCCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGac                                       <  1:1198277/232‑1 (MQ=255)                                                                                                                                                                                                 ccTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGaa                                                                                                                                                                                                          <  1:836043/76‑1 (MQ=255)                                                                                                                                                                                                 ccTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGaa                                                                                                                                                                                                          >  2:836043/1‑76 (MQ=255)                                                                                                                                                                                                    tggaaaaacCAG‑GCCCCCCCCCCCGCTCAA‑C‑CCCCC‑TTTACCCCTAACCGCCCAAACGCCCGACCTGCCGGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGacaac                                    <  2:986810/235‑1 (MQ=255)                                                                                                                                                                                                    ggataaacaCGG‑GCCCCCCTCCCCTGTCAG‑G‑CCCCC‑TAACCCCGCAACCGCCAAAACGCCAGGCCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGacaac                                    <  1:1186033/235‑1 (MQ=255)                                                                                                                                                                                                      tGAAAAACAG‑GCCTCGCTCCACTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTCCCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACCCGCTCCAGATCACAAATGCATTGAATTCCAATAATTAACCGTTTTAAGTTCATTTCATCCTTTTTTGCCAACTAACCCAATAATTTATTATAAcaaaaa                                  >  2:1138051/1‑239 (MQ=255)                                                                                                                                                                                                          aaaCAG‑GCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACaacaaac                                >  1:150292/1‑237 (MQ=255)                                                                                                                                                                                                               g‑gCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCCCATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAACACTGtt                          >  1:1180305/1‑238 (MQ=255)                                                                                                                                                                                                               g‑gCCTCGCTCCCCTGTCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAg                                                                                     >  1:541227/1‑179 (MQ=255)                                                                                                                                                                                                                               tCAG‑A‑CCACC‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCCCTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAAAACTGTTTCCCCCAAAatggcgt          >  1:501668/1‑239 (MQ=255)                                                                                                                                                                                                                                         cc‑TATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAATCTGTTATGACAACAAACACTGTTTCCCGCAAAGTGGCGTGGCTACgg  >  1:359365/1‑239 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           AGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAAC‑TGTTTGTAAAAAATGAAACACCGCCAAAAGTT‑ACTCACGAAGACCTGGTGAAAAACAG‑GCCTCGCTCCCCTGTCAG‑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 < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑