Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	976,458	G→A	A45T (GCG→ACG)	mukB →	chromosome partitioning protein MukB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	976,458	0	G	A	96.5%	163.9 / NA	59	A45T (GCG→ACG)	mukB	chromosome partitioning protein MukB
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (28/28);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGTGAAGCACAGCGTCGCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGGCGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTG  >  NC_000913/976211‑976703                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        cccGTGAAGCACAGCGTCGCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATgacgg                                                                                                                                                                                                                                                      <  2:338505/251‑1 (MQ=255)             aGCGTCGCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACgg                                                                                                                                                                                                                                          <  1:402944/251‑1 (MQ=255)                 tCGCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGGAAATCCCCCACCCTGACGGCGTTCGTTACGCCGCt                                                                                                                                                                                                                                      >  2:9477/1‑251 (MQ=255)                  cGCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCGTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCt                                                                                                                                                                                                                                      <  2:63254/250‑1 (MQ=255)                   gCCTGATTAGTGAGGGCGAAACAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGCGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATACAACACCATGACGCCGTTCGGTACggggggg                                                                                                                                                                                                                                     >  2:343134/1‑245 (MQ=255)                        aTTCGTGATGTCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTTGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATcccc                                                                                                                                                                                                                               <  1:50807/251‑1 (MQ=255)                        aTTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATcccc                                                                                                                                                                                                                               <  1:207395/251‑1 (MQ=255)                              gATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTg                                                                                                                                                                                                                         >  2:104338/1‑251 (MQ=255)                                   agaagaaATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCAAGATAGCGGAGAGGAAGAATAATGATTGAAAGCGGTCAATTTCGCTCACTGAAGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTAACGACGCTTTCTGGCGGTAACGGGGCGGGGAAAACCAACACCCTGACGTCGTTTGTTACGGCGCTAATCCCCGACCTGccccc                                                                                                                                                                                                                    >  2:387855/2‑250 (MQ=255)                                      aGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTAAGGCGCATAACCCCGACATGACCctgat                                                                                                                                                                                                                 >  2:334514/1‑249 (MQ=255)                                       gCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCctgctg                                                                                                                                                                                                                <  2:45217/251‑1 (MQ=255)                                                       atcatcTGCGACTCAACGATGAAACCGGAGAGAATCAGCCAGATAGCGGCGAGGAATAATAATGCTTGAAAGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCAAGACGGCGGTCGTTAAGGGGCTTATACACGGCCTGAGCCTGTCGCATTGCCGTAACACt                                                                                                                                                                                                 >  2:82291/1‑250 (MQ=255)                                                              gCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTGGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGTCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCGTTTCCGTATCACTACGGaa                                                                                                                                                                                           <  2:271228/249‑1 (MQ=255)                                                                aaCTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTAACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGGAAATCCCCCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGACTTTCCGTAACACTACGGAAGccc                                                                                                                                                                                       >  2:176511/1‑250 (MQ=255)                                                                     aaCGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGGGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCTGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGTGGGGAAATCCACCACGCTGACGGCGTTCGTTACGGCGCTGATCCCCGACCTTAACCTGATGCATTTCCGTAACACTACGGAAGCCGGGGcc                                                                                                                                                                                  >  1:52988/1‑251 (MQ=255)                                                                       cGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGcc‑ac                                                                                                                                                                               >  1:128428/1‑251 (MQ=255)                                                                        gATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGTCGCTGATTAACTGGAACGGCTTGTTTGCCCTAACTTTTGACCTTGACGAGCTGGTCACGACGCGTGCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGcc‑acc                                                                                                                                                                              <  1:41936/251‑1 (MQ=255)                                                                         aactaaCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGcc‑acca                                                                                                                                                                             <  2:67132/247‑1 (MQ=255)                                                                         aTGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGcc‑acc                                                                                                                                                                              >  1:301388/1‑250 (MQ=255)                                                                            aaaCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGGAAATCCACCAACATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGAACATGCTGCATTTCCGTAACACTACGGAAGGCGGGGCC‑ACCAGCg                                                                                                                                                                          >  2:53214/1‑251 (MQ=255)                                                                             aaCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGCCGGCGTTCGGTACGGCGCTGATCCCCGACCTGACCCTGCTGCAATTCCGTAACAACACGGAAGCCGGGGcc‑accc                                                                                                                                                                             >  1:67132/1‑246 (MQ=255)                                                                             aaCCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCCTTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCgg                                                                                                                                                                         >  1:308241/1‑251 (MQ=255)                                                                                 gAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTcg                                                                                                                                                                     <  1:73172/251‑1 (MQ=255)                                                                                    gagaATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCCTGACGGCGGTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCCTTTCCGGAAAACTACCGGAGCCGGGGTC‑CCCAACGGTTagcgc                                                                                                                                                                  >  2:358950/1‑251 (MQ=255)                                                                                    gagaATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTcgcgc                                                                                                                                                                  <  2:141471/251‑1 (MQ=255)                                                                                      gaATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATTCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCCCACC‑GCGGTTcgcgcgg                                                                                                                                                                >  2:59725/1‑250 (MQ=255)                                                                                                    gCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCGCCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACg                                                                                                                                                  <  1:210525/251‑1 (MQ=255)                                                                                                    gCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTAACTGACGCTGATTAAGTGGAGCGGCTTTTTTGCCCGAACTTGTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCGTGGCGGCGTTCGTTACGGCGCTGATCCCCGACCTGAGCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCaag                                                                                                                                                  <  2:202213/251‑3 (MQ=255)                                                                                                     cGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACgg                                                                                                                                                 <  2:122484/251‑1 (MQ=255)                                                                                                             aagaaTAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTg                                                                                                                                         >  1:215239/1‑250 (MQ=255)                                                                                                                 ataatGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGc                                                                                                                                    <  1:244109/251‑1 (MQ=255)                                                                                                                 ataatGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGATGGTCACGACGATTTCTGGCGGAAACGGGGCGGGGAAACCAACCACCATGACGGCGTTTGGTAAGGCGCTGAAACCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACAAGCGGTTAGCGCGATAAAGGGATGAAAGGT‑AAGCTGAAAGc                                                                                                                                    >  1:53274/1‑251 (MQ=255)                                                                                                                  taatGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGCAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCATCATGACGGCGTTCGTTTCGGCGCTGCTCCCCGACCCGACCCTGCCGCATCTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCCCGATAAAGGTCTGCACGGTTAAG‑TGAAAGCg                                                                                                                                   >  2:357925/1‑251 (MQ=255)                                                                                                                      gATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGgtg                                                                                                                               <  1:401466/251‑1 (MQ=255)                                                                                                                               cgGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTAtt                                                                                                                      <  2:332679/251‑1 (MQ=255)                                                                                                                                   aaaTTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AATCTGAAAGCGGGTGTCTGTTATTCGAt                                                                                                                  <  1:165773/251‑1 (MQ=255)                                                                                                                                        tCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTc                                                                                                              <  1:104338/250‑1 (MQ=255)                                                                                                                                             cACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGCAAATCCACCACCATGACGGCGTTCGTTAACGCGCTGCTCCCCGACCTGAACCTGATGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGCTCGCTCGCTAAAGGTCTGCACGGT‑CAACCGCAAGCGGGTGCCCGTTTTTCGACGCTCgccccc                                                                                                        >  2:442039/1‑246 (MQ=255)                                                                                                                                                   cGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGGAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGAACTGACACTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGACTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTTGATGCTCGACACCATTAAc                                                                                                  >  1:6672/1‑251 (MQ=255)                                                                                                                                                       gATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTAGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGTGTGTCTGTGATTCGATTCTCGACACCATTAACTcgc                                                                                              <  2:268566/251‑1 (MQ=255)                                                                                                                                                                gAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAgc                                                                                     <  1:365311/251‑1 (MQ=255)                                                                                                                                                                      cTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCgtgg                                                                               >  1:23757/1‑251 (MQ=255)                                                                                                                                                                                 gagcTTTTGACCTTGACTAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCCCCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGtgtg                                                                    <  2:845/248‑1 (MQ=255)                                                                                                                                                                                     ttttGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCCGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTc                                                                >  1:81534/1‑251 (MQ=255)                                                                                                                                                                                                        gTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCCAGACGGCGATCGTTACGGCGCTGATCCACGACCTGACCCTGCTGCATTTCCGTCACCCTAAGGAAGCCGGGGCC‑ACCAGCGGTTCGCACGATAAAGGTCTGCACGGT‑AAGCAGAAAGCGGGGGCCTGTTATTCGATGCTCCACACCCTTAACTCGCGCCCCCCGCGCGTGGTGGTGGGTGGGCGTCTGCAACATGTTGCCGGAcg                                             >  1:246772/1‑251 (MQ=255)                                                                                                                                                                                                             gccgCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGTTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGTTGGTCGGTTTGCGTCTGCAACAGGTTGCCGGACGCGATc                                        <  1:173801/249‑1 (MQ=255)                                                                                                                                                                                                              acgCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCCCCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCCGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACCGGTTGCCGGACGCGATCg                                       >  2:44104/1‑251 (MQ=255)                                                                                                                                                                                                               cgCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGt                                      <  1:375433/251‑1 (MQ=255)                                                                                                                                                                                                                   ttCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTaa                                    <  2:131161/249‑1 (MQ=255)                                                                                                                                                                                                                       ggCGGTAACGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGt                                 <  1:299487/248‑1 (MQ=255)                                                                                                                                                                                                                             aaCGGGGCGGGAAAATCCACCACCATGACCGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGGTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGGGTGCGTCTGCAACAGGTTGCCGGACGCGTTCGTAAAGTGGATATCa                         >  1:378716/1‑250 (MQ=255)                                                                                                                                                                                                                               cGGGGCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGc                      <  2:287268/251‑1 (MQ=255)                                                                                                                                                                                                                                 gggCGGGTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCg                    <  1:389943/251‑1 (MQ=255)                                                                                                                                                                                                                                       gTAAATCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAACGAGGCTATCAAGCCGGTTGcc              >  2:277937/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaaTCCGCCACCATGGCGGCGTTCGTTACGGCGCTGATCCCCGACCTGACACTGCGGCAGTTCCGTAACAATACGGAAGCCGGGGAC‑ACCAGCTGTTCGCGCGATAAAGGTCTGCTCGGT‑ATGCTGATAGCGGGTGTCTTTTATTTGATGCGCGACACCATTAACTCGCGCCGCCAGCGCGTGGTGGTCTCTGTGCATCGGCAACATGTTGTCGGGCGAGATCGTAAAGTGGatatgaatctgttggccat            <  2:63203/251‑16 (MQ=255)                                                                                                                                                                                                                                          aaTCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGGGGATAGCAAGCCGTTTGCCAtt           <  2:301388/251‑1 (MQ=255)                                                                                                                                                                                                                                           aTCCACCACCATGACGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGTCACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTc          >  2:302909/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCCACCACCATGACCGCGTTCGTTACCGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGAC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGCAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTCACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTc          >  2:214366/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ccaTGACGGCGTTCGTTACGGCGCTGCTCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGGGGATATCAAGCCGTTTGCCATTCCGGGACTg  >  2:441683/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ccaTGACGGCGTTCGGTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCGGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTg  <  2:149228/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        CCCGTGAAGCACAGCGTCGCCTGATTCGTGATGGCGAAGCAATGCCGATTGAAAATCATCTGCAACTCAACGATGAAACCGAAGAGAATCAGCCAGATAGCGGAGAGGAAGAATAATGATTGAACGCGGTAAATTTCGCTCACTGACGCTGATTAACTGGAACGGCTTTTTTGCCCGAACTTTTGACCTTGACGAGCTGGTCACGACGCTTTCTGGCGGTAACGGGGCGGGTAAATCCACCACCATGGCGGCGTTCGTTACGGCGCTGATCCCCGACCTGACCCTGCTGCATTTCCGTAACACTACGGAAGCCGGGGCC‑ACCAGCGGTTCGCGCGATAAAGGTCTGCACGGT‑AAGCTGAAAGCGGGTGTCTGTTATTCGATGCTCGACACCATTAACTCGCGCCACCAGCGCGTGGTGGTCGGTGTGCGTCTGCAACAGGTTGCCGGACGCGATCGTAAAGTGGATATCAAGCCGTTTGCCATTCAGGGACTG  >  NC_000913/976211‑976703

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