Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,340	0	A	C	45.2%	104.9 / 18.5	76	G61G (GGT→GGG)	ettA	energy‑dependent translational throttle A
Reads supporting (aligned to +/- strand):  ref base A (26/14);  new base C (0/33);  total (26/50)
Fisher's exact test for biased strand distribution p-value = 5.55e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCG  >  NC_000913/4630194‑4630490                                                                                                                                                    |                                                                                                                                                       aGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCtt                                                                                                                                                    <  1:253128/151‑1 (MQ=255)  gcgcATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTc                                                                                                                                                   <  2:33496/151‑1 (MQ=255)     cATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGAt                                                                                                                                                <  1:371925/151‑1 (MQ=255)        acacTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTc                                                                                                                                             >  1:259143/1‑151 (MQ=255)          acTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCtt                                                                                                                                           <  1:163685/151‑1 (MQ=255)               aTCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATc                                                                                                                                       <  1:519916/150‑1 (MQ=255)                 ccAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCaa                                                                                                                                     <  2:780166/150‑1 (MQ=255)                    ggCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGcc                                                                                                                                 >  2:495141/1‑151 (MQ=255)                      cGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCg                                                                                                                               >  1:846204/1‑151 (MQ=255)                        tttCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGcc                                                                                                                             >  1:341453/1‑151 (MQ=255)                            aGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATa                                                                                                                          >  1:562315/1‑150 (MQ=255)                             gggCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAAt                                                                                                                        >  2:1005449/1‑151 (MQ=255)                             gggCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAAt                                                                                                                        >  2:716401/1‑151 (MQ=255)                                    aaCCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCgca                                                                                                                   >  1:1009416/1‑149 (MQ=255)                                        aCTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAggg                                                                                                             >  2:330432/1‑151 (MQ=255)                                          ttCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTg                                                                                                           >  2:901627/1‑151 (MQ=255)                                             tGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGAc                                                                                                        >  2:633001/1‑151 (MQ=255)                                               aaccctttccttcAAGGGTCTCACGCACGGTGTGTTCCGGGTTCAGGTGCGTTGCCGGCGGGAGATAACCAATTTT‑GATGCCTGCCTGCGGGCGCTTTTCGCTTGCGTTGTCTTTATCATGGCCCGCCATAATGCGCAGCAGGGt                                                                                                            <  2:531329/136‑1 (MQ=255)                                               aacccGTTTCTCCAATGGACCCACGCGGGGTGTTTTCGGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCGGGCTGCGGACGCGCTTCCCTTTCGTTGTTTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACtt                                                                                                      <  2:612004/148‑1 (MQ=255)                                               aacccGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCATTCTT‑GATGTCTGGCTGCGGACGCGCTTCCCCTTCGATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACtt                                                                                                      <  1:479989/148‑1 (MQ=255)                                               aaaCCGTTCCTTCAAGGGCCCCACTCATTTTGGTTTCCGGGTCCGCCTGGGTTTCCGGCGGCAGAAAACAACTTTT‑TATTGCCGCCTGCGCAGGCGGTTCCCCTTAGGTGTTTTTATCAGTGCCCGCCATAATGCGCAGCAGGGTGGACtt                                                                                                      <  2:676515/151‑1 (MQ=255)                                                      ttcttcAGTGGCCCCCCGCGCGGTGTGTCCGGGTTTCAGCTGCGGTTCCGGCGGCAGAAACCCAATCTT‑GATGTCGGGCTGCGGCCGCGTTTCCCCTTCGTTGTCTTTACAAATCCCCGCCATAATGCGCAGCAGGGTGGACTTAccc                                                                                                  <  2:28948/148‑1 (MQ=255)                                                       tcttcAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGCGGACTTACCCgcgc                                                                                              >  2:470719/1‑151 (MQ=255)                                                          tcAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCAt                                                                                           >  2:728682/1‑151 (MQ=255)                                                            aagggACCCACGCGGGGTGTGTGCCGGGTTCAGCTGCGTTTCCTGCGGCAGATAACCAATTTT‑GTTGTCTGGCGGCGGCCGCGCTCCCCCTTCGTTGTCTTTAACAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCAt                                                                                           <  2:886770/146‑1 (MQ=255)                                                            aaTGGACTCACGCATGGTGTGTTCCGGTTCCAGCGGCGGTTCCTGCGGAAGATACCAATTTTT‑GATGTCGGGCGGCGGACGCGCTTCCCCTTCGTTTTCTTAATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCAt                                                                                           <  2:755796/149‑1 (MQ=255)                                                               gtccccccgcgcgGTTTGTCCGGGGTTCGGCTGCGGTTCCTGCGGCAGATACCCAATTTT‑GATGTCGGGCTGCGGCCGCGCTTCCCCTTCGATGTTTTTATCAAGGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAg                                                                                       <  2:194443/143‑1 (MQ=255)                                                               ggccccacgccggGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAg                                                                                       <  2:314704/145‑1 (MQ=255)                                                               ggACCCCCGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAAAAACCCAATTTT‑GATGTCTGGCTGCGGACGGCCTTCCCTTTCGATTTCTTTATCAAGCCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAg                                                                                       <  2:368500/150‑1 (MQ=255)                                                                gACTCCCGCACGGTGTGTTCCGGGTTCAGCGGCGGTTCCTGCGGCAAATACCCATTCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATCCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCa                                                                                        <  2:221494/148‑1 (MQ=255)                                                                gACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGAc                                                                                     >  2:918163/1‑151 (MQ=255)                                                                  cccacgcccgGTGTGTTCCGGTTTCTGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCCCTTCACCTTCAATGTCTTTATCAGTGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGAccc                                                                                   <  1:887006/149‑1 (MQ=255)                                                                  cTCACGCCCGGTTTGTTCCGTGTTCAGCTGCGTTTCCTGCGGCAGAAAACCAATCTT‑GTTGTCTGGCTGCGGACGCGCTTCCCCTTCGTTGTCTTTATCATTGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGAcc                                                                                    <  2:660734/150‑1 (MQ=255)                                                                        cacgGTGTGTGCCGGGTTCAGCTGCGGTTCCTGCGGCAAATAACCAATCTT‑GTTGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGa                                                                               <  2:824388/149‑1 (MQ=255)                                                                            gtgtgtTCGGGGTTCAGCTGCGTTCCCGGCGGAAGATAACAATTCTT‑GATGTCGGGCGGCGGCCGCGCTTCCCCTGCGATGTCTTTAACAACGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACAcc                                                                           <  2:9158/149‑1 (MQ=255)                                                                                 tgccggatgacgCTGCGTTCCCTGCGCCAGAT‑ACCCATCTTGGGTGTCTGGCCGCGGCCGCGCTTCCCCTGTGATTTTTTTATCATTCTCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAAtttt                                                                     <  2:285991/139‑1 (MQ=255)                                                                                      ggTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGcccc                                                               >  1:917855/1‑151 (MQ=255)                                                                                          cAGCTGCGTTCCCTGCGGCAGATAACCAATTTT‑GGTGTCTGGCGGCGGACGCGCTTCTCCTTCGATGTCTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAg                                                             <  2:175978/149‑1 (MQ=255)                                                                                          cAGCTGCGGTTCCGGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCCCTTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAg                                                             <  2:587571/149‑1 (MQ=255)                                                                                                cgtttCCTGCGGAAGATAACCAATTTT‑GTTGTCTGCCTGGGGCCGCTCTTCCCCTTCGTTTTTTTTACAAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGgaaga                                                       <  2:490231/146‑1 (MQ=255)                                                                                                cGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGgaagaa                                                      >  2:207852/1‑150 (MQ=255)                                                                                                 ggTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAAc                                                    >  1:54230/1‑151 (MQ=255)                                                                                                 ggTTCCGGCGGCAGATAACCAATCTT‑GATGTCGGGCTGCGGACGCTCTTCCCCTTGGATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGgaaga                                                       <  1:440106/148‑1 (MQ=255)                                                                                                     ccTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTc                                                  >  2:639757/1‑149 (MQ=255)                                                                                                       tGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCag                                                >  1:236147/1‑149 (MQ=255)                                                                                                        gCGGCAGATAACCATTCTT‑GATGTCGGGCTGCGGCCGCGCTTCCCCTTCGATTTCTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagag                                              <  2:131518/150‑1 (MQ=255)                                                                                                         cGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGCCCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagag                                              >  2:952930/1‑149 (MQ=255)                                                                                                          ggaaGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCACCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagag                                              <  1:513728/145‑1 (MQ=255)                                                                                                          ggCAGATACCCATTTTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagagag                                            <  1:561610/150‑1 (MQ=255)                                                                                                          ggCAGATAACCAATTTT‑GATTTCGGGCGGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCag                                                <  2:283098/146‑1 (MQ=255)                                                                                                          ggCAAATAACCATTCTT‑GATGTCTGGCTGCGGACGCGTTTCCCTTTCGTTTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagagag                                            <  2:492544/150‑1 (MQ=255)                                                                                                          gacAGATACCCAATCTT‑GGTGTGGGGCTGGGGCCGGGCTTCCCCTTTGATTTTTTTATAATTGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCagagag                                            <  2:900358/148‑1 (MQ=255)                                                                                                             aGATAACCAATTTT‑GATGTCTGGCTGCGGACGCGCTTCACTTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATg                                         <  1:120657/150‑1 (MQ=255)                                                                                                               aTACCCACTCTT‑GGTGTCGGGCTGGGGCCGCGTTCCCCTTTCGTTTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                      <  2:71193/151‑1 (MQ=255)                                                                                                                 aaCCATTCTT‑GATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGtttt                                     <  1:158559/150‑1 (MQ=255)                                                                                                                 aaCAAATCTT‑GGTGTCTGGCGGCGGACGCGTTTCCCCTTCGTTTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttttt                                    <  2:972200/151‑1 (MQ=255)                                                                                                                       tCTT‑GATTTCTGGCTGCGGACGCGTTTCCCCTTCGATGTCTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCaaaa                               <  1:466065/150‑1 (MQ=255)                                                                                                                       tCTT‑GATTTCTGCTTGCGGCCGCGCTTCTCCTGCGCTGTCTTAACCAATTCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCaaaa                               <  1:495147/150‑1 (MQ=255)                                                                                                                       tCTT‑GATGTCTGGCTGCGGCCGCGCTTCCCCTTCGATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCaaaa                               <  2:375598/150‑1 (MQ=255)                                                                                                                       tCTT‑GATGTCTGGCTGCGGACGCGCTTCCCCTTCGTTTTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCaaa                                <  1:233538/149‑1 (MQ=255)                                                                                                                         tt‑GATGTCTGGCTGCGGACGCGCTTCCCCTTCGTTTTCTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGtt                                       <  2:990748/140‑1 (MQ=255)                                                                                                                         tt‑GATGTCTGGCTGCGGACGCGCTTCCCCTTCGATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAat                              <  2:701137/149‑1 (MQ=255)                                                                                                                          t‑GATGTCGGGCTGCGGACGCGCTTCCCCTTCGATGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAatat                            <  2:386955/150‑1 (MQ=255)                                                                                                                             aTGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc                         >  1:463877/1‑151 (MQ=255)                                                                                                                             aTGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc                         >  2:339832/1‑151 (MQ=255)                                                                                                                                 cTGGCGGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGAc                         <  1:242164/147‑1 (MQ=255)                                                                                                                                  tGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTc                    >  2:308067/1‑151 (MQ=255)                                                                                                                                   ggCTGCGGACGCGCTTCACCTTCGATTTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGttt                     <  1:932022/149‑1 (MQ=255)                                                                                                                                       gCGGACGCGCTTCCCCTTCGTTTTCTTTATCAATCCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcggcg                <  2:111667/150‑1 (MQ=255)                                                                                                                                       gCGGACGCGCTCCCCTTTGGATGTTTTTACCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcggcg                <  1:697808/150‑1 (MQ=255)                                                                                                                                         ggccccgcTTCCCCTTCGATTTTTTTATAAATCCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcggc                 <  2:662380/142‑1 (MQ=255)                                                                                                                                         ggACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcggcgg               <  2:397539/149‑1 (MQ=255)                                                                                                                                           aCGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATTCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGaac            >  1:505410/1‑150 (MQ=255)                                                                                                                                            cgcgCTTCCCCTTCGATGTTTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGaac            <  2:770524/149‑1 (MQ=255)                                                                                                                                             gcgcTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGaacaac         >  1:560847/1‑151 (MQ=255)                                                                                                                                                 ttCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTg     >  2:71086/1‑151 (MQ=255)                                                                                                                                                    aCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCg  >  1:480783/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       AGCGCATACACTTCATCCAGGCGTTTCAGGGCGTTAACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTT‑GATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCG  >  NC_000913/4630194‑4630490

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: ‑