Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG41.2% 125.1 / 29.1 80Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (14/33);  new base G (33/0);  total (47/33)
Fisher's exact test for biased strand distribution p-value = 1.10e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATAT  >  NC_000913/2434738‑2435032
                                                                                                                                                 |                                                                                                                                                          
ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                                                      <  2:731644/151‑1 (MQ=255)
ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                                                      <  1:443320/151‑1 (MQ=255)
  atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAcc                                                                                                                                                   <  2:113028/151‑1 (MQ=255)
   tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCAcc                                                                                                                                                   >  1:29451/1‑150 (MQ=255)
    atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGGTAATACCGGGCCGGGGAAATCAcc                                                                                                                                                   >  2:446179/1‑149 (MQ=255)
    atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACcg                                                                                                                                                  <  1:561358/150‑1 (MQ=255)
      ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTGAATTCTACCGGTCCGGGGAATCACCcgc                                                                                                                                                 >  1:488979/1‑149 (MQ=255)
       tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTAAACCCGGCAGGGGaaac                                                                                                                                                        >  1:762699/1‑140 (MQ=255)
            aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGTAATCCACCCCGAcc                                                                                                                                             >  2:361682/1‑147 (MQ=255)
                   gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGACCCGAAACAGccc                                                                                                                                >  1:409147/1‑151 (MQ=255)
                      cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGGAGATCATAACGGGCAGGGGAAAACACCGCGACCAAAAAAAGCCCGTc                                                                                                                             >  1:678083/1‑151 (MQ=255)
                       tCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCCGGGGAATCCACCCCGACCAGATACAGCCCTTcc                                                                                                                            >  1:352638/1‑151 (MQ=255)
                        cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTc                                                                                                                             <  1:768000/149‑1 (MQ=255)
                        cGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCg                                                                                                                           <  2:96026/151‑1 (MQ=255)
                          tCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATTCACCGCGACCAGATACAGCCCTTTCGc                                                                                                                          >  2:695082/1‑150 (MQ=255)
                           cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtt                                                                                                                        <  2:608700/151‑1 (MQ=255)
                                ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGAACATACCCGGCCGGGGAAATCACCCCGCCCCGATACAGCCCCTCCCCttttttc                                                                                                                   >  2:662719/1‑148 (MQ=255)
                                ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCtttt                                                                                                                      <  1:81906/148‑1 (MQ=255)
                                ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCttt                                                                                                                       <  2:87483/147‑1 (MQ=255)
                                   aGGCTGTAACGACGTGGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCTGAAATAGCGGGCCCATTGGCGGTTTTGGAAGGTCATACCGGTCGGGGTAATCCACCGCGACCAGATACAGAACTTCCGCTTTTGCCGtt                                                                                                                <  1:616399/151‑1 (MQ=255)
                                        gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgctg                                                                                                            <  1:485943/150‑1 (MQ=255)
                                        gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtgct                                                                                                             <  1:453217/149‑1 (MQ=255)
                                        gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCACTTCCGCTTTTGCCGTtgct                                                                                                             <  1:92625/149‑1 (MQ=255)
                                        gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGAAACAGCCCTTCCCCTTTTTcc                                                                                                                   >  2:724278/1‑143 (MQ=255)
                                         tAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGACCCGATACAGCCCCTCCGCCTTTTCCGtt                                                                                                                >  1:737954/1‑145 (MQ=255)
                                          aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCATACCGGTCAGGGGAATACACCGCGACCCGATACAGCCCTTCCGCTTTTTCCGTTTCTgcc                                                                                                          >  1:107949/1‑150 (MQ=255)
                                            cTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCCCCTTTTCCCTTTCTgccgcc                                                                                                       >  2:215824/1‑151 (MQ=255)
                                                gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCCCCCCGACCAGAAACAGCCCTTCCGCTTTTTCCGTTGCTgccgcc                                                                                                       >  1:758115/1‑147 (MQ=255)
                                                  tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAACCGGTCCGGGGAATCCACCGCGACCAGATAAAAGCCCTCCCCCTTTTTCGTTGCTgccgcc                                                                                                       >  1:206384/1‑145 (MQ=255)
                                                  tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGt                                                                                                  <  1:160552/150‑1 (MQ=255)
                                                   ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATACCGGGCAGGGGAAACCCCCGCGACCCGATACCGCCCCTCCGCCTTTGCCGTTGCTTCCGCCcgcggt                                                                                                 >  1:629076/1‑148 (MQ=255)
                                                   ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGGAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCCCCcgcg                                                                                                   >  1:159714/1‑148 (MQ=255)
                                                   ggCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATTATACCCGTCCGGGGAATCCACCGCGACCCGATACAACCCTTCCGCTTTTTCCGTTGCTGCCCCCAGCGt                                                                                                  >  2:229872/1‑149 (MQ=255)
                                                    gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCt                                                                                               <  2:136618/151‑1 (MQ=255)
                                                     cTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGTCAGGGGAATCCACCCCGACCCGATACAACCCCTCCGCTTTTTCCGGTGCTTCCGCCcgcggt                                                                                                 >  2:301606/1‑146 (MQ=255)
                                                          gCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCtt                                                                                        <  2:240887/151‑1 (MQ=255)
                                                           ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                       <  2:110064/151‑1 (MQ=255)
                                                           ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCAATTGGCGGTTTTGGAAGATCATACTGGTCAGGGTAATCCAGCGCGACCAGATACAGCCCTTCCTCAGTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                       <  1:543473/151‑1 (MQ=255)
                                                               aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccg                                                                                   <  1:199532/151‑1 (MQ=255)
                                                                 aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                                                <  2:737258/121‑1 (MQ=255)
                                                                 aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccgcc                                                                                 <  2:82813/151‑1 (MQ=255)
                                                                 aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGATTTTGAAGGTCATTCCCGGCAGGGGAATCCCCCCCGACCCGACACAGCCCCTTCCGTTTTGCCGGTGGTGCCGCCCGCGTTCTGTCCTTTgccgcc                                                                                 >  2:872057/1‑151 (MQ=255)
                                                                  tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATAAAAACCGTCCGGGGGAACCACCGCGCACACATACAGCGCCTTCGCTTTTTGCGTTGCCTGCGCCCGCGTTCTTTTCTTTtcccgc                                                                                 >  1:775350/1‑147 (MQ=255)
                                                                   gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                                <  2:362413/150‑1 (MQ=255)
                                                                   gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                                <  2:505180/150‑1 (MQ=255)
                                                                     tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCCGATACAACCCCTCCGCTTTTTCCGGTGCTGCCGCCCGCGGTCTGTCCTTTGccgcccgc                                                                              >  1:850570/1‑150 (MQ=255)
                                                                      ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCACCCCGGCCAGATACAGCCCCTCCGCTTTTTCCGTTGCTGCCCCCCGCGGTGTGTCCtttttcccc                                                                                 >  1:111810/1‑140 (MQ=255)
                                                                        gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATAATAACGGTCAGGGGAATCCACCGCGACAAGATACAGCCCTTTCGCGTTTGTCGGTGGTGCCGCCAGCGGTCTGTCCTTTTCCGCCAGCaa                                                                            >  2:29960/1‑149 (MQ=255)
                                                                            attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATACACCGCGAACAAATAAAGCCCTTCCCCTTTTTTCGTTGCTGCCCCCAGCGGTCTGTCCTTTTCCGCCCGCAActct                                                                        >  1:691539/1‑149 (MQ=255)
                                                                             ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATTCACCGCGAACAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCCGCAActc                                                                         >  1:601176/1‑147 (MQ=255)
                                                                                ttCTCGTGAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTCTAGATCATCCCCGTGCGGGGGATTACACCCGGCACGAGAAAACCCCTCCGCCTTTTTCGTTGCTCTCGCCACCGCGCTGTCCTTTGTCGCCAGCAACTTTTCTATc                                                                  >  1:492410/1‑151 (MQ=255)
                                                                                 tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCCCCAGCGTTCTGTCCTTTTCCCCCAGCAACTCTGCttt                                                                   >  2:541705/1‑147 (MQ=255)
                                                                                 tctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTGAAGATCATAACGGGCAGGGTAATTCACCCCGCCCAGATACAGCCCCTTCCCCTTTGCCGTTTCTCCCGCCCGCGTTCTTTCCCTTTCCGCCAGCaa                                                                            >  1:681980/1‑140 (MQ=255)
                                                                                  ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAAACCACCCCGACCAGATACAGCCCTTTCCCCTTTGCCGTTGCTGCCCCCCGCGTTCTGTCCTTTTCCCCCCGCAACTCTGCTAt                                                                   >  1:690888/1‑148 (MQ=255)
                                                                                    cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGTAGGGGAATACACCCCCGACCAATACAGCCCTTCCCCTTTCCCCGGTGGTGCCCGCCGCGTTCTGATCTTTGCCGGCCACAACTCTGCTATc                                                                  >  1:793589/1‑147 (MQ=255)
                                                                                         gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATAAAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCCGCAACTCTGCTCTCCAACGCTcc                                                         >  1:395142/1‑151 (MQ=255)
                                                                                           ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCgg                                                       <  1:441124/151‑1 (MQ=255)
                                                                                           ccGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCg                                                        <  2:468400/150‑1 (MQ=255)
                                                                                           ccGCCAGAAATAGCGGGCCCATTGGCGGATTTTGAAAAACATACCCGGCAGGGTAATACACCGCGACCAAACAAACCCCCTCCCGTTTTGGCGTTGTTGCCCCCAACGGTCTGTGCCTTTTCCCCCGCAAATCTCCTATccccgtcctc                                                         >  2:755716/1‑141 (MQ=255)
                                                                                            cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                                      <  1:649580/151‑1 (MQ=255)
                                                                                                  aaaTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCCGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCCGCTgg                                                   >  2:499742/1‑148 (MQ=255)
                                                                                                    aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgtg                                              <  1:29773/151‑1 (MQ=255)
                                                                                                     tAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCACCAATCTGCTATCCAGCTCTCCCGCTGGGTGTgg                                             >  1:480328/1‑151 (MQ=255)
                                                                                                      aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTggg                                            <  1:720264/151‑1 (MQ=255)
                                                                                                          ggCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGTAATCCACCGCGACCCGATACAGCCCTTTCGCTTTTTCCCTTGCTGCCCCCCGCGTTCTGTCCTTTTCCCCCCGCCACACTGCTATCCAGCTCTCCGGCTGGGTGTGGGCAcc                                        >  1:530659/1‑151 (MQ=255)
                                                                                                                 tGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTcc                                 >  1:722736/1‑151 (MQ=255)
                                                                                                                   gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                               >  2:40258/1‑151 (MQ=255)
                                                                                                                     ggTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCa                             <  1:831089/151‑1 (MQ=255)
                                                                                                                       ttttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                           >  2:378391/1‑151 (MQ=255)
                                                                                                                       ttttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                           >  1:678171/1‑151 (MQ=255)
                                                                                                                           ggAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                        <  2:21877/150‑1 (MQ=255)
                                                                                                                             aaGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                     >  1:511105/1‑151 (MQ=255)
                                                                                                                               gATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGAc                   >  2:805656/1‑151 (MQ=255)
                                                                                                                                   aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat                >  2:199385/1‑150 (MQ=255)
                                                                                                                                   aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat                >  1:720100/1‑150 (MQ=255)
                                                                                                                                   aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGATTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAata               <  2:36886/151‑1 (MQ=255)
                                                                                                                                         gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg         >  2:14074/1‑151 (MQ=255)
                                                                                                                                            cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc       <  2:485972/151‑1 (MQ=255)
                                                                                                                                            cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc       <  2:105990/151‑1 (MQ=255)
                                                                                                                                                gTAATGCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat    >  1:86999/1‑150 (MQ=255)
                                                                                                                                                 tAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCatat  >  2:205329/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                          
CCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATAT  >  NC_000913/2434738‑2435032

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