Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 61,967 A→C 19.3% L433R (CTG→CGG)  rapA ← RNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927361,9670AC19.3% 351.3 / 113.8 192L433R (CTG→CGG) rapARNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase
Reads supporting (aligned to +/- strand):  ref base A (84/71);  new base C (21/16);  total (105/87)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 193 total aligned reads displayed.

GCGTTATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAACA  >  CP009273/61820‑62111
                                                                                                                                                   |                                                                                                                                                
gCGTTATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCAATAATGCAGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCa                                                                                                                                                  >  1:1033288/1‑148 (MQ=255)
  gTTATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGc                                                                                                                                                <  2:265699/148‑1 (MQ=255)
  gTTATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGc                                                                                                                                                <  2:2194438/148‑1 (MQ=255)
  gTTATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGc                                                                                                                                                <  2:667417/148‑1 (MQ=255)
    tATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTc                                                                                                                                              >  1:1181551/1‑148 (MQ=255)
    tATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGACCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTc                                                                                                                                              >  1:3087273/1‑148 (MQ=255)
     aTCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTcg                                                                                                                                             >  1:1543699/1‑148 (MQ=255)
         ccTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGtt                                                                                                                                         >  2:2726454/1‑148 (MQ=255)
         ccTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGtt                                                                                                                                         >  2:269954/1‑148 (MQ=255)
                 ttCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAAt                                                                                                                                 >  1:2660327/1‑148 (MQ=255)
                    cTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCt                                                                                                                              >  1:2410851/1‑148 (MQ=255)
                     tGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCtt                                                                                                                             >  1:2396418/1‑148 (MQ=255)
                         aaaTACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTcac                                                                                                                         <  1:3121024/148‑1 (MQ=255)
                          aaTACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTcaca                                                                                                                        >  1:2189857/1‑148 (MQ=255)
                           aTACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCacac                                                                                                                       >  2:2503092/1‑148 (MQ=255)
                           aTACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTGTCGGGAATCCTTTCacac                                                                                                                       >  1:2631974/1‑148 (MQ=255)
                            tACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACAcc                                                                                                                      >  2:904096/1‑148 (MQ=255)
                            tACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACAcc                                                                                                                      >  1:2224907/1‑148 (MQ=255)
                            tACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACAcc                                                                                                                      >  1:3627287/1‑148 (MQ=255)
                            tACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACAcc                                                                                                                      >  1:2646856/1‑148 (MQ=255)
                             aCGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCg                                                                                                                     <  1:2419408/148‑1 (MQ=255)
                             aCGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCg                                                                                                                     >  1:2604052/1‑148 (MQ=255)
                             aCGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCg                                                                                                                     <  2:3690040/148‑1 (MQ=255)
                                cTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTa                                                                                                                  <  2:3164760/148‑1 (MQ=255)
                                  ccGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTAcg                                                                                                                >  2:3437783/1‑148 (MQ=255)
                                    gggTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTAcgcg                                                                                                              <  1:130908/148‑1 (MQ=255)
                                         gagCATCTCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCAGAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTa                                                                                                         <  2:3605157/148‑1 (MQ=255)
                                          agCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTAc                                                                                                        >  2:2705519/1‑148 (MQ=255)
                                          agCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTAc                                                                                                        >  1:639172/1‑148 (MQ=255)
                                           gCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACg                                                                                                       >  1:2996141/1‑148 (MQ=255)
                                              tatCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGaa                                                                                                    <  2:3160439/148‑1 (MQ=255)
                                               atCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAAc                                                                                                   >  1:2021818/1‑148 (MQ=255)
                                                tCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACa                                                                                                  <  2:3299130/148‑1 (MQ=255)
                                                 cgcgAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAg                                                                                                 >  2:1487646/1‑148 (MQ=255)
                                                       cACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCAcgcg                                                                                           >  2:3252096/1‑148 (MQ=255)
                                                       cACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCAcgcg                                                                                           >  1:2848352/1‑148 (MQ=255)
                                                       cACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCAcgcg                                                                                           >  1:936750/1‑148 (MQ=255)
                                                         cGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGc                                                                                         <  2:507566/148‑1 (MQ=255)
                                                           aTCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTg                                                                                       >  1:245028/1‑148 (MQ=255)
                                                            tCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTgg                                                                                      <  2:2786509/148‑1 (MQ=255)
                                                            tCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTgg                                                                                      >  2:1045170/1‑148 (MQ=255)
                                                            tCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTGTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTgg                                                                                      >  2:3656765/1‑148 (MQ=255)
                                                             cTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGt                                                                                     >  2:3757295/1‑148 (MQ=255)
                                                               tCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGc                                                                                   <  2:114764/148‑1 (MQ=255)
                                                                 cGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCg                                                                                 <  1:2260216/148‑1 (MQ=255)
                                                                 cGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCg                                                                                 <  2:3632209/148‑1 (MQ=255)
                                                                 cGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCg                                                                                 <  1:2794604/148‑1 (MQ=255)
                                                                   cACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTg                                                                               >  2:630910/1‑148 (MQ=255)
                                                                       tttACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGa                                                                           <  1:3767027/148‑1 (MQ=255)
                                                                        ttACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGAt                                                                          <  1:3437783/148‑1 (MQ=255)
                                                                         tACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATc                                                                         <  1:707688/148‑1 (MQ=255)
                                                                          aCGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATcc                                                                        <  1:3818418/148‑1 (MQ=255)
                                                                                 cccATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCa                                                                 >  2:2210394/1‑148 (MQ=255)
                                                                                  ccATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCAt                                                                >  2:3189878/1‑148 (MQ=255)
                                                                                  ccATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCAt                                                                >  2:1785383/1‑148 (MQ=255)
                                                                                     tAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGa                                                             >  2:1396563/1‑148 (MQ=255)
                                                                                      aaTGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGaa                                                            <  2:1878027/148‑1 (MQ=255)
                                                                                       aTGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGaaa                                                           >  2:3093513/1‑148 (MQ=255)
                                                                                       aTGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGaaa                                                           <  1:1247247/148‑1 (MQ=255)
                                                                                        tGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAAc                                                          <  1:1259012/148‑1 (MQ=255)
                                                                                           cGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAg                                                       <  2:2742971/148‑1 (MQ=255)
                                                                                            ggAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTTCAGCTCGCGTTTCGGGAATCCTGTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGc                                                      >  2:832526/1‑148 (MQ=255)
                                                                                             gagaCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCt                                                     <  2:2415479/148‑1 (MQ=255)
                                                                                              agaCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTc                                                    >  1:596744/1‑148 (MQ=255)
                                                                                              agaCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGGGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTc                                                    >  1:1611987/1‑148 (MQ=255)
                                                                                                aCTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCt                                                  >  2:3558217/1‑148 (MQ=255)
                                                                                                 cTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTg                                                 <  1:1106242/148‑1 (MQ=255)
                                                                                                 cTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTg                                                 <  2:1225741/148‑1 (MQ=255)
                                                                                                 cTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTg                                                 <  1:239440/148‑1 (MQ=255)
                                                                                                   ttAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGAc                                               >  2:1673649/1‑148 (MQ=255)
                                                                                                        aGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGgcg                                          <  1:1847281/148‑1 (MQ=255)
                                                                                                        aGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGgcg                                          <  1:3650858/148‑1 (MQ=255)
                                                                                                             tCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGTCGGGCGCTCTg                                     >  2:2183836/1‑148 (MQ=255)
                                                                                                             tCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTg                                     <  2:1946605/148‑1 (MQ=255)
                                                                                                              cTGATACTGCGTCGGTAGCGGAAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTgg                                    >  2:3567127/1‑148 (MQ=255)
                                                                                                               tGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTggg                                   <  2:2983403/148‑1 (MQ=255)
                                                                                                                 aTACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCa                                 <  1:1001360/148‑1 (MQ=255)
                                                                                                                   aCTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATc                               >  1:3176661/1‑148 (MQ=255)
                                                                                                                    cTGCGTCGGTAGCTGCAGCTTAATGGTGTGCAGCTCGCGTGTCGGTAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGAGATCCATCAGCATCGAAACCAGGTCCTGACGGGCGCTCTGGGCATCt                              <  2:2755576/148‑1 (MQ=255)
                                                                                                                     tGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCtt                             <  1:630910/148‑1 (MQ=255)
                                                                                                                       cGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtcg                           <  1:1399035/148‑1 (MQ=255)
                                                                                                                          cGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtcgctg                        <  2:2145422/148‑1 (MQ=255)
                                                                                                                               gCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTtcgctgtcgct                   >  1:2888564/1‑148 (MQ=255)
                                                                                                                                 ggCAGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTcgctgtcgctgt                 <  1:3252096/148‑1 (MQ=255)
                                                                                                                                 ggCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTcgctgtcgctgt                 <  1:269954/148‑1 (MQ=255)
                                                                                                                                   cAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGttt               >  2:1921795/1‑148 (MQ=255)
                                                                                                                                    aGCTTAATGGTGTGCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtg              <  1:1487646/148‑1 (MQ=255)
                                                                                                                                     gCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgc             >  1:1720898/1‑148 (MQ=255)
                                                                                                                                      cTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgct            >  2:1294630/1‑148 (MQ=255)
                                                                                                                                      cTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgct            >  2:2022902/1‑148 (MQ=255)
                                                                                                                                      cTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgct            >  2:2857163/1‑148 (MQ=255)
                                                                                                                                       ttAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgctg           <  2:995434/148‑1 (MQ=255)
                                                                                                                                        tAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTtgctgc          <  2:1146394/148‑1 (MQ=255)
                                                                                                                                         aaTGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACTGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcc         <  2:3586855/148‑1 (MQ=255)
                                                                                                                                           tGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcctg       >  2:1431442/1‑148 (MQ=255)
                                                                                                                                           tGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcctg       >  2:292849/1‑148 (MQ=255)
                                                                                                                                              tgtgCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGctgcctgcca    >  2:1198676/1‑146 (MQ=255)
                                                                                                                                               gtgCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAAc   >  2:489541/1‑148 (MQ=255)
                                                                                                                                                tgCCGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAACa  <  2:691638/148‑1 (MQ=255)
                                                                                                                                                tgCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAACa  <  1:2503092/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGTTATCACCTTCAAATTCCTGATAAATACGCTCCGGGTAGAGCATATCGCGAGCACGATCTTCCGCACTTTTACGTGCGCCCATAATGCCGGAGACTTTAATAGCCGTCTGATACTGCGTCGGTAGCGGCAGCTTAATGGTGTGCAGCTCGCGTTTCGGGAATCCTTTCACACCGTTACGCGTGTTACGGAACAGCACGCGGCTGGTGCCGTGGCGATCCATCAGCATCGAAACCAGCTCCTGACGGGCGCTCTGGGCATCTTCGCTGTCGCTGTTTGCTGCCTGCAACA  >  CP009273/61820‑62111

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: