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

CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  NC_000913/2434733‑2435028
                                                                                                                                                       |                                                                                                                                                    
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGt                                                                                                                                                      <  2:951734/150‑1 (MQ=255)
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTa                                                                                                                                                     <  2:86925/151‑1 (MQ=255)
  ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGTCCGGGGaa                                                                                                                                                    >  2:1278192/1‑150 (MQ=255)
   gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATc                                                                                                                                                  <  1:955876/151‑1 (MQ=255)
    tCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGTCCGGGGAAAcc                                                                                                                                                 >  1:365424/1‑151 (MQ=255)
       atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCAcc                                                                                                                                              >  1:213292/1‑151 (MQ=255)
       atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAcc                                                                                                                                              <  2:439372/151‑1 (MQ=255)
       atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAcc                                                                                                                                              <  2:1079874/151‑1 (MQ=255)
             gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  1:888804/1‑151 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATCCACCGCGAcc                                                                                                                                        >  1:304629/1‑147 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGa                                                                                                                                     <  2:371826/150‑1 (MQ=255)
                    ccTTGTACTCGTCATCTTTCAGGTTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACa                                                                                                                                <  1:1328923/151‑1 (MQ=255)
                       tgaaCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGTAGATCATACCGGGCAGGGGAATCCACCCGGACCAGAAACAAcc                                                                                                                            >  1:64856/4‑151 (MQ=255)
                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGGAAATCATAACCGTCAAGGGAATCCACCCCCACCAGAAACAccccc                                                                                                                          >  1:849193/1‑150 (MQ=255)
                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGACCAGAAACAGcccc                                                                                                                          >  1:510791/1‑150 (MQ=255)
                            tCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCCGGGGAAACCACCCCGACCAGATACAACCCCTcc                                                                                                                       >  2:969420/1‑151 (MQ=255)
                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTcccc                                                                                                                     >  1:411881/1‑149 (MQ=255)
                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGc                                                                                                                     <  1:158475/151‑1 (MQ=255)
                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAAACCACCGCGACCCGATACAACCCCTCCGc                                                                                                                     >  1:303544/1‑151 (MQ=255)
                              gTCATCTTTCAGGCTGTAACGACGTTGGCTACAGAATGACATGTTCCGAGCATTATTCTCGTAAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTAGAAGATCATACCGGTCAGGGTAAAACACCGCGACCAGAGACAGCCCTTCCGc                                                                                                                     <  2:755100/151‑1 (MQ=255)
                                cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCt                                                                                                                    <  2:166069/150‑1 (MQ=255)
                                 aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATAACGGGCAGGGGAATCCACCGCGGCCCAATACAACCCCTCCGCttt                                                                                                                  >  1:502733/1‑151 (MQ=255)
                                     ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCttt                                                                                                                  <  1:492742/147‑1 (MQ=255)
                                      tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGGCCCGATACAACCCCTTCCCCTTTTccc                                                                                                              >  2:883628/1‑151 (MQ=255)
                                        aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                                           <  2:102281/151‑1 (MQ=255)
                                           cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATAATACCGGGCAGGGGAATACACCGCGACCCGATAAAACCCTTCCCCTTTTTCCGTtgct                                                                                                        >  2:1338391/1‑151 (MQ=255)
                                            tGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATTCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGCtgccc                                                                                                       >  1:1083127/1‑149 (MQ=255)
                                             gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgctg                                                                                                       <  1:1334507/150‑1 (MQ=255)
                                             gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGAAACAGCCCTTCCGCTTTTGCCGTtgctgc                                                                                                      >  1:1131285/1‑151 (MQ=255)
                                              tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCCCGACCAGAAAAAGCCCCTCCCCTTTTTCCGGTGCTgcc                                                                                                     >  1:305332/1‑151 (MQ=255)
                                                     gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCCGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCcgcg                                                                                              >  1:974284/1‑151 (MQ=255)
                                                        ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGtt                                                                                            <  2:693894/150‑1 (MQ=255)
                                                           tACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCCCTTTTTCCGGTGCTGCCCCCCGCGGTCttt                                                                                        >  1:1114837/1‑149 (MQ=255)
                                                             cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCCCCCGCGTTCTGTcc                                                                                      >  1:6217/1‑151 (MQ=255)
                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGGCAGGGTAAACCACCCCGACCAGATAAAACCCTTCCCCTTTTGCCGTTGCTGCCCCCCGCGGTCTGTCCttttc                                                                                 >  2:291553/1‑149 (MQ=255)
                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATAACCGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGGTGCTGCCCCCAGCGTTCTGTCCTTTgt                                                                                 >  2:490682/1‑150 (MQ=255)
                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGGGGATCATAACCGGCAGGGGAAACCACCGCGGCCCGATACACGCCTTCCGCCTTTTCCGTTGCTGCCCCCCGCGTTCTTTCCttttt                                                                                 >  2:564578/1‑149 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTg                                                                                  <  2:1263875/149‑1 (MQ=255)
                                                                    aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCCGGGGAAACCACCCCGACCAGATACAGCCCCTCCGCCTTTTCCGTTGCTGCCGCCCGCGGTCTTTCCTTTtccc                                                                               >  2:1077854/1‑150 (MQ=255)
                                                                    aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCCCCCGCGGTTTGTCCTTTtccg                                                                               >  1:924980/1‑151 (MQ=255)
                                                                        gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAg                                                                           <  1:837010/151‑1 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACAGTCAGGGGAATCCACCGCCACCCGAAACAGCCCCTCCACTTTTTGCGTTGCTCGCGCCACCGCTTTTTCCTTTtccccccacc                                                                         >  1:760829/1‑147 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCCGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCCGTGCTGGCCCCCGCGTTCTTTCCTTTTccgcccgc                                                                          >  1:1359808/1‑150 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGACCCGATACAGCCCTTCCCCTTTTGTCGTTGCTGCCGCCCGCGTTTTGTCCTTTTCCGCCAGCAACTCTGc                                                                  >  2:808354/1‑151 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGGAGATTATAACGGGCCGGGGAATCCACCCCGACCAGATACAGCCCCTCCGCTTTTTCCGTTTCTGCCCCCAGCGGTCTTTCCTTTTCCGCCAGCAActcttc                                                                  >  2:914346/1‑149 (MQ=255)
                                                                                   tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAAACCACCGCGGCCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCtt                                                                >  1:1293150/1‑150 (MQ=255)
                                                                                           ttAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCAGCAACTCTGCTATCCAGctc                                                        >  1:504075/1‑151 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCCCCTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCCACTCTGCTATCCAGctcc                                                       >  2:1025552/1‑150 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAATCCACCCCGACCCGATACCGCCCCTCCGCTTTTTCCGTTGCTGCCCCCCGCGGTCTTTCCTTTGCCGCCCGCCACTCTTCTTTCCCGctct                                                       >  2:1266805/1‑151 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGAAGATGATACCGGGCAGGGGAAACCACCGCGACCAGATACAGGCCCTTCCCTTTTTCCGGTTGTTCCGCCAGCGTGCTGTCCTTTTTCGCCAGCAACACTGCTATCCAGctct                                                       >  2:572767/1‑151 (MQ=255)
                                                                                             aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATACCCGGCCAGGGAATCCACCCCCACCAGATACAGCCCCTCCCCTTTTTTCGGTGCTGGCCCCCGCGGTTTTTCCTTTTCCGCCCGCAACCCTTCTAttccaccctc                                                      >  2:1052635/1‑142 (MQ=255)
                                                                                             aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGctct                                                       <  1:1069914/150‑1 (MQ=255)
                                                                                              gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCGCTTTTTCCGTTGCTGCCGCCAGCGGTCTGTCCTTTGCCGCCCGCAACTCTGCTATCCAGCTCTcc                                                     >  1:1014447/1‑151 (MQ=255)
                                                                                              gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                                     <  2:464382/151‑1 (MQ=255)
                                                                                                ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCgg                                                   <  2:484169/151‑1 (MQ=255)
                                                                                                       aaaTAGCGGGCCCATTGGCGGTTTTGGAAGATTATACCGGGCAGGGGAATCCACCCCGACCAGAAAAAACCCCTCCGCCTTTGGCGTTTCTGCCGCCCGCGTTCTGGCCTTTTTCGCCAGCAACCCTGCTATCCAGCTCTCCGGCTgggtt                                            >  1:1305407/1‑148 (MQ=255)
                                                                                                               ggCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCAcc                                    <  2:213282/151‑1 (MQ=255)
                                                                                                                gCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTTCTGCCGCCAGCGGTCTGTCCTTTGCCGCCCGCAACCCTGCTATCCAGCTCTCCCGCTGGTTGTGGGCAccc                                   >  1:132188/1‑150 (MQ=255)
                                                                                                                  ccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc                                 <  1:43918/151‑1 (MQ=255)
                                                                                                                  ccATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGAc                                 <  2:6217/151‑1 (MQ=255)
                                                                                                                    aTTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACtt                               <  2:1337344/151‑1 (MQ=255)
                                                                                                                       ggCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                             >  1:1294298/1‑150 (MQ=255)
                                                                                                                           gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                        <  1:976351/151‑1 (MQ=255)
                                                                                                                           gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                        <  1:1365316/151‑1 (MQ=255)
                                                                                                                                 ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                    <  2:1005506/150‑1 (MQ=255)
                                                                                                                                 ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGc                   >  1:1099102/1‑151 (MQ=255)
                                                                                                                                      aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              <  1:1172836/151‑1 (MQ=255)
                                                                                                                                       tCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa              >  1:749090/1‑150 (MQ=255)
                                                                                                                                        cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACaa             <  1:1343761/150‑1 (MQ=255)
                                                                                                                                        cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat            <  1:808423/151‑1 (MQ=255)
                                                                                                                                        cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat            >  1:639213/1‑151 (MQ=255)
                                                                                                                                               ggTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt      >  1:903570/1‑151 (MQ=255)
                                                                                                                                                   aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc  <  1:269339/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACC  >  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: