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.5%	123.0 / 13.1	65	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (7/31);  new base G (27/0);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 3.74e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATAT  >  NC_000913/2434733‑2435032                                                                                                                                                       |                                                                                                                                                       cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGt                                                                                                                                                        <  1:806689/150‑1 (MQ=255)   ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAAt                                                                                                                                                     <  1:311448/151‑1 (MQ=255)      ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCGGGCAGGGGAATcc                                                                                                                                                   >  1:148404/1‑150 (MQ=255)      ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                                                  <  1:394510/151‑1 (MQ=255)      ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                                                  <  1:814863/151‑1 (MQ=255)        atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGCAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCGGGCAGGGGAATCCAcc                                                                                                                                                >  1:260097/1‑151 (MQ=255)          atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACcgc                                                                                                                                              <  2:803942/151‑1 (MQ=255)          atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGGAATCCAcccc                                                                                                                                              >  2:139299/1‑149 (MQ=255)            ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACcgcg                                                                                                                                             <  2:476424/150‑1 (MQ=255)            ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACcgcg                                                                                                                                             <  1:836949/150‑1 (MQ=255)            ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAAAAGCGGGCCCATTGGCGG‑TTTTTGAAGAGCATAACCGTCAAGGGAATCCAcccc                                                                                                                                              >  1:628496/1‑147 (MQ=255)             tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGa                                                                                                                                            <  1:933230/150‑1 (MQ=255)               ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGAcc                                                                                                                                          >  1:372467/1‑150 (MQ=255)                cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGGTCATACCGGTCGGGGGAAAACAcccccgcc                                                                                                                                          >  1:18711/1‑143 (MQ=255)                  aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                      <  2:186643/151‑1 (MQ=255)                          tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAAccc                                                                                                                               >  1:107393/1‑150 (MQ=255)                            cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATAACGGGCAGGGGAATCCACCCCGACCCGATACAGCCCtt                                                                                                                             >  1:77178/1‑150 (MQ=255)                              cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTcc                                                                                                                           <  1:885607/150‑1 (MQ=255)                                              gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGACCCGATACAGCCCCTCCCCTTTTTCCGTtgtt                                                                                                            >  1:360387/1‑147 (MQ=255)                                              gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgat                                                                                                            <  2:66495/149‑3 (MQ=255)                                                aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACCTTACCGG‑CAGGGGAAATCCCCCCGACCCGAAACAACCCTTCCCCTTTTTCCGTTGCTgcc                                                                                                         >  1:823009/1‑150 (MQ=255)                                                       ttGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAg                                                                                                    >  2:712185/1‑148 (MQ=255)                                                         ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGtt                                                                                                <  2:163474/150‑1 (MQ=255)                                                         ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGtt                                                                                                >  2:40642/1‑150 (MQ=255)                                                          gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATACCGGGCAGGGGAATACACCCCGCCCAGATACAACCCTTCCGCTTTTTTCGTTGCTgccgcc                                                                                                      >  1:701589/1‑143 (MQ=255)                                                                   tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCCGCCGCCCGCGGTCTGTCCttt                                                                                       >  2:882771/1‑149 (MQ=255)                                                                   tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑GTTTTGGAGATCATACCGGGCAGGGTAAACCACCCCGACCCGATACAGCCCTTCCCCTTTTGCCGGTGCTGCCCCCAGCGGTCTGTCCttt                                                                                       >  1:30681/1‑149 (MQ=255)                                                                    gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCCCTTTTTCCGTTGCTGCCCCCCGCGGTCTGTCCttt                                                                                       >  1:860250/1‑148 (MQ=255)                                                                    gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgc                                                                                     <  1:19973/150‑1 (MQ=255)                                                                     aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCTGGCAGGGGAATCCACCGCGACCAGAAACAACCCTTCCGCTTTTTGCGTTTGTGCCCCCCGCGTTCTGTCCTTTtcc                                                                                    >  1:527412/1‑150 (MQ=255)                                                                         gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAAAACATACCAGTCAGGGGAATCCACCCCCACCAGATACAGCCCTTCCACGTTTTCCGTGGCTGCCCCCCGCGTTCTTTTCCTTgccggc                                                                                 >  1:830599/1‑147 (MQ=255)                                                                         gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                                <  1:349540/150‑1 (MQ=255)                                                                          ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATACCGGGCCGGGGAAACCACCCCGGCCAGATAAAGCCCTTCCCCTTTTTTCGTTTCTTCCCCCAGCCGTCTTTTCTTTTCCCCCCGc                                                                              >  1:328490/1‑151 (MQ=255)                                                                            ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                                             <  2:256035/150‑1 (MQ=255)                                                                                  attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGc                                                                      <  1:195650/151‑1 (MQ=255)                                                                                     attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCCGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCCCCAGCGTTCTGTCCTTTGCCGCCCGCAActc                                                                         >  1:611844/1‑145 (MQ=255)                                                                                     attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATCCTAACGGTCCGGGGAAACCACCCCGGCCAGATACACACCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCCACTCTCCTAt                                                                   >  2:18239/1‑151 (MQ=255)                                                                                      ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTAt                                                                   <  1:283323/150‑1 (MQ=255)                                                                                         tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATACCGGGCAGGGGAAAACACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGGTGGTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATcc                                                                 >  1:404949/1‑149 (MQ=255)                                                                                          cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCa                                                                <  2:117784/149‑1 (MQ=255)                                                                                            ttAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATTATACCGGGCAGGGGAATCCAACCCGACCAGATAAAGCCCTTCCGCTTTTTCCGTTGTTGCCCCCCGCGTTCTGTCCTTTTTCCCCCGCAACTCTGCTCTcc                                                                 >  1:628693/1‑146 (MQ=255)                                                                                            ttAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCAGCGGTTTGTCCTTTGCCGCCCGCAACTCTGCTAttc                                                                 >  2:787006/1‑144 (MQ=255)                                                                                                 ccGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                                        <  1:38188/150‑1 (MQ=255)                                                                                                 ccGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                                        <  1:475473/150‑1 (MQ=255)                                                                                                  cGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                                      <  2:933375/151‑1 (MQ=255)                                                                                                    ccAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                                                    <  2:237618/151‑1 (MQ=255)                                                                                                      aGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAACCCTTCCCCTTTTTTCGTTGCTGCCCCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGGTTTCCAGCTCTCCCGCCgg                                                   >  2:269352/1‑150 (MQ=255)                                                                                                        aaaTAGCGGGCCCATTGGCGG‑TTTTGGGAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGCTCTCCCGCTgg                                                   >  2:387291/1‑148 (MQ=255)                                                                                                        aaaTAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtg                                                <  2:803679/151‑1 (MQ=255)                                                                                                                 gCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCCCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCCGCTCTCCGGCTGGGTGTGGGGAcc                                        >  1:563836/1‑150 (MQ=255)                                                                                                                  cccATTGGCGG‑TTTTTGAAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCCGCCGCCCGCGTTTTGTCCTTTGCCGCCAGCAACTCTGCTTTCCAGCTCTCCGGCCGGGTGTGGGCaac                                        >  1:164046/1‑147 (MQ=255)                                                                                                                            g‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                              >  2:615097/1‑149 (MQ=255)                                                                                                                               tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                           <  2:87167/150‑1 (MQ=255)                                                                                                                               tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                           <  2:473731/150‑1 (MQ=255)                                                                                                                                ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc                          <  2:679515/150‑1 (MQ=255)                                                                                                                                      gATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGAc                   <  1:823784/151‑1 (MQ=255)                                                                                                                                       aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                  <  1:805350/151‑1 (MQ=255)                                                                                                                                        tCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                  >  1:570188/1‑150 (MQ=255)                                                                                                                                          aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat                >  1:600259/1‑150 (MQ=255)                                                                                                                                          aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat                >  1:667337/1‑150 (MQ=255)                                                                                                                                            aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt             <  1:135533/151‑1 (MQ=255)                                                                                                                                               ggTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt          >  1:729745/1‑151 (MQ=255)                                                                                                                                                  cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGa        <  2:715709/150‑1 (MQ=255)                                                                                                                                                   aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc      >  1:37847/1‑151 (MQ=255)                                                                                                                                                    gggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCa     <  2:5520/151‑1 (MQ=255)                                                                                                                                                     ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat    <  2:107395/151‑1 (MQ=255)                                                                                                                                                       tAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCatat  >  2:57049/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATAT  >  NC_000913/2434733‑2435032

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