Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,413	A→G	V370A (GTG→GCG)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,413	0	A	G	96.9%	172.4 / ‑4.1	66	V370A (GTG→GCG)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (56/7);  total (59/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

GCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGC  >  NC_000913/73175‑73659                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         gcgcgACCAGCCTTCAATGCCCAGCGACTGACACAACATACGGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCa                                                                                                                                                                                                                                              <  1:409091/247‑1 (MQ=255)  cgcgACCAGCCTTTAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATAGCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAg                                                                                                                                                                                                                                             <  1:487243/248‑1 (MQ=255)         aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                                                                                                                      >  1:15299/1‑251 (MQ=255)         aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                                                                                                                      >  1:701057/1‑251 (MQ=255)             ttCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACAAGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGcc                                                                                                                                                                                                                                  >  2:675116/1‑251 (MQ=255)                aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCcgcc                                                                                                                                                                                                                               >  1:180463/5‑250 (MQ=255)                     ccAGCGACTGACACAACATACTGTAGCGGGCGGTGATCTCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGCAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAACA‑GCCGGCAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCACATCTCAAGCGCCTGACCCGCCCGCATTTTCTGcccaccccc                                                                                                                                                                                                                           >  2:655913/1‑243 (MQ=255)                          gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc                                                                                                                                                                                                                     >  1:628589/1‑251 (MQ=255)                           aCTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCCCCCGAAGTTcg                                                                                                                                                                                                                    >  2:322371/1‑251 (MQ=255)                                   aaCATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCCGCAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCCCCTGACCCCACAGCCTTTTCTGCCGCGccctaaattctcgactgcc                                                                                                                                                                                                            >  2:142582/1‑235 (MQ=255)                                      atacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGCCTGCTcc                                                                                                                                                                                                          >  1:130781/4‑250 (MQ=255)                                       tacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATAAACCTGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCCCGCCCGAAGTTCGCGACTGCTccct                                                                                                                                                                                                        >  2:71379/3‑249 (MQ=255)                                           gTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGGTCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCACTGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGACCAA‑TCCCCGGCATGGCGAGGTTCAACATGCCGCTCATCTCAAGCGCCTGACCCGCCAGCATTTTCAGCCGCGCCCGAAGTTCGCGACGGCTACAtcac                                                                                                                                                                                                     >  2:27616/1‑247 (MQ=255)                                                gggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGACCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATGGCGGCAGGCCAATGGTG‑TTGTTGAGCATAAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                                                                                                                                                                                               <  1:251620/251‑1 (MQ=255)                                                       gATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAAGGATAGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAAGGGTG‑TTGGTGAGCAGTAAAAAGAA‑GCCGGTCGCC‑AGGACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCCGACCAGCCAGCATTTTCTACAGAGCCAGGAGTACGCGACGGCTCCATAGCAGAATCATGGTc                                                                                                                                                                                         >  2:368314/1‑250 (MQ=255)                                                          atCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                                                                                                                      >  1:108139/1‑250 (MQ=255)                                                             gcgcATCGGGTTTTCCAGTACTTTCAGCGCCTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCCCAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGTAAAACGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGCGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCGCGCCCGCACTTCGCGCCCGCTCCCTAGCCGCCTCCTGGTCAGccac                                                                                                                                                                                   >  2:77243/1‑247 (MQ=255)                                                              cgcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTAc                                                                                                                                                                                >  1:168520/1‑251 (MQ=255)                                                                   cGGGGTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGAAAAAAGAA‑GCCGGGAGCC‑AGCACAA‑TCCCCGGCATGGCGGGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CCACGcac                                                                                                                                                                            >  1:124368/1‑250 (MQ=255)                                                                      gTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCCGCA‑CTACGCACAAt                                                                                                                                                                         >  1:492711/1‑250 (MQ=255)                                                                        tttCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCCGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGATCCCTAGCAGCATCATGGTCAGCA‑CTACGCCCAATac                                                                                                                                                                       >  1:710645/1‑250 (MQ=255)                                                                        ttcccagttcttacggcgcATAAGGGATCGC‑ATATAACGTATGGGTAAAAATCACAGTGCCTGCAGCAGATGGCGGCAGCCCAATGGTG‑TTGTTGAGCATAAAAAAGAG‑GCCGGTAGCC‑AG‑ACAACTCCCCGGCATGGCGAGTATCAACATGCCTCCCCTCTCCAGCGCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATac                                                                                                                                                                       <  2:560012/235‑1 (MQ=255)                                                                            cagtaCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCAGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTATGCCGCGCCCGAAGTTCGGGACTGCTCCATAGCAGCATCATGGACAGCA‑CTACGCACAATACACCtt                                                                                                                                                                  >  2:498160/5‑250 (MQ=255)                                                                               taCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTgccc                                                                                                                                                               >  1:170675/2‑250 (MQ=255)                                                                                     cAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGCG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAgcc                                                                                                                                                         >  2:110044/1‑250 (MQ=255)                                                                                        cgaaTAAGGGATCCCCAT‑TACTCCATTGGTAAAAATCACAATGCCGCCAGCAGCTTGCGGC‑GGGCAATGGGGGTGGTTGAACAG‑AAAAAGAAAGCGGGTGCCC‑ACCACAA‑CCCCCGGCATGGCGAGACCCAACATGCCGCTCATCTCCAGCGCTGGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CGACGCACAATACACCTGCCGCCAGCGCaa                                                                                                                                                      <  2:31981/248‑1 (MQ=255)                                                                                         gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTTGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCCCAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                                                                                                                     >  1:315422/1‑251 (MQ=255)                                                                                         gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                                                                                                                     >  1:150998/1‑251 (MQ=255)                                                                                           aTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAc                                                                                                                                                   >  1:534879/1‑251 (MQ=255)                                                                                                 gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACg                                                                                                                                            >  1:378357/1‑251 (MQ=255)                                                                                                 gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACg                                                                                                                                            >  1:361226/1‑251 (MQ=255)                                                                                                 gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACg                                                                                                                                            >  1:172415/1‑251 (MQ=255)                                                                                                             aCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAgc                                                                                                                                 >  1:91573/1‑251 (MQ=255)                                                                                                               gCATTGGTGAAAATAACTATGCCGTCAGCAGATTGCGGCAGGCCCATGGTG‑TTGTTGCGCAGTAAACAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCACTCATCTCCACCGCCTGACCCGCCAGCAATTTATGCCGCGCCCGCAGCTCGCGCATGCTCCCTAACCGCATCATGGTCACCAACAACAC‑CAACACCCCTGACGCCAGCGCCATCC‑CCAACGAGGCCcacacccc                                                                                                                               >  2:618329/1‑247 (MQ=255)                                                                                                                 aTTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCt                                                                                                                             >  1:78648/1‑251 (MQ=255)                                                                                                                    ggTGAAAATCACAATGCCGTCAGCAGATTGCGGGAGGCCAATGGTG‑TTGTAGAGCAGAAAAAAGAA‑GCGGGTAGCC‑AGCACAA‑CCCCCGGCAGGGCGGGGATCAACAGGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGcc                                                                                                                          <  1:289090/251‑1 (MQ=255)                                                                                                                     gagAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCCTGGGCAGCA‑CTACGCACAATACACCTCACGCCCGCGCCACAC‑GCCACGCGGTCCACAgcccccgcc                                                                                                                          >  2:28280/3‑243 (MQ=255)                                                                                                                       gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCACCGAGGGCCACAGCGCCTGCcacc                                                                                                                       >  2:18106/1‑250 (MQ=255)                                                                                                                       gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTg                                                                                                                                                                  >  1:292238/1‑209 (MQ=255)                                                                                                                        aaaaTCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGGCCAAAGCGCCTGCcac                                                                                                                        >  1:39739/1‑249 (MQ=255)                                                                                                                                             gATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCACTAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCAACGGATGTGCCCGCACTTc                                                                                                 >  2:555823/1‑251 (MQ=255)                                                                                                                                                   gggcaggccagtggtgttttggAGCAGTAAAAAGAA‑GCCGG‑AGCCCAGCACAA‑CCCCCGGCATGGGGAGGAAAAACAGGCCCCTCTTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCa                                                                                            <  1:555823/230‑1 (MQ=255)                                                                                                                                                      aGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATCCACATGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAAcgg                                                                                        >  2:156659/1‑249 (MQ=255)                                                                                                                                                         ccaatggtg‑tTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTACCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCCTGGTCAGCA‑CTACGCACCATACACCTGCCGCCAGCGCCATAC‑GCAACGAGGTCCACAGCGCATTCCACACCCCCGGTTGAGCCAGCACTTCCGACAACCGGCg                                                                                     >  2:32313/7‑251 (MQ=255)                                                                                                                                                                            aGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGACGATTTAGCCCATCGACGAt                                                                   >  1:90608/1‑251 (MQ=255)                                                                                                                                                                            aGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGAATAGCGAGGAAAAACATGCCGCTCATCTCCAGCGCCTGACAAGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCCGCGCAATAC‑GCAACGAGGTCAACAGCGCATGCCACAGCACCGGCTGTGCCAGCCATTCCGGCAACTGGCGATTTATCCCATCGACGAt                                                                   >  1:308805/1‑251 (MQ=255)                                                                                                                                                                              tAAAAAGCA‑GCCGAAAGCC‑AGCACAA‑TCCACGTCAGGGCGAGGATAAAAATGCCGCTCATCTCCAACGCCTGACCCGCCAGCATTTTAAGCCGCGCCCGACGATCGCGACTGCTCCTTAGCAGCATCCTGGTCAGCA‑CTACGCACAATACACCTGCCGCCCGCGCAA‑ACAGAAACGAGGTCCACACCGCCGGCAACAGCAACGGTTGAGCCAGCACCTGCGGCAACTGGCCACTTAGCCCATccacgtacc                                                                 >  2:639821/1‑247 (MQ=255)                                                                                                                                                                               aaaaaGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCAc                                                                >  1:253242/1‑251 (MQ=255)                                                                                                                                                                                  aagaa‑gCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGc                                                             >  1:429579/1‑251 (MQ=255)                                                                                                                                                                                             tAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATAAACATGCCGCACATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCTCCCGACGCTCGCGACTGCTCCCTAGCAGAATCATGGTCAGCA‑CTACGCACAATCCACCTGCCGCCCGCGCAATCC‑GCAACGAGCTCCACAGCGCCTGCCACCGCACCGGTTGTGCCAGCACTTCCTGCAAAGGGCGATTTAGCCCACCGACGATCCCCCCCAGCAACGGt                                                   >  2:433258/1‑251 (MQ=255)                                                                                                                                                                                                 c‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGgc                                                >  1:442278/1‑250 (MQ=255)                                                                                                                                                                                                 c‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATAGACGATCACCGCCAGCAACGGTGgc                                                >  2:300624/1‑250 (MQ=255)                                                                                                                                                                                                     cacaA‑TCACCGGCCAGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACACGCCAGCATTTTCTGCCGCGCCAGAAGTTCGCGACTGCACCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACAACTGCCGCCAGCGCAACAC‑GCAACGAGGACCAAAGCGCCGGCCAACGCACCGGTTGTGCCAGCACTTCCGGCCACCGGCGCTTTAGCCCATCAACAATCACCGCCAgacacaggtgaaga                                             >  2:340289/1‑237 (MQ=255)                                                                                                                                                                                                      acaA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGACAGCAACGGTGgcagca                                            >  2:72938/1‑250 (MQ=255)                                                                                                                                                                                                         a‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGACAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAAc                                          >  1:680889/1‑249 (MQ=255)                                                                                                                                                                                                                 gCATGGAGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCAGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGAACTTCCGGCAACGGGCGATTTAGCCCATCGAAGATCACCGCCAGCAACGGAGGCAGCGACagcagcagc                                 >  2:525529/1‑251 (MQ=255)                                                                                                                                                                                                                           gATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGcaca                        >  1:534243/1‑250 (MQ=255)                                                                                                                                                                                                                                   tGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGAGGCCGAAATAGCCGCAGCGCCAGCACAATTAACAcc               >  2:485717/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCg              >  1:220389/1‑251 (MQ=255)                                                                                                                                                                                                                                      cGCTCATCTCCAGCGCCTGACCCGACAGCATTTACTGCCGCGCCCGACGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCCACGGTGGAAGCAACAGCAGCAGCGCCAGCACAATTAACACCgtg            >  1:242181/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCCACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCACCAGCAGCAGCGCCAGCACAATTaaacccgtgc           >  2:591415/1‑243 (MQ=255)                                                                                                                                                                                                                                        cTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCAc                                                                                                    >  1:287573/1‑161 (MQ=255)                                                                                                                                                                                                                                         tCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCCTGGTCAGCA‑CTACGCACACTACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGTTTGTGCCAGCACTTCCGGCAACTGGCGATTTAACCCATCGACGATCACCGCCAGCACCGGTGGCAGCAACCGCAGCAGCGCAAGcgaaactagacccgtgtcc         >  2:496099/1‑233 (MQ=255)                                                                                                                                                                                                                                           aTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCACCGCCGATCACCGCCAGCAACGGTGGCAGCCACAGCAGCAGCGCCAGCACACTTAACACCGAGTcccc       >  2:647953/1‑248 (MQ=255)                                                                                                                                                                                                                                            tctcCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGCTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCGTCGACGCTCACCGCCAGCAACGGAGGCAGCACCAGCAGCAGCGCCCGCCCAATTAACACCGTGTCGCa       >  2:97435/1‑250 (MQ=255)                                                                                                                                                                                                                                                cAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTAGCCCATCGACGATCACCGCCAGCAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGc  >  1:505423/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         GCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCC‑AGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAC‑GCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTCCGGCAACTGGCGATTTACCCCATCGACGATCACCGCCAGTAACGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGC  >  NC_000913/73175‑73659

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