Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	46,634	G→T	G276G (GGG→GGT)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	46,634	0	G	T	91.5%	125.7 / ‑0.4	48	G276G (GGG→GGT)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (21/22);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTG  >  NC_000913/46385‑46868                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           gCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTgggg                                                                                                                                                                                                                                            <  1:216809/250‑1 (MQ=255)     gAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTgg                                                                                                                                                                                                                                       >  1:357670/1‑251 (MQ=255)      aaTCTCCCCTCTGGTTAGTACGAAAAGGGCGAGTAACAGAGTGCGAAGAGAGGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGTAGCCGCAGAAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGc                                                                                                                                                                                                                                      <  2:392263/251‑1 (MQ=255)           ccccGCTGGTTATGACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACTCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCaaaaa                                                                                                                                                                                                                                 <  1:143145/251‑1 (MQ=255)                     tattaCGCAAAGGGCGAGGAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGCTCACAATGTTCGCCATTTACACCTTTGGCCCGCAATTCGTTGGTTTGTTGGGTGTGGGTGTTGGCAAAAATGCGGcgctg                                                                                                                                                                                                                       >  2:126571/1‑247 (MQ=255)                      attaCGCAAAGGGCGAGTAAAAGAGAGCGAAGAGATGATACTAATACTGTTTGGCGAACCGAGGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTAGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGCCCTGCCAGGTGATCCCAATGTTCGCCATTTACACCATTGGCCAGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCACAAACGCGGcgctgg                                                                                                                                                                                                                      >  1:279181/1‑246 (MQ=255)                          cGCAAAGGGGGAGTAAAAGAGTGCGAATAGATGATGTGAAAACTGTGTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGTTCTGTTTAATCGCCGCCATTTTCCTTTTTTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTTATCCCAATGTTCGCCATTTACACCTTTGGCTCGCAAATCGTTGGTTTGGTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGaa                                                                                                                                                                                                                  <  2:195933/251‑1 (MQ=255)                              acgggCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGATTTCGATGAAGAGCAGCCGAAGCAAACCCGTTTTCGCGATATGTTTAATCGCAGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGCTCCCAATGTTCGCCATTTACCCCTTTGGCCCGCAAAACGTTGGATTGTTGGGTTTGGGTGTTGGCCAAAACGCGGCGCTGGGGAAcgtg                                                                                                                                                                                                              >  2:458736/3‑251 (MQ=255)                                    aaGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGTTTGGGTGTTGGCAAAAACGCGGcgctgggggacggggtgatt                                                                                                                                                                                                        >  1:266714/2‑232 (MQ=255)                                      gTAAAAGAGTGCGAAGAGATGATGATAAACCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAg                                                                                                                                                                                                      <  1:156947/251‑1 (MQ=255)                                        aaaaGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGCGTTGGCACACACGCGGCGCTGGGGAACGTGGTGATTAGcc                                                                                                                                                                                                    >  2:132070/1‑251 (MQ=255)                                             agTGCGCACAGATGATGATAAAACTGATTGGCAAACCGGAAGCGTTCGATGAAGAGCAGCCGCAGCAAACCAGATTTCCAGCTATGTTTAATCGCCACAATTTTCCTCTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATACCAATGTTCGCCATTTACACCGTTGGCCCGTAAATCGTTGGTTTTTTGGGATTGGGTGGTGGCGAAAACGCGGCGCTGGTGAGCGTGGCGATTAGCCTGttt                                                                                                                                                                                               >  2:42263/1‑250 (MQ=255)                                                gCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCtt                                                                                                                                                                                             <  1:458258/250‑1 (MQ=255)                                                              ataaaaCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTAt                                                                                                                                                                              >  2:157666/4‑251 (MQ=255)                                                              aaaaaaCGGTTGGGCGAACCGGGGGCTTGCGATGAAGACCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTTCACCTTTGGCCCGCAAATCGTTGGTTTTTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTAt                                                                                                                                                                              <  2:54306/248‑1 (MQ=255)                                                                      gTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGCACGTGGTGCTTAGCCTGTTCTTTATGCTCGGCTGTATTccgccgg                                                                                                                                                                      >  2:384537/1‑250 (MQ=255)                                                                                   gTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                          >  1:467654/1‑218 (MQ=255)                                                                                                  gagCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCg                                                                                                                                          <  1:384537/251‑1 (MQ=255)                                                                                                            agcaATCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGGGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCAttgttg                                                                                                                                >  1:390676/1‑251 (MQ=255)                                                                                                             ccaAATCCGTTTTCGGGTTTTGTCTAATCGCCGCCATTTTCCTTTGGTGCTTGTTGTTGCCGCCATCTGGGCCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTTCATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGa                                                                                                                               <  2:457891/250‑1 (MQ=255)                                                                                                                aaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTg                                                                                                                            <  1:114213/251‑1 (MQ=255)                                                                                                                     gTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGc                                                                                                                       >  1:182811/1‑251 (MQ=255)                                                                                                                                   tttAATCGCCGCCATTCTCCTTTTGTTCTGTTTGTTGCCGCCATCTGTCCCTGCCAGGTGATCCCAATGTTCTCCATTTCCACCTTTGGCCCGCACATCGTTGGTTTGTTGGGCTTGGGTGTTGGCCCAAACGCGGCGCTGGGGAACGTGGTGATTCGCCTGTTCTTTATGCTCGGCTGTCTTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCTTCCATTGTTGATTGGCTGCTTTGCCCTGATGAc                                                                                                         >  2:27094/1‑251 (MQ=255)                                                                                                                                           cctccATTTTCCTTTTGTTCTGTTTGTTGCCACCATCTGGACCTGCCAGGTGCTCCCAATGTTCGCCATTTACACCTTTTGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTAGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGCTGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTTGCAGCTTGTCCATGATGAcgctggcgc                                                                                                 >  1:5975/4‑251 (MQ=255)                                                                                                                                            cgccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGt                                                                                                                                                                                                 >  2:121519/1‑154 (MQ=255)                                                                                                                                            cgccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGt                                                                                                                                                                                                 <  1:121519/154‑1 (MQ=255)                                                                                                                                                     ccTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGtt                                                                                                                                                            <  1:175816/182‑1 (MQ=255)                                                                                                                                                     ccTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGtt                                                                                                                                                            >  2:175816/1‑182 (MQ=255)                                                                                                                                                                gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                          <  2:311505/141‑1 (MQ=255)                                                                                                                                                                gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                          >  1:311505/1‑141 (MQ=255)                                                                                                                                                                gtttgtgGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTGGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCTATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCtgat                                                                              <  1:2297/249‑4 (MQ=255)                                                                                                                                                                        ccgccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTaa                                                                                                                                                          <  1:197827/165‑1 (MQ=255)                                                                                                                                                                        ccgccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTaa                                                                                                                                                          >  2:197827/1‑165 (MQ=255)                                                                                                                                                                             aTCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGGTTTgggggtggttccggatattggg                                                                 >  2:196207/1‑231 (MQ=255)                                                                                                                                                                                    ccTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGACTGGCAGCTTTGcc                                                                                                                 <  2:413892/194‑1 (MQ=255)                                                                                                                                                                                    ccTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGACTGGCAGCTTTGcc                                                                                                                 >  1:413892/1‑194 (MQ=255)                                                                                                                                                                                    catggcaggggttCCCAAGGTTCGTCATTTACACCTTTGGCCCGCAAATCGTGGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGAACGGGGTGATTAGCCTGTTCTTTAGGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATctggctg                                                        <  2:350564/239‑1 (MQ=255)                                                                                                                                                                                                    aaTGTTCGCCATTTACACCTTTGGCCCGCAAATCGTGGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAGGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATGTGGGGATCTGGCTGGTAGTGATGGCCTTTg                                        <  2:390676/251‑1 (MQ=255)                                                                                                                                                                                                      tGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGCTTTGTTGGGATTGGGTGTTGGCCAAAACGCGGCGCTGGGGACCGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGCGGCTGGGGTTAAACACTGCCGGACGGAGTCCATTGTTGATTGGCAGCTTTGCCATGATGCCGCTGCCGCGGGCGGTTTTGGGGGTGGTCCCGGGTGTGGGGAttggggtggggtgggtgggctttgtg                                      >  2:37219/1‑226 (MQ=255)                                                                                                                                                                                                          cGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGgtgt                                  <  1:157666/251‑1 (MQ=255)                                                                                                                                                                                                                ttacacctgtgcggcGCAAATCGTTGGTTTGTTGGGATTGGGTGTTAGCAATAACGGGGCGTTGGGGAGCGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGGATGCCGCCGATGCTCTGGTTAACCACTGGCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGATGGCGGTTTTGTGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTGGCGGTGTATGcc                             <  1:475279/236‑1 (MQ=255)                                                                                                                                                                                                                  acacCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTGATGCTCGGCTGTATTCCTCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCttt                          <  1:390018/251‑1 (MQ=255)                                                                                                                                                                                                                            cccGCAAATTGTCGGATTGGTGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGGACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTATGCCATGATGACGCTGGCGCTGGCGGATTTGGGGGTGAGCCCGGGTATGGGGACCTGGTGGGGGGGGAGTGCCTTTGTGGTGTATGCCtttttgttctgcg                >  2:83036/1‑243 (MQ=255)                                                                                                                                                                                                                                     cGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGAt                                                                                                                              >  2:140602/1‑132 (MQ=255)                                                                                                                                                                                                                                     cGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGAt                                                                                                                              <  1:140602/132‑1 (MQ=255)                                                                                                                                                                                                                                      gTTGGTTTGTTGGGATTGTGTGTTGGCAAAAACGCGTGGCTGGGGAAGGTGGTGAGTAGCCTGTTCTTTATTCTCGGCTGTATTCCGCCGATGCGTTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCTGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTaa      <  2:388897/251‑1 (MQ=255)                                                                                                                                                                                                                                        tgggttgttgGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATgg                                              <  1:368741/205‑1 (MQ=255)                                                                                                                                                                                                                                        tGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTTGGGCTGATCCCGGGTATGGGGGTCTGGCTGGTGGTGATGGGCTTTTCGGTGTAAGCCTTTTTTTCTGGCGGGGCGGGGAAtt    >  1:136036/1‑251 (MQ=255)                                                                                                                                                                                                                                          gTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTg  <  1:403237/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                           GCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTG  >  NC_000913/46385‑46868

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