Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,950,201	0	T	G	30.3%	135.8 / 10.0	66	G320G (GGT→GGG)	argA	N‑acetylglutamate synthase
Reads supporting (aligned to +/- strand):  ref base T (15/31);  new base G (20/0);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 8.20e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCAAATTGTTTGTGCTGACCACGCGC  >  NC_000913/2949953‑2950450                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           gAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTAt                                                                                                                                                                                                                                                          <  2:171288/251‑1 (MQ=255)        aaGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTTATTCGCCCACTGGAGCAGCAAGGGATTCTGGt                                                                                                                                                                                                                                                    >  1:16281/1‑250 (MQ=255)             gATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGGTTGATTCGCCCCCTGGGGCCGCAAGGTATTTTGGGAC‑Gccc                                                                                                                                                                                                                                             >  1:288697/1‑250 (MQ=255)             gATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCCCTGGAGCAGCAAGGGATTTTGGTAC‑Gccc                                                                                                                                                                                                                                             >  1:80073/1‑250 (MQ=255)                   aaCTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTTTGGGTTTATTCGCCCACTGGAGCAGCAAGGGATTCTGGTAC‑GCCCTtctcc                                                                                                                                                                                                                                       >  1:72074/1‑250 (MQ=255)                   aaCTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGGTTTGTTCGCCCACTGGAGCAGCAAGGGAATTTGGGAC‑GCCGTtctcc                                                                                                                                                                                                                                       >  2:194555/1‑250 (MQ=255)                             cGGTGCGCTTTTCGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTg                                                                                                                                                                                                                             <  1:493239/251‑1 (MQ=255)                               gTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGGGTGATTACCCCACCGGGGGAGCAAGGGCTTCTgggaaccccttcttgcgagcaagtggg                                                                                                                                                                                                                            >  1:235763/1‑226 (MQ=255)                                      ttttGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATgg                                                                                                                                                                                                                      <  2:86178/249‑1 (MQ=255)                                                 gcAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAAtt                                                                                                                                                                                                         <  2:216667/251‑1 (MQ=255)                                                       aaaGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTcgc                                                                                                                                                                                                                                      <  1:93758/216‑1 (MQ=255)                                                                  cAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTTGGGGTTGATTTGCCCACTGGGGCAGGAAGGGATTCTGGGAC‑GCCGGTCTCGCGAGCAGCTGGGGATGGGAAACGACAAAATTACCATTATTTAGcgcg                                                                                                                                                                                        >  1:188428/1‑251 (MQ=255)                                                                       gcgtgcgtCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGAt                                                                                                                                                                                      <  1:337147/248‑1 (MQ=255)                                                                       gcgtgcgtCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCTGGAGTTTATTCGCCCACTGGAGCAGCAAGGGGTTTTTGTACCGCGGTT‑TCGCGCGCCGGTGGGGATGGGAATCGACAAAATCCCCATTATTCAGcgcgcatg                                                                                                                                                                                    >  2:337147/1‑246 (MQ=255)                                                                             gtCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATc                                                                                                                                                                                                                 <  2:414389/215‑1 (MQ=255)                                                                               cGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTAt                                                                                                                                                                           <  1:487128/251‑1 (MQ=255)                                                                                 cTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGGTGATTCGCCCACTGGAGCAGCAAGGGATTTTGGTAC‑GCCGGTTTCGCGAGCAGCCGGGGATGGAAATCGACCAATTCACCATTATTCAGCGCGATAACACGACTAttt                                                                                                                                                                         >  2:557207/1‑250 (MQ=255)                                                                                   gTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGc                                                                                                                                                                        <  2:96696/250‑1 (MQ=255)                                                                                         tAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTgcgc                                                                                                                                                                  <  1:91698/250‑1 (MQ=255)                                                                                          aaTCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGGTTGATTCGCCCACTGGGGCAGCAAGGTATTTTTGTAC‑GCCCTTCTCGCGAGCCGCTGGAGATGGAAATCGACAAATTCAACATTATTTAGCGCGATAACACGACTATTGCCTGcgccg                                                                                                                                                                >  1:25836/1‑251 (MQ=255)                                                                                              aGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTTGGAGTTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTAC‑GCCGTTTTCTCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGGTAAAACGACAATTGCCTGCGCCgcgcc                                                                                                                                                            >  2:413152/1‑250 (MQ=255)                                                                                                   tCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATc                                                                                                                                                       >  1:27184/1‑251 (MQ=255)                                                                                                    cAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATcc                                                                                                                                                      <  2:80073/251‑1 (MQ=255)                                                                                                          gATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACCAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCCCGCTCTATCCGTTccc                                                                                                                                                >  1:447768/1‑251 (MQ=255)                                                                                                              gcgcgcTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGgaa                                                                                                                                            <  1:194555/251‑1 (MQ=255)                                                                                                                 cgcTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCgggag                                                                                                                                           >  1:13123/1‑246 (MQ=255)                                                                                                                      ttGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCCTTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATACCTTCCCGgaagaga                                                                                                                                        >  1:439825/1‑247 (MQ=255)                                                                                                                        gCAAGAGTTGTTCTCACGCGCCGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCTGGAGGTGATTCGCACACTGGAGCAGCAAGGGAATATGGTAA‑GCCGTGCTCGCCGGCAGCTGGAGATGGGAATCGACAAATTCACCATTATTCTGCCCGATAAAACGACTATTGCCTGCGGCGCCCCCTATTCGTTCCCGgcagagccgattt                                                                                                                                  >  2:180588/1‑244 (MQ=255)                                                                                                                         cAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTAC‑GCCGTTCTCGCGAGCCGCTGGGGGTGGAAAACGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGCAGAGGAGATTgg                                                                                                                                 >  2:139293/1‑251 (MQ=255)                                                                                                                          aaGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACCGGAGCAGCAAGGGATTCTGGTAC‑GCCGCTCTCGCGCGCAGCTGGAGATGGAAAACGACAAATTTACCCTTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATTCGTTCCCGGAAGAGAAGATTggg                                                                                                                                >  2:413287/1‑251 (MQ=255)                                                                                                                            gagTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGa                                                                                                                              <  1:479764/251‑1 (MQ=255)                                                                                                                                tgttCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGGGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGAGTGGGGAAATg                                                                                                                          >  1:336918/1‑251 (MQ=255)                                                                                                                                    ctcACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGGATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGGAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCt                                                                                                                      >  1:133122/1‑251 (MQ=255)                                                                                                                                      cACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGcc                                                                                                                       <  1:462339/248‑1 (MQ=255)                                                                                                                                         gcgACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCtgtgtg                                                                                                                 <  2:288697/251‑1 (MQ=255)                                                                                                                                          cgACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTgg                                                                                                                <  2:45774/251‑1 (MQ=255)                                                                                                                                                    ggTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGGGTTGATTTGCCCCCCGGGGCAGCAAGGGATTCTGGGAC‑GCCGTTCTCGCGCGCCGCCGGGGATGGGAAACGACAAAATAAACATTATTCAGCGGGATAACACCACTATTGCCTGCGCCGCGCGCTCTCCGTTCCCGGACGAGGAGAGTGGGGGAATGGcccggggtgcaggtccccc                                                                                                      >  2:166637/1‑234 (MQ=255)                                                                                                                                                    ggTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCTGGGGTTGGTTCGCCCCCTGGAGCAGCAAGGGATTCTTGTAC‑GCCGTTCTCCCCAGCAGCTGGAGATGGAAATTGACAAAATCACCATTATTTAGCGCGACAAAACGACTATTGCCCGCGCCGCGCTCTCTCCGTTCCCGGGAGAGAAGGTTGGGGAAATGGCCTGGGTGGCAGtttcacc                                                                                                      >  2:66878/1‑246 (MQ=255)                                                                                                                                                       aCGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTTGCCCACTCGAGCAGCAAGGTATTCTGGTAC‑GCCGTTTTCGCGAGCAGCTGGAGATGGAAATCGGCAAATTCACCATTATTCAGCGCGATAACACAACTATTGACTGCGCCCCGCTCTATCCGTTCCCCGAAGAGAAGATTGGGGAAATGGTCCGGGTGGCCGTTCACACgg                                                                                                    >  1:45452/1‑250 (MQ=255)                                                                                                                                                                gtgATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACTCGAATTACCGCa                                                                                           <  1:386622/250‑1 (MQ=255)                                                                                                                                                                        aaGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGAtt                                                                                                 >  1:502345/1‑236 (MQ=255)                                                                                                                                                                        aaGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGAtt                                                                                                 <  2:502345/236‑1 (MQ=255)                                                                                                                                                                            gcCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTTTTCTGGAGTTGATTCGACCACTGGAGGAGCCAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGGTGGGGAGGGGAATCGACAAATTCCCCATTATTCAGCGCGATAACACGACAATTGCCTGCGCCGCGCTCTATCCGTTTCCGGAAGGGAAGATTGGGGGAATGGGCTGTGTGGCAGTTCAACCGGATTACCGCAGTTCATCAAgggg                                                                              >  1:235233/1‑251 (MQ=255)                                                                                                                                                                             ccGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGt                                                                             <  1:139293/251‑1 (MQ=255)                                                                                                                                                                                   aGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGtt                                                                       <  2:416858/251‑1 (MQ=255)                                                                                                                                                                                        cgtcgCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGATCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGTGTTAAGTTctgc                                                                   <  2:102621/250‑1 (MQ=255)                                                                                                                                                                                            gcgcAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTg                                                                                                                                                                     <  1:6424/148‑1 (MQ=255)                                                                                                                                                                                            gcgcAACAATCAACGATATTGGCGGTATTCTGGAGGTTATTCGCCCCCTGGAGCAGCAAGGGATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGGGATGGAAATCGACAAATTCACCATTTTTCAGCGCGATAACACGACTATTGCCTg                                                                                                                                                                     >  2:6424/1‑148 (MQ=255)                                                                                                                                                                                              gcAACAATCAACGATATTGGCGGGATTCTTGAGTTGATTCGCCCACACGAGCAGCAAGGGATTCTTGTAA‑GCCGTTTTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAAAACGACTATTGCCTGCGCCGCGCTCTATCCGGTCCCGGGAGAGAAGGTTGGGGGAATGGCCGGTGTGGCAGTTCCCCCGGGTTACCGCAGTTCATCAAGGGGGGAAGGTTTGCTGGAAc                                                             >  2:393630/1‑250 (MQ=255)                                                                                                                                                                                                 acaaTCAACGATAGTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCAt                                                         <  2:188428/251‑1 (MQ=255)                                                                                                                                                                                                  caaTCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCCCTGGAGCAGCAAGGGAGTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCCGTGTGGCAGTTTCCCCCGATTACCGCAGTTCATCAAGGGGGGAAGTTCTGGTGGAACGGAtt                                                        >  2:116057/1‑251 (MQ=255)                                                                                                                                                                                                       aaCGATATTGGCGGTATTCTGGGGTTGATTCGCCCAATGGAACAGCAAAGGATTCTGGTAA‑GCCGTTCTCGCGGGCAGCTGGGGGTGGAAATTGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCGCTATCCGGTCCCGGGAGAGAAGATTGGGGAAATGGCCCGTGTGGCAGTTCACCCGGGGTAACGCAGTTTATCAAGGGGGGAAGGTCTGCTGGAACGCCTTGcccc                                                   >  2:301401/1‑249 (MQ=255)                                                                                                                                                                                                           atatTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGATCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAg                                               <  2:350012/251‑1 (MQ=255)                                                                                                                                                                                                               tGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCg                                            <  2:173110/250‑1 (MQ=255)                                                                                                                                                                                                                  cGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGc                                        <  2:27184/251‑1 (MQ=255)                                                                                                                                                                                                                    gTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGc                                        >  2:131337/1‑249 (MQ=255)                                                                                                                                                                                                                         cTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGtt                                                                                                                                                   <  2:42706/137‑1 (MQ=255)                                                                                                                                                                                                                         cTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTAACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGtt                                                                                                                                                   >  1:42706/1‑137 (MQ=255)                                                                                                                                                                                                                           ggAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGTTAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCGCAGGCGAAGCAGAGCGGc                                <  2:439752/250‑1 (MQ=255)                                                                                                                                                                                                                           ggAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCg                                  <  2:31721/248‑1 (MQ=255)                                                                                                                                                                                                                            gAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGGAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCtt                              >  2:7974/1‑251 (MQ=255)                                                                                                                                                                                                                            gAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCtt                              <  1:301401/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCa                         >  1:475000/1‑250 (MQ=255)                                                                                                                                                                                                                                    tCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCaaa                       <  2:16281/250‑1 (MQ=255)                                                                                                                                                                                                                                               agcagcAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCAAATTGTTTGTGCTg           >  2:198555/1‑251 (MQ=255)                                                                                                                                                                                                                                                        gTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCAAATTGTTTGTGCTGACCACgcgc  >  1:36380/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           GAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTGGAGTTGATTCGCCCACTGGAGCAGCAAGGTATTCTGGTAC‑GCCGTTCTCGCGAGCAGCTGGAGATGGAAATCGACAAATTCACCATTATTCAGCGCGATAACACGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCCTGTGTGGCAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTTAAGCAAATTGTTTGTGCTGACCACGCGC  >  NC_000913/2949953‑2950450

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