Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,140,834 A→T 6.7% T147S (ACC→TCC)  rluC → 23S rRNA pseudouridine(955,2504,2580) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,140,8340AT 6.7% 635.9 / 28.3 210T147S (ACC→TCC) rluC23S rRNA pseudouridine(955,2504,2580) synthase
Reads supporting (aligned to +/- strand):  ref base A (100/96);  new base T (5/9);  total (105/105)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 216 total aligned reads displayed.

GAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTCAAGCACCGTTA  >  CP009273/1140687‑1140980
                                                                                                                                                   |                                                                                                                                                  
gAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGaca                                                                                                                                                    <  2:1682042/148‑1 (MQ=255)
    atgatCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTc                                                                                                                                                <  2:1904387/148‑1 (MQ=255)
    atgatCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTc                                                                                                                                                <  2:110164/148‑1 (MQ=255)
    atgatCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACAGCGCGGCAGTGGATAAAGCTTAGGCGGTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTc                                                                                                                                                >  1:3489820/1‑148 (MQ=255)
     tgatCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCa                                                                                                                                               >  1:3490839/1‑148 (MQ=255)
     tgatCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCa                                                                                                                                               >  2:1367143/1‑148 (MQ=255)
       atCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACTCCTCAgg                                                                                                                                             <  1:912926/148‑1 (MQ=255)
       atCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAgg                                                                                                                                             <  2:144889/148‑1 (MQ=255)
           caTCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGtt                                                                                                                                         <  1:2331672/148‑1 (MQ=255)
            aTCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACTCCTCAGGTGttt                                                                                                                                        <  1:2709391/148‑1 (MQ=255)
            aTCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGttt                                                                                                                                        <  2:1376069/148‑1 (MQ=255)
             tCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGtttt                                                                                                                                       >  1:786903/1‑148 (MQ=255)
              ccTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTg                                                                                                                                      >  1:2092032/1‑148 (MQ=255)
                tGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGGTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCt                                                                                                                                    >  1:1704225/1‑148 (MQ=255)
                            aaCCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGaaaaa                                                                                                                        >  2:1/1‑148 (MQ=255)
                            aaCCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGaaaaa                                                                                                                        >  2:24/1‑148 (MQ=255)
                             aCCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGaaaaaa                                                                                                                       <  1:219137/148‑1 (MQ=255)
                                 tCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCt                                                                                                                   <  1:425534/148‑1 (MQ=255)
                                 tCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCt                                                                                                                   >  1:81603/1‑148 (MQ=255)
                                      tACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGgcg                                                                                                              <  1:2902687/148‑1 (MQ=255)
                                       aCGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGgcgt                                                                                                             <  2:615449/148‑1 (MQ=255)
                                         ggcggTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGgcgttg                                                                                                           <  2:2728347/148‑1 (MQ=255)
                                          gcggTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGgcgttgc                                                                                                          >  2:3343285/1‑148 (MQ=255)
                                             gTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGgcgttgcgtt                                                                                                       <  1:1999260/148‑1 (MQ=255)
                                              tACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTtc                                                                                                      <  1:2000497/148‑1 (MQ=255)
                                              tACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTtc                                                                                                      >  2:2487777/1‑148 (MQ=255)
                                              tACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTtc                                                                                                      >  2:352345/1‑148 (MQ=255)
                                               aCATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTtct                                                                                                     <  2:2080165/148‑1 (MQ=255)
                                                  tGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTg                                                                                                  <  2:2273430/148‑1 (MQ=255)
                                                  tGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAAAGCTCGGCGTTGCGTTCTCTg                                                                                                  >  1:675642/1‑148 (MQ=255)
                                                       gcAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGa                                                                                             >  1:192252/1‑148 (MQ=255)
                                                          gTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACTCCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCa                                                                                          <  1:1578604/148‑1 (MQ=255)
                                                           tGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCaa                                                                                         <  2:1608724/148‑1 (MQ=255)
                                                           tGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCaa                                                                                         >  1:851303/1‑148 (MQ=255)
                                                           tGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCaa                                                                                         <  2:1735630/148‑1 (MQ=255)
                                                            ggTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCCGAAGCGCGGTTACTTGAAATGGTTAATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAAAGCTCGGCCTTGCGTTCACAGCATGATCAAt                                                                                        >  1:1828604/1‑148 (MQ=255)
                                                              tttAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTa                                                                                      <  1:559525/148‑1 (MQ=255)
                                                                tAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACg                                                                                    <  1:3059512/148‑1 (MQ=255)
                                                                 aaGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGt                                                                                   >  1:2377802/1‑148 (MQ=255)
                                                                 aaGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGt                                                                                   >  1:593358/1‑148 (MQ=255)
                                                                  aGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCACTTACGTg                                                                                  >  1:3379543/1‑148 (MQ=255)
                                                                  aGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCACCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTg                                                                                  <  2:2612149/148‑1 (MQ=255)
                                                                   gCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACTCCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGa                                                                                 <  1:2979456/148‑1 (MQ=255)
                                                                   gCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGa                                                                                 >  2:1729314/1‑148 (MQ=255)
                                                                   gCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAAAGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGa                                                                                 >  1:1598841/1‑148 (MQ=255)
                                                                    cTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGTCCGGGACTCCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGaa                                                                                >  1:967091/1‑148 (MQ=255)
                                                                    cTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGGTCCTTGAACAGGTTCAACGTCTTGACCGGGACTCCTCAGGTGTTTTGCTGGTAGCGAAAAAAAGCTCGGCGTTGCGTTCTCTGCATGAGCAATTAAGTGaa                                                                                >  1:1207680/1‑148 (MQ=255)
                                                                      tCGGCGTTATTGAAGGTTTGCGGGCGGTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGaaaa                                                                              <  2:2995233/148‑1 (MQ=255)
                                                                       cGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTAGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGaaaaa                                                                             >  2:290339/1‑148 (MQ=255)
                                                                        ggCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAg                                                                            <  2:1402212/148‑1 (MQ=255)
                                                                        ggCGTTATTGAAGGTGTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTATGCATGAGCAATTACGTGAAAAAg                                                                            <  2:2889473/148‑1 (MQ=255)
                                                                            ttATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGAt                                                                        >  1:1992529/1‑148 (MQ=255)
                                                                            ttATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTCCTCTGCATGAGCAATTACGTGAAAAAGGGAt                                                                        <  1:2574619/148‑1 (MQ=255)
                                                                                tGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCaa                                                                    >  2:1873492/1‑148 (MQ=255)
                                                                                  aaGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCaaaa                                                                  >  1:2400988/1‑148 (MQ=255)
                                                                                   aGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACTCCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCaaaaa                                                                 <  2:233367/148‑1 (MQ=255)
                                                                                   aGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCaaaaa                                                                 <  2:43693/148‑1 (MQ=255)
                                                                                    ggTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAg                                                                <  2:2426391/148‑1 (MQ=255)
                                                                                    ggTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAg                                                                >  1:1109605/1‑148 (MQ=255)
                                                                                    ggTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAg                                                                <  1:3184353/148‑1 (MQ=255)
                                                                                    ggTTTGCGGGCGGTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAg                                                                <  2:3144823/148‑1 (MQ=255)
                                                                                     gTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAGCGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGa                                                                         >  1:1820153/1‑138 (MQ=255)
                                                                                      tttGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGAt                                                              <  2:304976/148‑1 (MQ=255)
                                                                                       ttGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGAtt                                                             >  1:1983608/1‑148 (MQ=255)
                                                                                        tGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTa                                                            <  2:1396806/148‑1 (MQ=255)
                                                                                        tGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTa                                                            <  1:246559/148‑1 (MQ=255)
                                                                                         gCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACTCCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTAc                                                           <  2:3299/148‑1 (MQ=255)
                                                                                         gCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTAc                                                           <  2:878729/148‑1 (MQ=255)
                                                                                                ttGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGgcgc                                                    >  1:732629/1‑148 (MQ=255)
                                                                                                ttGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGgcgc                                                    >  1:2857896/1‑148 (MQ=255)
                                                                                                 tGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACTCCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCt                                                   >  1:1005285/1‑148 (MQ=255)
                                                                                                   cgcCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTgg                                                 <  2:1704225/148‑1 (MQ=255)
                                                                                                       cGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTgcg                                             >  1:1145549/1‑148 (MQ=255)
                                                                                                       cGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTgcg                                             >  1:1285993/1‑148 (MQ=255)
                                                                                                            gcgcGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTc                                        <  1:199334/148‑1 (MQ=255)
                                                                                                             cgcgGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCa                                       >  1:1167964/1‑148 (MQ=255)
                                                                                                             cgcgGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCa                                       >  1:1967670/1‑148 (MQ=255)
                                                                                                                 gTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTgg                                   <  1:357604/148‑1 (MQ=255)
                                                                                                                  ttCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAACGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGc                                  >  2:2242630/1‑148 (MQ=255)
                                                                                                                  ttCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGc                                  >  1:2692895/1‑148 (MQ=255)
                                                                                                                   tCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCa                                 <  1:3316703/148‑1 (MQ=255)
                                                                                                                     cTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGt                               <  2:638839/148‑1 (MQ=255)
                                                                                                                     cTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGt                               <  2:3051191/148‑1 (MQ=255)
                                                                                                                     cTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGt                               <  2:2015586/148‑1 (MQ=255)
                                                                                                                       tGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCg                             >  2:1726856/1‑148 (MQ=255)
                                                                                                                       tGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCg                             >  1:2872385/1‑148 (MQ=255)
                                                                                                                           cTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCAtg                         <  2:786903/148‑1 (MQ=255)
                                                                                                                              gTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGa                      >  1:3005856/1‑148 (MQ=255)
                                                                                                                               ttCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGaa                     <  2:286195/148‑1 (MQ=255)
                                                                                                                               ttCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGaa                     <  2:81603/148‑1 (MQ=255)
                                                                                                                                 cATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAaga                   <  2:3436567/148‑1 (MQ=255)
                                                                                                                                  aTCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAagag                  >  2:524539/1‑148 (MQ=255)
                                                                                                                                     gTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGt               <  1:1690494/148‑1 (MQ=255)
                                                                                                                                       cTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTc             >  2:1187659/1‑148 (MQ=255)
                                                                                                                                       cTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTc             >  2:1525319/1‑148 (MQ=255)
                                                                                                                                        ttGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTCa            >  2:1302592/1‑148 (MQ=255)
                                                                                                                                           aCCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTCAAGc         >  2:1349588/1‑148 (MQ=255)
                                                                                                                                               ggACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTCAAGCACCg     >  2:1488016/1‑148 (MQ=255)
                                                                                                                                               ggACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTCAAGCACCg     <  2:3093313/148‑1 (MQ=255)
                                                                                                                                               ggACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTCAAGCACCg     >  1:1407402/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAGATGATCACATCCTGGTGCTGAATAAACCTTCCGGTACGGCGGTACATGGCGGCAGTGGTTTAAGCTTCGGCGTTATTGAAGGTTTGCGGGCGTTGCGCCCGGAAGCGCGGTTCCTTGAACTGGTTCATCGTCTTGACCGGGACACCTCAGGTGTTTTGCTGGTAGCGAAAAAACGCTCGGCGTTGCGTTCTCTGCATGAGCAATTACGTGAAAAAGGGATGCAAAAAGATTACCTGGCGCTGGTGCGCGGTCAGTGGCAGTCGCATGTGAAGAGCGTTCAAGCACCGTTA  >  CP009273/1140687‑1140980

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