Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG50.7% 100.3 / 130.0 72Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (36/0);  total (43/29)
Fisher's exact test for biased strand distribution p-value = 2.37e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-13

CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  NC_000913/2434733‑2435028
                                                                                                                                                       |                                                                                                                                                    
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGt                                                                                                                                                      <  2:951734/150‑1 (MQ=255)
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTa                                                                                                                                                     <  2:86925/151‑1 (MQ=255)
  ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑ACGGTCCGGGGaa                                                                                                                                                    >  2:1278192/1‑150 (MQ=255)
   gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATc                                                                                                                                                  <  1:955876/151‑1 (MQ=255)
    tCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATCATA‑ACGGTCCGGGGAAAcc                                                                                                                                                 >  1:365424/1‑151 (MQ=255)
       atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATA‑CCGGGCAGGGGAATCCAcc                                                                                                                                              >  1:213292/1‑151 (MQ=255)
       atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCAcc                                                                                                                                              <  2:439372/151‑1 (MQ=255)
       atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCAcc                                                                                                                                              <  2:1079874/151‑1 (MQ=255)
             gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATCATA‑CCGGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  1:888804/1‑151 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCCGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  1:304629/1‑147 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGa                                                                                                                                     <  2:371826/150‑1 (MQ=255)
                    ccTTGTACTCGTCATCTTTCAGGTTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑a                                                                                                                                <  1:1328923/151‑1 (MQ=255)
                       tgaaCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGTAGATCATA‑CCGGGCAGGGGAATCCACCCGGACCAGAAAC‑A‑Acc                                                                                                                            >  1:64856/4‑151 (MQ=255)
                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGGAA‑ATCATA‑ACCGTCAAGGGAATCCACCCCCACCAGAAAC‑A‑ccccc                                                                                                                          >  1:849193/1‑150 (MQ=255)
                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCGGGCAGGGGAATCCACCCCGACCAGAAAC‑A‑Gcccc                                                                                                                          >  1:510791/1‑150 (MQ=255)
                            tCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCGGGCCGGGGAAACCACCCCGACCAGATAC‑A‑ACCCCTcc                                                                                                                       >  2:969420/1‑151 (MQ=255)
                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCAGATAC‑A‑GCCCCTcccc                                                                                                                     >  1:411881/1‑149 (MQ=255)
                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGc                                                                                                                     <  1:158475/151‑1 (MQ=255)
                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGGCAGGGGAAACCACCGCGACCCGATAC‑A‑ACCCCTCCGc                                                                                                                     >  1:303544/1‑151 (MQ=255)
                              gTCATCTTTCAGGCTGTAACGACGTTGGCTACAGAATGACATGTTCCGAGCATTATTCTCGTAAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTAGAAGATCATA‑CCGGTCAGGGTAAAACACCGCGACCAGAGAC‑A‑GCCCTTCCGc                                                                                                                     <  2:755100/151‑1 (MQ=255)
                                cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCt                                                                                                                    <  2:166069/150‑1 (MQ=255)
                                 aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATTATA‑ACGGGCAGGGGAATCCACCGCGGCCCAATAC‑A‑ACCCCTCCGCttt                                                                                                                  >  1:502733/1‑151 (MQ=255)
                                     ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCttt                                                                                                                  <  1:492742/147‑1 (MQ=255)
                                      tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCGGGCAGGGGAATCCACCCCGGCCCGATACAA‑CCCCTTCCCCTTTTccc                                                                                                              >  2:883628/1‑151 (MQ=255)
                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGtt                                                                                                           <  2:102281/151‑1 (MQ=255)
                                           cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATAATA‑CCGGGCAGGGGAATACACCGCGACCCGATAA‑A‑ACCCTTCCCCTTTTTCCGTtgct                                                                                                        >  2:1338391/1‑151 (MQ=255)
                                            tGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGGAATTCACCGCGACCAGATAC‑A‑GCCCTTCCCCTTTTTCCGCtgccc                                                                                                       >  1:1083127/1‑149 (MQ=255)
                                             gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTtgctg                                                                                                       <  1:1334507/150‑1 (MQ=255)
                                             gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGAAAC‑A‑GCCCTTCCGCTTTTGCCGTtgctgc                                                                                                      >  1:1131285/1‑151 (MQ=255)
                                              tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATCATA‑CCGGTCAGGGGAATCCACCCCGACCAGAAAA‑A‑GCCCCTCCCCTTTTTCCGGTGCTgcc                                                                                                     >  1:305332/1‑151 (MQ=255)
                                                     gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCCGATAC‑A‑GCCCTTCCGCTTTTTCCGTTGCTGCCGCCcgcg                                                                                              >  1:974284/1‑151 (MQ=255)
                                                        ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGtt                                                                                            <  2:693894/150‑1 (MQ=255)
                                                           tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCAGATAC‑A‑ACCCTTCCCCTTTTTCCGGTGCTGCCCCCCGCGGTCttt                                                                                        >  1:1114837/1‑149 (MQ=255)
                                                             cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCGGTCAGGGGAATCCACCCCGACCAGATAC‑A‑GCCCTTCCCCTTTTTCCGTTGCTGCCCCCCGCGTTCTGTcc                                                                                      >  1:6217/1‑151 (MQ=255)
                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑ACCGGCAGGGTAAACCACCCCGACCAGATAA‑A‑ACCCTTCCCCTTTTGCCGTTGCTGCCCCCCGCGGTCTGTCCttttc                                                                                 >  2:291553/1‑149 (MQ=255)
                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGAACATA‑ACCGTCAGGGGAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGGTGCTGCCCCCAGCGTTCTGTCCTTTgt                                                                                 >  2:490682/1‑150 (MQ=255)
                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑GTTTGGGGGATCATA‑ACCGGCAGGGGAAACCACCGCGGCCCGATAC‑A‑CGCCTTCCGCCTTTTCCGTTGCTGCCCCCCGCGTTCTTTCCttttt                                                                                 >  2:564578/1‑149 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTg                                                                                  <  2:1263875/149‑1 (MQ=255)
                                                                    aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCCGGCCGGGGAAACCACCCCGACCAGATAC‑A‑GCCCCTCCGCCTTTTCCGTTGCTGCCGCCCGCGGTCTTTCCTTTtccc                                                                               >  2:1077854/1‑150 (MQ=255)
                                                                    aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTTCCGTTGCTGCCCCCCGCGGTTTGTCCTTTtccg                                                                               >  1:924980/1‑151 (MQ=255)
                                                                        gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAg                                                                           <  1:837010/151‑1 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑ACAGTCAGGGGAATCCACCGCCACCCGAAAC‑A‑GCCCCTCCACTTTTTGCGTTGCTCGCGCCACCGCTTTTTCCTTTtccccccacc                                                                         >  1:760829/1‑147 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGGAGATCATA‑CCGGTCCGGGGAATCCACCGCGACCAGATAC‑A‑ACCCTTCCGCTTTTTCCCGTGCTGGCCCCCGCGTTCTTTCCTTTTccgcccgc                                                                          >  1:1359808/1‑150 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCCGATAC‑A‑GCCCTTCCCCTTTTGTCGTTGCTGCCGCCCGCGTTTTGTCCTTTTCCGCCAGCAACTCTGc                                                                  >  2:808354/1‑151 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑ATTTTGGAGATTATA‑ACGGGCCGGGGAATCCACCCCGACCAGATAC‑A‑GCCCCTCCGCTTTTTCCGTTTCTGCCCCCAGCGGTCTTTCCTTTTCCGCCAGCAActcttc                                                                  >  2:914346/1‑149 (MQ=255)
                                                                                   tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATA‑CCGGGCAGGGGAAACCACCGCGGCCAGATAC‑A‑GCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCtt                                                                >  1:1293150/1‑150 (MQ=255)
                                                                                           ttAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGctc                                                        >  1:504075/1‑151 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑ACGGGCAGGGGAATCCACCCCGACCAGATAC‑A‑GCCCCTCCCCCTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCCACTCTGCTATCCAGctcc                                                       >  2:1025552/1‑150 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGAAGATCATA‑CCCGTCAGGGGAATCCACCCCGACCCGATAC‑C‑GCCCCTCCGCTTTTTCCGTTGCTGCCCCCCGCGGTCTTTCCTTTGCCGCCCGCCACTCTTCTTTCCCGctct                                                       >  2:1266805/1‑151 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑GTTTTGAAGATGATA‑CCGGGCAGGGGAAACCACCGCGACCAGATAC‑AGGCCCTTCC‑CTTTTTCCGGTTGTTCCGCCAGCGTGCTGTCCTTTTTCGCCAGCAACACTGCTATCCAGctct                                                       >  2:572767/1‑151 (MQ=255)
                                                                                             aGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATTATACCCGGCCAGGG‑AATCCACCCCCACCAGATAC‑A‑GCCCCTCCCCTTTTTTCGGTGCTGGCCCCCGCGGTTTTTCCTTTTCCGCCCGCAACCCTTCTAttccaccctc                                                      >  2:1052635/1‑142 (MQ=255)
                                                                                             aGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctct                                                       <  1:1069914/150‑1 (MQ=255)
                                                                                              gTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTTGGAGATCATA‑CCGGGCAGGGGAATCCACCGCGACCAGATAC‑A‑GCCCCTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCCGCAACTCTGCTATCCAGCTCTcc                                                     >  1:1014447/1‑151 (MQ=255)
                                                                                              gTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                                     <  2:464382/151‑1 (MQ=255)
                                                                                                ccGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCgg                                                   <  2:484169/151‑1 (MQ=255)
                                                                                                       aaaTAGCGGGCCCATTGGCGG‑TTTTGGAAGATTATA‑CCGGGCAGGGGAATCCACCCCGACCAGAAAA‑A‑ACCCCTCCGCCTTTGGCGTTTCTGCCGCCCGCGTTCTGGCCTTTTTCGCCAGCAACCCTGCTATCCAGCTCTCCGGCTgggtt                                            >  1:1305407/1‑148 (MQ=255)
                                                                                                               ggCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCAcc                                    <  2:213282/151‑1 (MQ=255)
                                                                                                                gCCCATTGGCGG‑TTTTTGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTTCTGCCGCCAGCGGTCTGTCCTTTGCCGCCCGCAACCCTGCTATCCAGCTCTCCCGCTGGTTGTGGGCAccc                                   >  1:132188/1‑150 (MQ=255)
                                                                                                                  ccATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc                                 <  1:43918/151‑1 (MQ=255)
                                                                                                                  ccATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc                                 <  2:6217/151‑1 (MQ=255)
                                                                                                                    aTTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACtt                               <  2:1337344/151‑1 (MQ=255)
                                                                                                                       ggCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                             >  1:1294298/1‑150 (MQ=255)
                                                                                                                           g‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                        <  1:976351/151‑1 (MQ=255)
                                                                                                                           g‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                        <  1:1365316/151‑1 (MQ=255)
                                                                                                                                 ggAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                    <  2:1005506/150‑1 (MQ=255)
                                                                                                                                 ggAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGc                   >  1:1099102/1‑151 (MQ=255)
                                                                                                                                      aTCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              <  1:1172836/151‑1 (MQ=255)
                                                                                                                                       tCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              >  1:749090/1‑150 (MQ=255)
                                                                                                                                        cATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACaa             <  1:1343761/150‑1 (MQ=255)
                                                                                                                                        cATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat            <  1:808423/151‑1 (MQ=255)
                                                                                                                                        cATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat            >  1:639213/1‑151 (MQ=255)
                                                                                                                                               ggTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt      >  1:903570/1‑151 (MQ=255)
                                                                                                                                                   aGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc  <  1:269339/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGG‑TTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATAC‑A‑GCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  NC_000913/2434733‑2435028

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