Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,8820TG43.2% 150.5 / 35.7 96Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (12/42);  new base G (41/0);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 3.88e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434733‑2435031
                                                                                                                                                     |                                                                                                                                                      
cAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGt                                                                                                                                                        <  2:407988/150‑1 (MQ=255)
     ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACATACCGGTCAGGGGAATcc                                                                                                                                                   >  2:85775/1‑150 (MQ=255)
     ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCa                                                                                                                                                  <  2:757332/151‑1 (MQ=255)
      cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGTAATCCAc                                                                                                                                                 >  1:655055/1‑151 (MQ=255)
      cATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCAc                                                                                                                                                 <  1:488932/151‑1 (MQ=255)
         atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAAACCACcgc                                                                                                                                              >  2:317559/1‑151 (MQ=255)
           ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGAGAGCATACAGGCCCGGGGAAACcccc                                                                                                                                                >  2:464885/1‑147 (MQ=255)
            tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGAGAAAACCGGGATGGGCAATACAcc                                                                                                                                                >  1:7240/1‑146 (MQ=255)
            tGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGa                                                                                                                                            <  2:279601/150‑1 (MQ=255)
             gCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAAACCACCCCGAcc                                                                                                                                          >  1:563354/1‑151 (MQ=255)
              ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATACACcgc                                                                                                                                              >  2:580255/1‑146 (MQ=255)
               cTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGGTGATTCTGGGCAGGGTAATCCCCCGCGAc                                                                                                                                           >  2:320865/1‑148 (MQ=255)
                tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGGCAGGGGAATCCACCGCGAcc                                                                                                                                          >  1:753853/1‑148 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGa                                                                                                                                       <  1:415409/150‑1 (MQ=255)
                 aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                      <  1:470912/151‑1 (MQ=255)
                   gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATAc                                                                                                                                   <  1:221943/151‑1 (MQ=255)
                    ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACa                                                                                                                                  <  2:299076/151‑1 (MQ=255)
                      ttGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTCTACCGGGGAGGGGAATCCACCCCGACCAGATACa                                                                                                                                  >  2:242199/1‑149 (MQ=255)
                        gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGccc                                                                                                                              >  2:268032/1‑151 (MQ=255)
                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATTCACCGCGACCAGATAAAAccc                                                                                                                              >  2:81814/1‑150 (MQ=255)
                         tACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGGTCATACCGGGCCGGGGTAACCAACGCGACCAGATACCGccc                                                                                                                              >  1:493362/1‑150 (MQ=255)
                              gTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCCGGGGAAACCACCCCCACCCGATAAAGCCCTTcccc                                                                                                                        >  2:684223/1‑149 (MQ=255)
                                cATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCCCGACCCGATACAGCCCCTCCCCtt                                                                                                                      >  2:577656/1‑151 (MQ=255)
                                 aTCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGAACAAACCCGGCAGGGGAAACCACCGCCACCAGATACCACCCCTCCCCttt                                                                                                                     >  2:417620/1‑151 (MQ=255)
                                   cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCCGGGGAATCCACCCCGACCCGATAACACCCCTCCCCCtttt                                                                                                                    >  2:537566/1‑151 (MQ=255)
                                   cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCCCCCCGACCCGATACAGCCCCTCCCCtttt                                                                                                                    >  2:45265/1‑150 (MQ=255)
                                     ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                                                   <  1:566541/149‑1 (MQ=255)
                                       cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCCTCCCCTTTTTcc                                                                                                                 >  1:185108/1‑149 (MQ=255)
                                       cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCg                                                                                                                <  2:292115/150‑1 (MQ=255)
                                       cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGAAGATCATACCGGGCAAGGGGAACCACCGCGACCAGATACAGCCCCTCCGCTTTTGcc                                                                                                                 >  2:557074/1‑149 (MQ=255)
                                          gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTTTCGCTTTTTCCGtt                                                                                                              >  2:397929/1‑149 (MQ=255)
                                           cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATTCACCCCGGCCAGATAAAGCCCCTCCGCTTTTTCCGtt                                                                                                              >  1:387691/1‑148 (MQ=255)
                                            tGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATACCCCGCGACCACATAAAGCCCTTCCGCTTTTTCCGTtgct                                                                                                           >  2:472342/1‑150 (MQ=255)
                                               aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCCGGGGAATCCACCCCGACCCGATACAACCCCTCCCCTTTTTCCCGTGCTtcc                                                                                                        >  2:376069/1‑150 (MQ=255)
                                               aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgccg                                                                                                       <  2:143309/151‑1 (MQ=255)
                                               aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTgcc                                                                                                        <  1:252001/150‑1 (MQ=255)
                                                   aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCa                                                                                                    <  1:186186/150‑1 (MQ=255)
                                                       tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAACCCCTCCCCTTTTTCCCTTGCTGCCCCCAGCGtt                                                                                               >  1:391427/1‑151 (MQ=255)
                                                       tGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGt                                                                                                <  1:40962/150‑1 (MQ=255)
                                                              aGCCTTAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCt                                                                                        <  1:483952/151‑1 (MQ=255)
                                                              aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCt                                                                                        <  2:407654/151‑1 (MQ=255)
                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                      <  2:315290/151‑1 (MQ=255)
                                                                ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCttt                                                                                      <  2:557025/151‑1 (MQ=255)
                                                                 cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGCCCCTCCCCTTTTTCCCTTGCTTCCGCCCGCGTtcttttcttt                                                                                      >  2:424563/1‑150 (MQ=255)
                                                                  tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGCTGCTGCCCCCAGCGTTCTGTCCttt                                                                                      >  2:718050/1‑149 (MQ=255)
                                                                   gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAACCCCTCCCCTTTTTCCGTTGCTGCCGCCCGCGGTCTGTCCTTTtcc                                                                                   >  2:649976/1‑151 (MQ=255)
                                                                    aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAAACCACCGCGACCAGAAACAACCCTTCCGCTTTTTCCGTTGCTGCCCCCAGCGTTCTGTCCTTTtcc                                                                                   >  1:112412/1‑150 (MQ=255)
                                                                    aaaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCCCGACCAGATACAACCCCTCCCCCTTTTCCGTTGCTGCCGCCCGCGGTCTGTTCttt                                                                                      >  2:496435/1‑147 (MQ=255)
                                                                     aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCCCCCCGACCAGATACAGCCCTTCCGCCTTTTCCGTTGCTGCCGCCCGCGTTCTGTCCTTTgcccc                                                                                 >  1:782572/1‑149 (MQ=255)
                                                                       tGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGTCAGGGGAAATCACCGCGGCCAGATACAGCCCTTCCGCTTTTTCCGTTTCTGCCCCCCGCGGTCTGTCCTTTtccgcc                                                                                >  2:84123/1‑150 (MQ=255)
                                                                         ttCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCATAACGGTCAGGGGAATCCACCCCCACCCGATAAAACCCTTCCGCTTTTTCCGTTGGTGGCCCCCGCGGTCTGTCCTTTtccgcc                                                                                >  2:479726/1‑148 (MQ=255)
                                                                          tCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGTCAGGGGAAACCACCCCGACCCGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCCGCGGTCTTTCCTTTtccgcc                                                                                >  2:289254/1‑147 (MQ=255)
                                                                           ccGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCa                                                                            <  2:702230/150‑1 (MQ=255)
                                                                             gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAAct                                                                         <  1:65900/151‑1 (MQ=255)
                                                                                 attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATACCGGGCAAGGTAAACCACCCCGGCCCGAAACAGCCCTTCCCCTTTTTCCCTTGCCGCCCCCAGCGTTCTTTCCtttttcccccgcaa                                                                           >  2:765112/1‑134 (MQ=255)
                                                                                   tattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCt                                                                    <  1:804719/150‑1 (MQ=255)
                                                                                       ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATTCTACCGGGCAGGGGAAACCACCGCGACCAGAGAAAGAACTTTCGCTTTTTCCGGTGGGGCCGGCAGGGGTCTTTTCTTTTCCGCCAGCaaaccctccctccc                                                               >  2:469093/1‑139 (MQ=255)
                                                                                         cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCa                                                               <  1:795286/149‑1 (MQ=255)
                                                                                         cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCa                                                               <  2:603321/149‑1 (MQ=255)
                                                                                          gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGGCAGGGGAAAACACCCCGACCCGAAACAGCCCTTTCCCTTTTTCCGGTTCTGGCGCCAGCGCTCTGTCCCTTTCCGCCAGCAAACCTGCTTTCCCGc                                                             >  1:381353/1‑150 (MQ=255)
                                                                                            tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGGAGATTAAACCGCCCCGGGTAATCCACCCCGACCCAATACAGCCCTTCCGCTTTTTTCCGTGGTGCCCCCCCAGTTCTTTCCTTTGCCCCCCACCACTCTGCTATCCaaacc                                                           >  2:474803/1‑146 (MQ=255)
                                                                                             aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCCGCAACTCTGCTATCCAGCtctc                                                         >  1:128503/1‑151 (MQ=255)
                                                                                              gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCCGCGTTCTTTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTcc                                                        >  1:291719/1‑151 (MQ=255)
                                                                                              gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                                        <  1:31126/151‑1 (MQ=255)
                                                                                                 cGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGc                                                     <  2:296962/151‑1 (MQ=255)
                                                                                                  gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCt                                                    <  2:773065/151‑1 (MQ=255)
                                                                                                  gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCt                                                    <  1:150150/151‑1 (MQ=255)
                                                                                                  gCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGTGAGGGGAATCCACCGCGGCCAGATACAGGCCTTTCGCTTTTGCCGGTGCTGCCGCCCGCGGTCTGTCCTTTGCCGCCCGCCACTCTGCTTTCCCGCTCCCCgg                                                      >  2:202018/1‑149 (MQ=255)
                                                                                                      gAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATACACCGCGACCCGATACAGCCCCTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTCCGGCTgg                                                  >  2:688139/1‑149 (MQ=255)
                                                                                                         aTAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCCCCCGCGCTCTGTCCTTTTCCGCCCGCAAATCTTCCATCCCGCTCTCCGGCTGGGtg                                               >  2:172987/1‑149 (MQ=255)
                                                                                                         aTAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgtg                                             <  2:523005/151‑1 (MQ=255)
                                                                                                          tAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCCGATACAGCCCCTCCCCTTTTGCCGGTGCTGCCGCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCCCCGGCTGGTTGTgg                                            >  1:33554/1‑151 (MQ=255)
                                                                                                          tAGCGGGCCCATTGGCGGTTTTGGAAGAACATAACGGTCAGGGGAATCAACCCCGGCCAGATAAAGCCCTTCCGCTTTTTTCGCTGCTTCCGCCAGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATTCCGCCCTCCCGCCgg                                                  >  2:650312/1‑145 (MQ=255)
                                                                                                             cGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCa                                         <  1:686786/151‑1 (MQ=255)
                                                                                                                   cATTGGCGGTTTTTGAAGGTCATACCGGGCAGGGGAATCCACCCCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCCCCCGCGTTCTGTCCTTTTCCGCCCGCAACTCTGCTATCCAGCTCTCCCGCTGGGTGTGGGGAcc                                       >  2:245540/1‑147 (MQ=255)
                                                                                                                        gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                              >  1:326463/1‑151 (MQ=255)
                                                                                                                        gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                              <  1:389256/151‑1 (MQ=255)
                                                                                                                        gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                              >  2:461675/1‑151 (MQ=255)
                                                                                                                           gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                           <  1:730044/151‑1 (MQ=255)
                                                                                                                              ttGGAAGATTATACCGGTCAGGGTAATCCACCGCGACCAGACACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGCTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGTGCACCGACTTCCATCAg                           >  2:374153/1‑148 (MQ=255)
                                                                                                                              ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                        >  1:105031/1‑151 (MQ=255)
                                                                                                                               tGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCt                        <  1:327579/150‑1 (MQ=255)
                                                                                                                                  aaGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                    <  1:563847/151‑1 (MQ=255)
                                                                                                                                   aGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                   <  2:427798/151‑1 (MQ=255)
                                                                                                                                     aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                 <  2:585514/151‑1 (MQ=255)
                                                                                                                                     aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                 <  1:471510/151‑1 (MQ=255)
                                                                                                                                      tCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                 >  2:781018/1‑150 (MQ=255)
                                                                                                                                       cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat               <  2:692770/151‑1 (MQ=255)
                                                                                                                                          aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATAtt            <  2:714643/151‑1 (MQ=255)
                                                                                                                                            cGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc          <  1:155893/151‑1 (MQ=255)
                                                                                                                                              gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg        >  1:575476/1‑151 (MQ=255)
                                                                                                                                              gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGACTGCCGACAATATTCCTg        <  2:234771/151‑1 (MQ=255)
                                                                                                                                                cAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAc      >  1:85366/1‑151 (MQ=255)
                                                                                                                                                 aGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc     >  1:492364/1‑151 (MQ=255)
                                                                                                                                                   ggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat   <  2:67183/151‑1 (MQ=255)
                                                                                                                                                    gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  2:671537/151‑1 (MQ=255)
                                                                                                                                                    gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  1:409869/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434733‑2435031

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: