Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,882	0	T	G	41.6%	129.0 / 34.8	89	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (18/34);  new base G (37/0);  total (55/34)
Fisher's exact test for biased strand distribution p-value = 1.03e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434739‑2435031                                                                                                                                                |                                                                                                                                                       cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATTCa                                                                                                                                                  >  1:994132/1‑150 (MQ=255)  atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAACCCGTCAGGGGAAaccacc                                                                                                                                                >  1:120619/1‑151 (MQ=255)  atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGATCATAACCGGCAGGGGGATTCAcc                                                                                                                                                >  2:1163266/1‑151 (MQ=255)    atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGGTCTTAACGGTCCGGGGGATCcccccc                                                                                                                                              >  1:282106/1‑149 (MQ=255)    atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc                                                                                                                                              <  2:75186/151‑1 (MQ=255)      ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCCGGGGAAaccacccc                                                                                                                                              >  1:214928/1‑147 (MQ=255)      ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCGGGCCGGGGAAAACACcg                                                                                                                                               >  1:930457/1‑148 (MQ=255)         ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCa                                                                                                                                         <  2:202134/151‑1 (MQ=255)          cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTAAACCGGGCAGGGGAATCCAccc                                                                                                                                               >  1:332668/1‑143 (MQ=255)           tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATAACGGGCCGGGGaaaacaacccgac                                                                                                                                           >  2:682347/1‑136 (MQ=255)           tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAGATCAAATCGGGCCGGGGAATCCACCGCGAcc                                                                                                                                          >  2:902936/1‑148 (MQ=255)            aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATAACGGTCAGGGGAATCCACCCCGccc                                                                                                                                          >  2:1101290/1‑147 (MQ=255)            aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                      <  2:233178/151‑1 (MQ=255)              gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTTGAGAACATTACTGTCAGGGGAAAACAACGCGACCAg                                                                                                                                        >  1:473558/1‑147 (MQ=255)                 ttGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCGCGGCCCGATACaac                                                                                                                                 >  2:52275/1‑149 (MQ=255)                 ttGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCCGGGGAAACCACCCCGACCAGATACa                                                                                                                                   >  2:659777/1‑149 (MQ=255)                    tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAAACCGGCAGGGGAATCCACCCCGAccccaaacagacc                                                                                                                               >  1:785035/1‑139 (MQ=255)                     aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCt                                                                                                                              <  2:71524/150‑1 (MQ=255)                      cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCCTACCGGGCAGGGTGATCCACCCCGACCCCCTACACCCCTTc                                                                                                                            >  1:379509/1‑151 (MQ=255)                      cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATTCACCCCGGCCAGATACAAccc                                                                                                                               >  2:239540/1‑148 (MQ=255)                          tCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGTCAGGGGAATCCACCGCGACCCGAAACAGCCCCTcccc                                                                                                                         >  1:1083987/1‑148 (MQ=255)                               tttCAGGCTGCAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCAGAAACAGCCCTTCCCCtttt                                                                                                                     >  2:838154/1‑149 (MQ=255)                                ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAcc                                                                                                                                                <  2:705738/121‑1 (MQ=255)                                 tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGGCGGGGAATCCACCCCGACCCGATACAACCCCTCCGCtttt                                                                                                                     >  1:227221/1‑147 (MQ=255)                                        gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgctg                                                                                                           <  2:49747/150‑1 (MQ=255)                                           aCTACGTTGGCTACAGCCTGAAATGTTCCGCGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgccgc                                                                                                       <  2:738420/151‑1 (MQ=255)                                              aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGCACCGATAAAACCCTTCCCCTTTTTCCGTTTCTgccgcc                                                                                                      >  2:1237675/1‑149 (MQ=255)                                              aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCa                                                                                                     <  1:959827/150‑1 (MQ=255)                                              aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAg                                                                                                    <  2:695936/151‑1 (MQ=255)                                               cGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGc                                                                                                   <  1:599859/151‑1 (MQ=255)                                                gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATAACGGGCAGGGGAAATCACCCCGCCCAGATAAAGCCCCTTCCCTTTTTCCGTTGCTTCCCCCcgcg                                                                                                  >  1:1188466/1‑151 (MQ=255)                                                gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGc                                                                                                   <  2:97101/150‑1 (MQ=255)                                                gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCg                                                                                                  <  2:1255653/151‑1 (MQ=255)                                                       aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGTAATCATAACCGTCGGGGGTATCCCCCCCGACCAACTTAACCCCCTTCGCTTTTTCCCTTGCTTCTGCCCGCGGTTTGt                                                                                            >  2:514533/1‑150 (MQ=255)                                                         aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCt                                                                                         <  2:1094303/151‑1 (MQ=255)                                                            cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAAATCACCCCGGCCAGAAAAAACCCCTCCCCCTTTTCCGGTTCTgccgcgaaccggagagcac                                                                                         >  2:639252/1‑134 (MQ=255)                                                            cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgccg                                                                                                        <  1:639233/133‑1 (MQ=255)                                                               aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGTGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccg                                                                                   <  2:1098613/151‑1 (MQ=255)                                                                      ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAGCAAAACGGGCAGGGGAAACCACCGCGACCAGATACAACCCCTCCGCTTTTTCCGTTGCTGCCGCCCGCGGTCTGTCCTTTtccgcc                                                                                 >  2:854490/1‑146 (MQ=255)                                                                        gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCAGGGGAATCCACCGCCACCAGAAACAACCCTTCCGCCTTTTCCGTTGCTGCCCCCCGCGTTTTGTCCTTTGCCGCCCGCAA‑cc                                                                         >  2:833375/1‑150 (MQ=255)                                                                          gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAGATCCTACCCGTCAGGGTAATCCACCGCGACCCAATACAGCCCTTCCCCCTTTTCCCTTTCTGCCGCCCACGTTCTGTCCCTTTCCGCCCGCaaaac                                                                         >  1:439754/1‑147 (MQ=255)                                                                           cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCAAACCGGGCCGGGGAATACACCGCGGCCAGATACAGGCCTTCCGCTTTTGCCCTTGCTGCCGCCCGCGGTTTGTCCTTTGCCGCCCGCAA‑CTCTg                                                                      >  2:269530/1‑151 (MQ=255)                                                                            attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGTCAGGGGAATCCACCCCGACCCGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTTTGTCCTTTGCCGCCCGCCA‑ccc                                                                        >  2:1144745/1‑146 (MQ=255)                                                                            attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAATATACCCGTCAGGGGAATCCCCCGCGCCCAGATACAACCCCTCCGCCTTTTCCGTTTCTGCCGCCCCCGTTCTGTCCTTTtcc                                                                                    >  1:877748/1‑137 (MQ=255)                                                                                 tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAA‑CCCTGCTATcc                                                                >  2:901555/1‑151 (MQ=255)                                                                                  ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGTAGATTATACCCGGCGCGGGAATACCCCCCGACCCGATACAACCCTTCCGCCTTTTTCGGTGTTTCCCCCACCGTTCTTTCCTCTTTCGCCCGCAA‑CTCTTCTATcc                                                                >  1:969580/1‑150 (MQ=255)                                                                                   tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTTCGCCCGCCA‑CTCTTCTAt                                                                  >  2:575986/1‑147 (MQ=255)                                                                                       tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGAACATACCGGTTAGGGGGATACACCCCGACCACATACCGACCTTCCGCTTTTTCCGTTGGTGCCCCCCGCGGTCTGTCCTTTTGCGCCCGCCA‑ATCTTCTAttccgc                                                             >  1:938886/1‑143 (MQ=255)                                                                                       tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAAACGGTCAGGGTAAACCACCCCGGCCAGATACAACCCTTCCGCTTTTTCCGGTGCTGCCCCCCGCGGTCTGTCCCTTTCCCCCAGCAA‑CCCCTCTATCCAGccc                                                           >  2:775541/1‑148 (MQ=255)                                                                                       tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGTTAGGGTAATCCACCCCGACCCGATACAGCCCCTCCGCTTTTGGCGTTGCTGCCCCCAGCGGTCTGTCCTTTGCCCCCCGCAA‑CTCTGCTTTCCa                                                               >  2:1019248/1‑146 (MQ=255)                                                                                       tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGTAGAACATAACGGGGAGGGGAATTCCCCCCGACCAGATACAGCCCCTCCGCCTTTTGCGTTGCTGCCCCCCGCGGTCTGGCCTTTTCCGCCCGCAA‑CCCCGCTTTCCAGct                                                            >  2:472445/1‑149 (MQ=255)                                                                                          tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAAACCACCCCGACCCGAAAAAACCCCTCCCCTTTTTCCCTTGCTGCCGCCCGCGGTCTGTTCTTTTCCCCCAGCAA‑CTCTGCTATCCAGCTCTcc                                                        >  1:1150116/1‑150 (MQ=255)                                                                                           ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCACACCCGGCAGGGGAATTCACCCCGGCCCGATAAAAACCCTCCGCTTTTTTCGCTTCTGCCGCCCGCGTTTTGGCCCTTTTCCCCCGCAA‑CTCTGCTTTCCCGctct                                                          >  1:1224907/1‑147 (MQ=255)                                                                                             gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCt                                                    <  1:1113890/151‑1 (MQ=255)                                                                                               cAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCCCGGCCAGATACAACCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑ATCTGCTATCCAGCTCTCCCGCTgg                                                  >  2:999153/1‑151 (MQ=255)                                                                                                aGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTTCTGCCGCCAGCGTTCTGTCCTTTTCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTgg                                                  >  2:2336/1‑150 (MQ=255)                                                                                                 gAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAAACCCCCCCGGCCCGATAAAGCCCCTCCGCTTTTTCCGCTGCTGCCGCCCCCGTTCTGTCCTTTTCCGCCCGCAA‑CTCTGCTATCCAGCCCTCCGGCTgg                                                  >  1:67661/1‑149 (MQ=255)                                                                                                    aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTtgtg                                             <  1:47007/151‑1 (MQ=255)                                                                                                             ccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc                                    <  2:263824/151‑1 (MQ=255)                                                                                                                  ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                                <  2:760427/150‑1 (MQ=255)                                                                                                                    cGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                             >  2:833501/1‑151 (MQ=255)                                                                                                                     ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa                            <  2:235156/151‑1 (MQ=255)                                                                                                                      gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                           <  2:422636/151‑1 (MQ=255)                                                                                                                         ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc                         <  1:641217/150‑1 (MQ=255)                                                                                                                         ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc                         <  1:891594/150‑1 (MQ=255)                                                                                                                           ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                       <  2:973227/150‑1 (MQ=255)                                                                                                                            gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                     >  1:359533/1‑151 (MQ=255)                                                                                                                            gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGcc                     >  1:233495/1‑151 (MQ=255)                                                                                                                              aGATCATACCGGTCAGGGTAATCCCCCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                   <  1:1307586/151‑1 (MQ=255)                                                                                                                              aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                    >  1:1007192/1‑150 (MQ=255)                                                                                                                              aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                   <  2:434495/151‑1 (MQ=255)                                                                                                                              aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                   <  1:138594/151‑1 (MQ=255)                                                                                                                              aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                   <  1:280514/151‑1 (MQ=255)                                                                                                                                  cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat               <  2:235063/151‑1 (MQ=255)                                                                                                                                   aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat               >  2:11885/1‑150 (MQ=255)                                                                                                                                     aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt            >  1:194673/1‑151 (MQ=255)                                                                                                                                     aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt            >  1:310024/1‑151 (MQ=255)                                                                                                                                        ggTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt         >  1:141219/1‑151 (MQ=255)                                                                                                                                         gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg        <  2:478073/151‑1 (MQ=255)                                                                                                                                          tCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGCCTTCCATCAGGCTGCCGACAATATTCCTg        >  2:743780/1‑150 (MQ=255)                                                                                                                                          tCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg        >  1:24289/1‑150 (MQ=255)                                                                                                                                           cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc      >  2:1205102/1‑151 (MQ=255)                                                                                                                                           cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc      >  1:60265/1‑151 (MQ=255)                                                                                                                                             gggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc     >  2:456372/1‑150 (MQ=255)                                                                                                                                             gggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc      <  2:955304/149‑1 (MQ=255)                                                                                                                                              ggTAATCCCCCGGAACCGAATCCAGCCCTCCCGTTTTGGCCTTTCCTGCCGCC‑GCCTTCTTGCCCTTGCCCCCCACAACCTCTGCTACCCAGCTCTCCGGCTGGTTGT‑GGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  1:52275/151‑1 (MQ=255)                                                                                                                                               gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  2:502077/151‑1 (MQ=255)                                                                                                                                               gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  1:184733/151‑1 (MQ=255)                                                                                                                                               gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat   >  1:836479/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                       CATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAA‑CTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434739‑2435031

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