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

CCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGA  >  NC_000913/2434738‑2435026                                                                                                                                                  |                                                                                                                                                    ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATcc                                                                                                                                                <  2:935008/150‑1 (MQ=255) ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATcc                                                                                                                                                <  1:1126829/150‑1 (MQ=255)  cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATAACCGGCAGGGGAATcc                                                                                                                                                >  1:279469/1‑149 (MQ=255)   atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCAcc                                                                                                                                             <  2:879105/151‑1 (MQ=255)       ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑ACAAACCCGTCAGGGGAATCCACcgc                                                                                                                                           >  2:911597/1‑149 (MQ=255)          ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATAACCGGCAGGGGAAACCACCCCGAcc                                                                                                                                       >  1:643556/1‑150 (MQ=255)                    gTACTCGTCATCTTTCAGGCTGTAACTACGTCGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCGGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCCCCGCGACCAGATAC‑AG‑ccc                                                                                                                          <  2:524387/151‑1 (MQ=255)                     tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAAA‑TCATACCGGGCAGGGGAATCCCCCGCGGCCCGATAC‑aagcc                                                                                                                           >  1:965911/1‑146 (MQ=255)                         cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGTGCCCATTGGCGGTTTTGGAAGG‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTGcc                                                                                                                      <  2:420072/150‑1 (MQ=255)                         cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑AAATACCGGTCCGGGGAAAACACCACGCCCCGATAC‑AG‑CCCtt                                                                                                                        >  2:1140604/1‑148 (MQ=255)                             aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TCATACCCGGCCGGGGAATCCACCCCGACCCAATAA‑AA‑CCCCTcc                                                                                                                      >  2:847791/1‑146 (MQ=255)                                 ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTg                                                                                                              <  1:700706/149‑1 (MQ=255)                                  tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TCATAACGGGCAGGGGAAATCACCGCGACCCGGTAC‑AA‑CCCCTCCC‑CCTTTGcc                                                                                                            >  1:66408/1‑150 (MQ=255)                                  tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TCATAACCGGCAGGGGAATTCCCCCCGGCCCGATAC‑CG‑TCCCTCCC‑Ctttt                                                                                                               >  2:989819/1‑147 (MQ=255)                                         gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTtgctgc                                                                                                    <  1:1273448/151‑1 (MQ=255)                                         gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTtgctg                                                                                                     <  2:585377/150‑1 (MQ=255)                                           aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATACCGGGCAGGGGAATCCACCGCGACCAGATAC‑AG‑CCCTTCCC‑CTTTTGCCGTTGCTgcc                                                                                                   >  1:177827/1‑150 (MQ=255)                                              tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATACCGGGCAGGGGAAACCCCCCCGACCCGATAC‑AG‑CCCTTCCG‑CTTTTTTCGTTGCTtccccc                                                                                                >  2:582924/1‑150 (MQ=255)                                                   tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑ACATACCGGGCCGGGGAATCCACCCCGACCAGATACAAA‑CCCTCCCC‑CTTTTTCC‑TTGCTgccccccccgt                                                                                           >  1:930580/1‑145 (MQ=255)                                                    ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATACCGGGCAGGGGAATCCACCGCGACCAGATAC‑AA‑CCCCTCCG‑CTTTTTCCGTTGCTGCCCCCcgcggt                                                                                          >  1:728825/1‑148 (MQ=255)                                                    ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGA‑TCATACCGGGCAGGGGAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTTCCGTTGCTGCCGCCAGCggt                                                                                          >  1:958254/1‑148 (MQ=255)                                                    ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGtt                                                                                          <  1:1267198/150‑1 (MQ=255)                                                    ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAAATTCAAACC‑GCCAGGGGAAACCAACGCGACCCGAAAA‑AA‑CCCCTCCC‑CTTTTTCCGGTGCTTccgcccgc                                                                                             >  1:1107181/1‑147 (MQ=255)                                                        aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTg                                                                                       <  2:138102/149‑1 (MQ=255)                                                        aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGt                                                                                      <  1:640787/150‑1 (MQ=255)                                                         cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑ACATACCGGGCAGGGGAAACCACCCCGACCAGATAC‑AG‑CCCTTCCC‑CTTTTTTCGTTGCTGCCGCCCGCGtt                                                                                          >  2:774501/1‑145 (MQ=255)                                                          aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTGGAGA‑TCATAACCGGCCGGGGAATACACCCCCACCCGATAC‑AA‑CCCTTCCC‑CTTTTTCCGTTGCTGCCGCCCGCGGTCttgcc                                                                                    >  2:612417/1‑146 (MQ=255)                                                            ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                 <  1:513805/151‑1 (MQ=255)                                                              tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TCATACCGGTCAGGGGAAACCACCCCCACCAGATAA‑AG‑CCCCTCCC‑CTTTTTCCGTTGCTGCCCCCCGCGTTCTTTCCttt                                                                                 >  2:772194/1‑149 (MQ=255)                                                               gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTg                                                                                <  1:911155/149‑1 (MQ=255)                                                                aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccg                                                                             <  2:579591/151‑1 (MQ=255)                                                                 aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccgc                                                                            <  2:822117/151‑1 (MQ=255)                                                                  aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGG‑TCATATCCGGCCGGGGAAAACACCGCGACCAAATAA‑CA‑CCCCTCCG‑CCTTTTTCCTTGCTGGCGCCAGCGTTCTTTCCtttttccc                                                                            >  1:861872/1‑145 (MQ=255)                                                                  aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATACCGGTCAGGGGAAACCACCGCGACCAGATAC‑AG‑CCCCTCCC‑CCTTTTCCGTTGCTGCCCCCCCCGGTCTGTCCTTTtccgcc                                                                           >  2:543498/1‑151 (MQ=255)                                                                   tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TCATACCGGGCAGGGGAATCCCCCCCGACCAGAAAC‑AA‑CCCCTCCC‑CTTTTTCCGTTGCTGCCGCCCGCGTTCTTTCCTTTgccccc                                                                           >  1:976663/1‑150 (MQ=255)                                                                    gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TTATACCCGGCAGGGGAATCCACCGCGACCAGAAAC‑AG‑CCCCTCCC‑CTTTTTCCCTTGCTGCCCCCAGCGTTCTGTCCTTTgccccc                                                                           >  2:1055720/1‑149 (MQ=255)                                                                      tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGc                                                                        <  1:493491/150‑1 (MQ=255)                                                                       ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCaa                                                                      <  2:807147/151‑1 (MQ=255)                                                                         gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAAct                                                                    <  1:224239/151‑1 (MQ=255)                                                                         gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATAACGGTCAGGGGAAACCACCGCGACCAGAAAC‑AG‑CCCCTCCC‑CCTTTTCCGTTGCTGCCCCCAGCGGTCTGTCCTTTTCCGCCAGCaa                                                                      >  1:372997/1‑149 (MQ=255)                                                                            cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTg                                                                 <  1:234299/151‑1 (MQ=255)                                                                             attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑ACATACCCGGCAGGGGAATCCACCGCCACCAGAGAC‑AG‑CCCTTCCC‑CTTTTTCCGCTTCGTCCGCCCGCGTTTTTTCCTTTTCCCCCCGCAAatc                                                                   >  2:305074/1‑148 (MQ=255)                                                                             attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGA‑TCATACCGGTCAGGGGAATCCACCGCGACCCGATAC‑AG‑CCCCTCCG‑CTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCCTTTCCGCCCGCAAc                                                                     >  1:679778/1‑146 (MQ=255)                                                                             attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTTGAGC‑TGATACCGGGCGAGGGAATACACCCCCCACCGCTTC‑AA‑TCCCTCCG‑CTTTTGCCGTTTCCTCCGCCCGCGCTCCGTCCttttttgcccgcaa                                                                      >  2:409270/1‑134 (MQ=255)                                                                               tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCt                                                               <  1:1115470/150‑1 (MQ=255)                                                                                     cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCa                                                          <  2:279462/149‑1 (MQ=255)                                                                                      gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TCATACCGGGCAGGGGAATCCACCCCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTTCCGTTGCTGCCGCCCGCGTTCTTTCCTTTTCCGCCCGCAACTCTTCTATCCAGc                                                        >  1:556850/1‑150 (MQ=255)                                                                                      gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGct                                                       <  1:134686/151‑1 (MQ=255)                                                                                          gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctct                                                     <  2:550069/149‑1 (MQ=255)                                                                                           tCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATAACGGGCAGGGGAATCCACCCCGGCCAGATAC‑AG‑CCCCTCCC‑CTTTTGGCGTTTGTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAAATCTGCTATTCAGctcttc                                                   >  1:637407/1‑148 (MQ=255)                                                                                             cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                                <  1:386018/151‑1 (MQ=255)                                                                                              gCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGA‑TCATAACCGGCAGGGGAATCCACCGCGGCCCGATAC‑AG‑CCCCTCCG‑CTTTTTCCGGTGCTGCCGCCAGCGTTCTTTCCTTTTCCGCCCGCCACTCTGCTCTTCAACTcccccg                                                 >  1:554436/1‑147 (MQ=255)                                                                                               ccAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGA‑TCATAACCGTCCGGGGAAACCACCCCGACCCGATAC‑AA‑CCCCTCCGCCTTTTGCC‑TTTCTGCCGCCCGCGGTCTGTCCTTTTCCGACAGCAAATCTGCTATCCCGCTCTcccg                                                 >  1:701676/1‑146 (MQ=255)                                                                                                cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGA‑TCATAACGGTCAGGGGAAACCCCCCCGGCCAGATAA‑AG‑CCCTTCCC‑CTTTTTCCGTTGGTGCCGCCCGCGTTCTGTCCTTTGCCGCCCCCAAACCTGCTATCCAACTCTcc                                                   >  2:236622/1‑145 (MQ=255)                                                                                                 aGAAATAGCGGGCCCATTGGCGGTTTTGGAAAA‑ACATACCGGTCCGGGGAAACCACCCCGACCCGATAC‑AA‑CCCTTCCC‑CCTTTTTCGTTGTTGCCGCCACCGTTTTGTTCTTTGTCCCCAGCAACCCTGCTTTCCAGCCCCCCCGCTg                                              >  2:851614/1‑149 (MQ=255)                                                                                                   aaaTAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtg                                          <  2:713889/151‑1 (MQ=255)                                                                                                     aTAGCGGGCCCATTGGCGGTTTTGGGAGG‑TCATACCCGGCAAGGTAATACCCCCCGACCCAATAC‑CG‑CCCTTCCG‑CTTTTTCCCTTGCTGCCCCCACCGGTCTGTCCTTTTCCCCCCCCCACTCTTCTCTCCCACTCCCCCGCTgg                                             >  1:216510/1‑146 (MQ=255)                                                                                                         cGGGCCCATTGGCGGTTTTGGAAGA‑TCATAACGGGCAGGGGTATCCACCCCGACCAGAAAA‑AGCCCCTTCCC‑CTTTT‑CCGCTGGCGCCCCCCGCGGTCTGTTCCTTTCCGCCCGCAAATCTTCTTTCCAGCTCTCCCGCTTGGtgtg                                        >  1:1062073/1‑147 (MQ=255)                                                                                                             cccATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGa                                <  1:371124/151‑1 (MQ=255)                                                                                                                   ggcggcTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑CG‑CCCTTCCG‑CTTTTGCCGGTGCTGCCGCCAGCGTTCTGTCCTTTACCGCCAGCAACTCTGCTATCCCGCTCTCCGGCTCGTTGTGGGCACCGACATc                            >  2:740976/1‑149 (MQ=255)                                                                                                                    gCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                         <  1:1212740/151‑1 (MQ=255)                                                                                                                     cGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                        <  2:965842/151‑1 (MQ=255)                                                                                                                        ttttGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                     >  2:297437/1‑151 (MQ=255)                                                                                                                          ttGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc                    <  1:1237423/150‑1 (MQ=255)                                                                                                                          ttGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                   >  2:733907/1‑151 (MQ=255)                                                                                                                           tGGAAGA‑TCATACCGGTCAGGGTAATCCCCCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                   <  2:800685/150‑1 (MQ=255)                                                                                                                           tGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                   <  1:1247626/150‑1 (MQ=255)                                                                                                                           tGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                   <  1:726383/150‑1 (MQ=255)                                                                                                                              aaGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg               >  1:1138891/1‑151 (MQ=255)                                                                                                                                 a‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa            <  1:1041743/151‑1 (MQ=255)                                                                                                                                       aCCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAata         <  2:414219/149‑1 (MQ=255)                                                                                                                                        ccGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGCCAATAtt       <  2:181289/150‑1 (MQ=255)                                                                                                                                        ccGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc      >  1:771327/1‑151 (MQ=255)                                                                                                                                         cGGTCAGGGTAATCCACCGCGACCAGATTC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc     <  1:384003/151‑1 (MQ=255)                                                                                                                                          ggTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt    >  2:72583/1‑151 (MQ=255)                                                                                                                                          ggTCAGGGCAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt    >  1:928968/1‑151 (MQ=255)                                                                                                                                             cAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTCCTCCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCTGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCTACAATATTCCTGa  <  2:129056/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGA‑TCATACCGGTCAGGGTAATCCACCGCGACCAGATAC‑AG‑CCCTTCCG‑CTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGA  >  NC_000913/2434738‑2435026

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