Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,616,800	G→A	23.4%	V346V (GTG→GTA)	yneK →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,800	0	G	A	23.4%	88.4 / 5.3	49	V346V (GTG→GTA)	yneK	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (27/10);  new base A (3/9);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 5.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTAGCAGAACAAGAATCTCCTCATTTTCAGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAA  >  NC_000913/1616574‑1617021                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              ttttAGCAGAACAAGAATCTCCTCATTTTCAGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACAGCC‑AAAAGTTACTCAC‑GAAGACCTGGTg                                                                                                                                                                                                                               >  2:38085/1‑227 (MQ=255)          aaCAAGAATCTCCTCATTTTCAGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGACAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTTAAACTGGTTGTAAAAAATGAAAAACAGCC‑AAAAGTTACCCAC‑CAAAACCTGGTGaaaaaaa                                                                                                                                                                                                                        >  1:565925/1‑223 (MQ=255)                atttttctctTTTTCAGAACTTTTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTAGGAATCCAAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCt                                                                                                                                                                                                               <  1:122631/222‑1 (MQ=255)                   ctcctcATTTTCAGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTccc                                                                                                                                                                                                            <  2:280392/228‑1 (MQ=255)                       tcATTTTCAGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGt                                                                                                                                                                                                        <  2:154795/228‑1 (MQ=255)                       tcATTTTCAGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAAAACCGCA‑AAAAGTTACTCAC‑GAAAACCTGGTAAAAAAAAAGCCTCAccccccag                                                                                                                                                                                                         >  1:280392/1‑225 (MQ=255)                             taaGAATTTGTTTTTTTTTTTCGCAAATATTGTATTTCACTTGCACCAATTTAAGAACCCAAATGAAAGTGATTTGAATGACGTGTTGAGGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAAGGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGacc                                                                                                                                                                                                  <  1:38085/226‑1 (MQ=255)                                gAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAAAACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAACCAGTCCCAGCTCCCCTTTCAGaccgc                                                                                                                                                                                                >  2:623933/1‑225 (MQ=255)                                     tGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAAAACCGCC‑ACAAGTTACTCAC‑GACGACCTGGTGAAAAACAGAACTCGCTCCCCGTTCAGacaccataaac                                                                                                                                                                                          >  1:547681/1‑219 (MQ=255)                                       tttttATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAAAACCGCC‑AAAATTTACTCAC‑GAAGCCCTGGTGAAAAACAGGCCTCGACACCCTGTCAGACCACCTATAacc                                                                                                                                                                                        >  2:248603/1‑228 (MQ=255)                                                               ttGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCCCCCCTGTCAACCCCCCTAAACCCGCAAACCCCAAAACCCCAAAcccc                                                                                                                                                                >  1:453189/1‑226 (MQ=255)                                                                 gCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACACCCTATACCCGCAACCGCCAAAACGCCAGACCTccc                                                                                                                                                              >  2:523307/1‑228 (MQ=255)                                                                  cACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCATACCTCCCt                                                                                                                                                             >  2:480162/1‑228 (MQ=255)                                                                    cTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGCCCACCTATACCCGCAACCGCCAAAACGCCAAACCTCCCTg                                                                                                                                                            >  2:464212/1‑227 (MQ=255)                                                                    cTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAACTTACTCAC‑GAAGACCTGGTAAAAAACAGGCCTCGCTCCCCTGTCAAACCACCTATACCCGCAACCCCAAAAACCCCAAAccccccca                                                                                                                                                            >  1:337619/1‑225 (MQ=255)                                                                     tttCACCAATTTAAGAACCCAAAAGAAAGTAATTTAAATAACGTGTGGAGGCCAGCATCATTAAGGGATAATAAAATTATAAAAAATAGGGCGCGCAGGCCCCAAAAACTGTTTGAAAAAAATGAAACCCCGCC‑AAAAGTTACCCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGaa                                                                                                                                                          <  1:40520/226‑1 (MQ=255)                                                                     tATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTcccc                                                                                                                                                                                                           <  2:73006/179‑1 (MQ=255)                                                                      aTCACCAATTTATGAATCCCTATGAAAGTGATTTGAATGACGTGTTGATGCCAGAATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGAAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCG‑AAAAGTTACTTAC‑GAAGACCTGGTGAAAAACAGGCCTCGTTCCCCTGTATGACCACCTAAACCCGCAACCGCCACAACGCCAGACCTCCCTGAAc                                                                                                                                                         >  2:549653/1‑228 (MQ=255)                                                                       tCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACg                                                                                                                                                        >  2:296351/1‑228 (MQ=255)                                                                       tCACCAATTTATGAATCCCAATGAAAGGGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACAACCTATACCCGCAACCGCCAAAAAGCCAGACTTCCCTGAACg                                                                                                                                                        >  2:118596/1‑228 (MQ=255)                                                                                           aagAAAGTTATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGAAAAAATTACCAACCATATGGCGCGAAGGCCCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATCCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTAc                                                                                                                                    <  2:453189/226‑1 (MQ=255)                                                                                            tGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAAGCCTCGCTCCCCTGTCAGACCACCAATACCCGCAACCGACAAAACGCCAGACCCCCCTGAACGTCATTAAACCGCGATGTTAcc                                                                                                                                   >  1:697654/1‑228 (MQ=255)                                                                                                aGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGTCGCGCAGGACCCTCAAACTATTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTAAC‑GAAGACCGTGTAAAAAACATGCCTCGCTCCCCGGTCAGACCACCTATACCCGCAACCGCCAAAACACCAGACCCCCCTGAACGACATTAAAACGAAATGTTACCgccc                                                                                                                               >  2:257710/1‑225 (MQ=255)                                                                                                    aTTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGAGGTTACCGActctct                                                                                                                           >  2:13422/1‑228 (MQ=255)                                                                                                                  ttGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGACCTCGCTCCCCTGTCAGACCACCTATACCCGCAACAGCCAAAAAGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACACCCCCTCTAACGCGCGAAAGaa                                                                                                             >  2:311559/1‑228 (MQ=255)                                                                                                                            catcatTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGa                                                                                                                                                                    >  2:444421/1‑163 (MQ=255)                                                                                                                                         tgatAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCCTGATGTTACCGACACTCTGACCCCTGAAAGAATAAGAGCCCGAGAACCGGaaaa                                                                                       >  2:271461/1‑227 (MQ=255)                                                                                                                                             aaaaTAATCAAACCTTGGCCACGCGAGACCCTCAAAATGTTGTAAAAAAAGGACCCACCCCC‑AAAAGTTACCCAC‑GAAACCCTGCTAAAAAACACGCCCCGCCCCCCTGTCACACCCCCTATACCCGCAACCGCCAAAACGCCAGGCCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGa                                                                                  <  1:686192/228‑1 (MQ=255)                                                                                                                                                  tataaaaaattgggCCCGCAGCACCACAAACCTGTTTGAAAAAAATGAAACCCCGCA‑AAAAGTT‑CACACAGAAGACCTGGTAAAAAACAGGCCTCCCTCCCCTGTCCGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAg                                                                             <  2:337619/216‑1 (MQ=255)                                                                                                                                                         catcgcccccaaaaccccccCAAATGGTTTGAAAAAAAACACACCCCCCC‑AAAACGACCACCC‑GAAGACCTGGTAAAAAACAGCCCTCCCCCCCCAGCCAAACCCCATATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGCCATTAAACCGCGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACaa                                                                      <  1:296351/212‑1 (MQ=255)                                                                                                                                                              ggCGCGCAGGACCCTCAAACTGTTTGTACAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCCCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACCCCAGACCTCCCTGAACGCCATTAAACCGTGATCTTCCCGCCTCTCTGACGCGTTACAGAATCAGCGACGGAGAAACGGAAAACCCCGTCCAACTCACaaacaca                                                                 >  1:19101/1‑224 (MQ=255)                                                                                                                                                              ggCGCGAAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACCCAC‑GAAGAACTGGGGAAAAAAAGGCCTCGACCCCCTGTCCGAACACCTCCCCCCGAAAACACCCCCCCGCCAGCCCTCCCTGAACGTCATTAAACACCGATGCTCCCGACCCTCTCACGCGGAAAACAACCCACCTCAGAAAAAAGCAAAACGAGAACCAGaaaccaaaacca                                                                 >  1:551222/1‑217 (MQ=255)                                                                                                                                                              gccgcgcATGCCCCTCAAACTG‑TTGTAAAAAACCAAACCCCCCCAAAAAGTACCCCAC‑GAACCCCGGGAGAAAAACCGCCCGCGCCCCCCTGCCACACCTCCTACACCCGCAGCCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCa                                                                 <  1:480746/226‑1 (MQ=255)                                                                                                                                                               gcgcgaggcaccccaaaaCTGTTAGAAAAAAATAAAACCCCCCC‑AAAAGTTACTCCA‑CAAGACCGGGTAAAAAACAGGCCTCGCTCCCCTGTAGGACCCCCTAACCCCGAAGCCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCAt                                                                <  1:547654/213‑1 (MQ=255)                                                                                                                                                                  cgcAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGACAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACCCGGGCAAGAATCAGCGCCAGAGAACCGGAAAACGCGATCCAGATCACAAATGAATTga                                                             >  1:656890/1‑227 (MQ=255)                                                                                                                                                                  cgcAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTACCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGCCATTAAACCGTGATGTTAcc                                                                                                                                   >  2:64236/1‑158 (MQ=255)                                                                                                                                                                      ggcacctaaaccagttttaaaaaaaTGAAACCCCGCA‑AAAAGTTCCCCCC‑GAAGCCCTGGAAAAAAACAGGCCCCGCTCCCCTGTCAGCCCACCTATCCCGGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTc                                                         <  2:656890/209‑1 (MQ=255)                                                                                                                                                                        ccccccaaacgtgttGAAAAAAATAAACCCCCCCC‑AAATGTTCCCCCC‑AAATGCCGGGCAAAAACCGGGCCGCCCCCCCCTGTAAACCCCCCTATACCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTcac                                                       <  1:464212/215‑1 (MQ=255)                                                                                                                                                                         accaaaatctgttgtgaaaaaaCGAACCACCCCC‑AAAATTTACCACC‑GAAGCCCTGGTAAAAAACAGGCCTCGCACCCCTGGCAGCCCCCCTACGCCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCGCGGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTca                                                        <  1:89685/210‑1 (MQ=255)                                                                                                                                                                                aaaaagagacagaaaTGAACCCCCGCC‑AAAAGTTACCCAC‑GAAGACCTGGTGAAAACCAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGAAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTaa                                               <  2:263939/216‑1 (MQ=255)                                                                                                                                                                                 cTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCAAGACCTCCCTGAACGCCATTAAACCGTTATGTTACCGACCCCCTGACGCGTGAAACAACAATCGTCACAGAACCGGAAAACCCGATCCAATTCACAAATACATAGCAttaaaaaaacaaac                                              >  1:699649/1‑216 (MQ=255)                                                                                                                                                                                  tGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCCAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGCCGCGTGAAAGAATCAGCGTCCGAGAAACGGACAACGCCTTCCAGATCCCAAATGCATTGTATTCACATCATTAAcc                                             >  2:57487/1‑228 (MQ=255)                                                                                                                                                                                   gTTTGTAAAAAATCAACCCCCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAACCAGGCCTCGCTCCCCTGTCAGCCCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTaa                                               <  2:482900/225‑1 (MQ=255)                                                                                                                                                                                    tttGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGACCCAGATCACAAATGCATTGGATTCCCATCATTAACCGt                                           >  2:20018/1‑228 (MQ=255)                                                                                                                                                                                    tttGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACCCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGGATTCACACCCTTAACCGt                                           >  2:451675/1‑228 (MQ=255)                                                                                                                                                                                      tGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGt                                                                                                                                            >  2:415296/1‑129 (MQ=255)                                                                                                                                                                                                       cGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAAAGCAATGAATTCACATCATTAACCGTTTTaaaacgatcccttaa                         >  1:263939/1‑214 (MQ=255)                                                                                                                                                                                                         cc‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGCGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCCCATCATTAACCGGTTTAAGATAATTTCTACACtt                      >  1:53807/1‑228 (MQ=255)                                                                                                                                                                                                                       cacaaaccccgggTAAAAACCGGGCCTGCCTCCCCTGTCAAACCACCAAACCCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAGATGCCTTGTTTTCCCATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCa          <  2:551222/217‑1 (MQ=255)                                                                                                                                                                                                                               ccTGGGAAAACACCGGCCTCGCTCCCCTGCCAGCCCCCCTATCCCCTCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATaa  <  2:53807/228‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              TTTTAGCAGAACAAGAATCTCCTCATTTTCAGAACTTGTTTTTATTTTTCGCCAATATGTTATTGCACTATCACCAATTTATGAATCCCAATGAAAGTGATTTGAATGACGTGTTGATGCCAGCATCATTAAGTGATGATAAAATTATCAAACATATGGCGCGCAGGACCCTCAAACTGTTTGTAAAAAATGAAACACCGCC‑AAAAGTTACTCAC‑GAAGACCTGGTGAAAAACAGGCCTCGCTCCCCTGTCAGACCACCTATACCCGCAACCGCCAAAACGCCAGACCTCCCTGAACGTCATTAAACCGTGATGTTACCGACTCTCTGACGCGTGAAAGAATCAGCGTCAGAGAAACGGAAAACGCGATCCAGATCACAAATGCATTGTATTCACATCATTAACCGTTTTAAGATCATTTCATCACTTTTTCGCAACTCACCCGATAA  >  NC_000913/1616574‑1617021

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