Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,367,944	G→A	A40T (GCG→ACG)	yhcG →	DUF1016 domain‑containing protein YhcG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,367,944	0	G	A	97.7%	135.8 / ‑4.0	43	A40T (GCG→ACG)	yhcG	DUF1016 domain‑containing protein YhcG
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (13/29);  total (14/29)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

CAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAG  >  NC_000913/3367701‑3368190                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       cAGGAAAGCTCTTCGCATAGCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTCCCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTa                                                                                                                                                                                                                                                 <  2:64238/251‑1 (MQ=255)   ggAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTAcc                                                                                                                                                                                                                                               >  1:133831/1‑251 (MQ=255)       aGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAgg                                                                                                                                                                                                                                            >  2:231857/1‑250 (MQ=255)              cGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTgg                                                                                                                                                                                                                                     <  2:157401/250‑1 (MQ=255)               gCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGcc                                                                                                                                                                                                                                   >  1:56209/1‑251 (MQ=255)                            aaaGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCg                                                                                                                                                                                                                       >  1:263973/1‑250 (MQ=255)                            aaaGTAAATGTAATACGCTGTTCTTAATTGTTCAGCGCTCACCACCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCg                                                                                                                                                                                                                       >  2:526147/1‑250 (MQ=255)                                          tcGCTGTTCTTAATTGTTCAGTTCTCACCCCCCAGAGTATATCTCATTCCATTATTCAACTTAGCTACATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCTGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACTTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGt                                                                                                                                                                                                         <  2:528841/249‑1 (MQ=255)                                               gTTCTTAATTGTTAAGTGCTCACCCCCAAGTGTATAGCTCATTCCATTATTCAACATAGCTCCATATAAAGGGACAGATAATGGAATCGCTCGCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATACATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGATGAATTGTCGAATTTGAACAAGGTGGCg                                                                                                                                                                                                     <  2:33685/249‑1 (MQ=255)                                                 tCTTAATTGTTCAGTGCTCCCCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCCCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAgg                                                                                                                                                                                                  <  2:443280/250‑1 (MQ=255)                                                    taagtggtCAGTTCTCACCCCCCTGAGTATAGCTCATTCCATGATTCAACTTAGCTCCAGATAAAGGGACAGATAATGGAAGCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCAGCTGGGCGATAGCGCCCGGACGGAAACGGTACGGAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTTGCGAGGcct                                                                                                                                                                                               <  2:575298/243‑2 (MQ=255)                                                     aaTTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAg                                                                                                                                                                                              <  1:98117/250‑1 (MQ=255)                                                                           agagTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTa                                                                                                                                                                        >  2:270928/1‑250 (MQ=255)                                                                              gTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGAGAGCGCCCGGACTGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCa                                                                                                                                                                    <  1:546176/251‑1 (MQ=255)                                                                                 taGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGAGTGCGTGTGGTGCGCAGCTAATCAAGc                                                                                                                                                                 >  2:395279/1‑251 (MQ=255)                                                                                      cATTCCATTATTCAACTTAGCTCCATATAAAGGGACATATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACt                                                                                                                                                             <  1:231857/250‑1 (MQ=255)                                                                                                  cAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTa                                                                                                                                                <  2:56209/251‑1 (MQ=255)                                                                                                         gCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTa                                                                                                                                          >  1:58963/1‑250 (MQ=255)                                                                                                                      ggACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGCTAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCAGACGAATTGTCGATTTTGAACAAGGAGGAAGGGCCGGAGCTGCGTATGGGGCGCAGCTAATCACGCGCCAATCAAAgatgggaagggcgccgtgtggggaa                                                                                                                              >  2:39427/1‑225 (MQ=255)                                                                                                                                tGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCAAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGc                                                                                                                   >  2:309019/1‑250 (MQ=255)                                                                                                                                   aaTCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCaaa                                                                                                                >  1:28750/1‑250 (MQ=255)                                                                                                                                      cGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCaaaaaa                                                                                                             <  2:263973/250‑1 (MQ=255)                                                                                                                                          ctcAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACtt                                                                                                          <  1:301875/249‑1 (MQ=255)                                                                                                                                           ccAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGTCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACg                                                                                                       <  1:397518/250‑1 (MQ=255)                                                                                                                                            cAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAAGCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACg                                                                                                       <  2:386974/250‑1 (MQ=255)                                                                                                                                               aatgAACCACAGCAGGCTACCAGCAATTCCACGACGTTCTTATTCATCTTGTCGATAGCGCCCGGACAGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGAGGCGAGGCCAGGGCTGCGTATGGAGCGCAGCTAATCAAGCGACTATCAAATGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCa                                                                                                    >  1:477735/4‑250 (MQ=255)                                                                                                                                                         agcatGCTACCAGCAAATCCACGACGGTATGATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGTCGAGGCCAGGGCGGCTTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGTTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCttt                                                                                         <  2:275271/251‑1 (MQ=255)                                                                                                                                                                 aCCAGCAAATCCACGACGGTATTATTCAGCTGGTCGAGAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTATAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACct                                                                                  <  2:116827/250‑1 (MQ=255)                                                                                                                                                                  ccAGCAAATCCACGCCGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGAGTTAACGCGTTAATGACCGCGACGTACTAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCACAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACctct                                                                                >  2:244675/1‑251 (MQ=255)                                                                                                                                                                                actgTATTATTCATCTGGTCGATAGCGACCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTGTTTCAACATGtt                                                                  <  2:387774/248‑1 (MQ=255)                                                                                                                                                                                     attattCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGATCACGACGTAACAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAAt                                                             <  1:147781/251‑1 (MQ=255)                                                                                                                                                                                     attattCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAAt                                                             <  2:530159/251‑1 (MQ=255)                                                                                                                                                                                     atgattCATCTGGTCGATAGCTCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGTCGAATTGTCGAATTTGAACATGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCTAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTATGCAAAAAACTTACGCCAAATGAGGCTTTTGTACCTCTTTTTTCAACATGTTGAAAt                                                             <  1:377403/248‑1 (MQ=255)                                                                                                                                                                                        attCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCTACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGATGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTc                                                           <  1:398506/250‑1 (MQ=255)                                                                                                                                                                                               gggTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACa                                                   <  2:148655/250‑1 (MQ=255)                                                                                                                                                                                                         cgcCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCg                                         <  2:133831/251‑1 (MQ=255)                                                                                                                                                                                                                    aaaCGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCa                              <  1:298096/251‑1 (MQ=255)                                                                                                                                                                                                                      aCGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACTCCAtt                            <  2:505709/251‑1 (MQ=255)                                                                                                                                                                                                                           aCGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGa                        >  1:103316/1‑250 (MQ=255)                                                                                                                                                                                                                                           gtAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAAGTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTGTTTACCTCTTTTTTCAACATGTTGAAAATCACCAGGCAATGTCTTGCGAATTAACACCATTGGGAAATTCGCGGACACCGTc       <  2:244899/250‑1 (MQ=255)                                                                                                                                                                                                                                                gACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACCCCGTCTGCa   <  2:488472/250‑1 (MQ=255)                                                                                                                                                                                                                                                gACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCa   <  1:295001/250‑1 (MQ=255)                                                                                                                                                                                                                                                gACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAg  <  1:395279/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       CAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAG  >  NC_000913/3367701‑3368190

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