Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,255,009 A→T 85.9% V256V (GTA→GTT)  dnaB → replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,255,0090AT85.9% 917.7 / 120.1 318V256V (GTA→GTT) dnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base T (146/127);  total (164/154)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 318 total aligned reads displayed.

AAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTC  >  CP009273/4254862‑4255156
                                                                                                                                                   |                                                                                                                                                   
 aaCACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCt                                                                                                                                                   >  1:402636/1‑149 (MQ=255)
  acacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCtt                                                                                                                                                  >  1:4259212/1‑149 (MQ=255)
  acacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACtt                                                                                                                                                  <  2:4341292/149‑1 (MQ=255)
   ccccGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTa                                                                                                                                                 >  1:4318629/3‑149 (MQ=255)
   cacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTa                                                                                                                                                 <  2:2204412/149‑1 (MQ=255)
    acCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAt                                                                                                                                                >  1:1227245/1‑149 (MQ=255)
     ccGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtc                                                                                                                                               <  1:1335597/149‑1 (MQ=255)
        gTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttc                                                                                                                                            <  1:3447626/149‑1 (MQ=255)
         ttATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTAtcttct                                                                                                                                           >  1:3239741/1‑149 (MQ=255)
          tATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTtctc                                                                                                                                          <  1:3550538/149‑1 (MQ=255)
             gACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTCTCTTCTCGCt                                                                                                                                       >  2:3520670/1‑149 (MQ=255)
             gACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCt                                                                                                                                       <  1:3089328/149‑1 (MQ=255)
              aagaactcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTTCGATGAACCTCGTCGAAAAAGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTg                                                                                                                                      >  2:857641/6‑149 (MQ=255)
                gATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGga                                                                                                                                    >  1:622630/1‑149 (MQ=255)
                  tctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgaga                                                                                                                                  <  2:1123877/149‑1 (MQ=255)
                   ctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGAt                                                                                                                                 >  1:1461885/1‑149 (MQ=255)
                      aacaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGcc                                                                                                                              >  1:1253804/1‑149 (MQ=255)
                         aaaaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGGTCTTATATTCTCGCTGGTGAAGCCCTc                                                                                                                           >  2:749089/1‑149 (MQ=255)
                          aaaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCa                                                                                                                          >  2:303799/1‑149 (MQ=255)
                             aaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGaa                                                                                                                       <  1:2397071/149‑1 (MQ=255)
                              aaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAAc                                                                                                                      <  1:1609493/149‑1 (MQ=255)
                                   cTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATc                                                                                                                 <  2:2692758/149‑1 (MQ=255)
                                   cTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATc                                                                                                                 <  2:2870350/149‑1 (MQ=255)
                                    tGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCa                                                                                                                >  2:2659488/1‑149 (MQ=255)
                                    tGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCa                                                                                                                >  2:4744096/1‑149 (MQ=255)
                                     ggCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCat                                                                                                               <  2:381854/149‑1 (MQ=255)
                                     ggCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCat                                                                                                               <  2:198357/149‑1 (MQ=255)
                                      gCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatg                                                                                                              >  1:3670050/1‑149 (MQ=255)
                                      gCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatg                                                                                                              >  1:781257/1‑149 (MQ=255)
                                          gCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGc                                                                                                          <  1:1190529/149‑1 (MQ=255)
                                           cAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCg                                                                                                         <  1:923006/149‑1 (MQ=255)
                                            aGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGt                                                                                                        <  1:3932762/149‑1 (MQ=255)
                                             gCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGtt                                                                                                       <  1:520437/149‑1 (MQ=255)
                                               cGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtct                                                                                                     <  1:1487592/149‑1 (MQ=255)
                                                gTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtctc                                                                                                    <  2:667877/149‑1 (MQ=255)
                                                gTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtctc                                                                                                    <  2:4650863/149‑1 (MQ=255)
                                                gTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtctc                                                                                                    <  1:250339/149‑1 (MQ=255)
                                                 tCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTctct                                                                                                   >  1:2753673/1‑149 (MQ=255)
                                                 tCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTctct                                                                                                   >  2:2460496/1‑149 (MQ=255)
                                                    gATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATCAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGc                                                                                                >  1:2595329/1‑149 (MQ=255)
                                                    gATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGc                                                                                                >  2:327380/1‑149 (MQ=255)
                                                    gATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGc                                                                                                >  2:2353824/1‑149 (MQ=255)
                                                    gATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGc                                                                                                >  2:1964325/1‑149 (MQ=255)
                                                    gATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGc                                                                                                >  2:303190/1‑149 (MQ=255)
                                                         gATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGc                                                                                           >  1:2289956/1‑149 (MQ=255)
                                                          atcatcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCt                                                                                          <  2:1852668/149‑1 (MQ=255)
                                                            catcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGt                                                                                        >  1:1001567/1‑149 (MQ=255)
                                                            catcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGt                                                                                        <  2:2160474/149‑1 (MQ=255)
                                                             atcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTc                                                                                       <  2:2200517/149‑1 (MQ=255)
                                                              tcgtcgCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTcg                                                                                      <  2:274715/149‑1 (MQ=255)
                                                                   gccgcGCGTCCGTCGATGTGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGt                                                                                 <  2:2617148/149‑1 (MQ=255)
                                                                      gcgcgTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGa                                                                              >  2:12132/1‑149 (MQ=255)
                                                                       cgcgTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGAc                                                                             <  1:4363592/149‑1 (MQ=255)
                                                                         cgtccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCa                                                                           <  1:4768875/149‑1 (MQ=255)
                                                                           tccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGa                                                                         >  2:101896/1‑149 (MQ=255)
                                                                             cgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACt                                                                       >  2:4581474/1‑149 (MQ=255)
                                                                               tcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaa                                                                     <  1:1539457/149‑1 (MQ=255)
                                                                                  aTGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAAt                                                                  <  2:451707/149‑1 (MQ=255)
                                                                                   tGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATc                                                                 >  2:2491572/1‑149 (MQ=255)
                                                                                    gggTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATcc                                                                >  1:2186019/1‑149 (MQ=255)
                                                                                    gggTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATcc                                                                >  2:3386387/1‑149 (MQ=255)
                                                                                     ggTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCg                                                               <  1:4546394/149‑1 (MQ=255)
                                                                                        aaaaCAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTAc                                                            >  1:3267311/1‑149 (MQ=255)
                                                                                          aacaacATTTGCGATGAACCTCGTCGAAAACTCGGCGATGTTGCAGGATAAACCGGTAATTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCTTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCg                                                          >  2:4358293/1‑149 (MQ=255)
                                                                                          aacaacATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCg                                                          >  2:593102/1‑149 (MQ=255)
                                                                                           acaacaTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCgg                                                         <  1:81273/149‑1 (MQ=255)
                                                                                            caacaTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCggg                                                        >  1:3788235/1‑149 (MQ=255)
                                                                                              acaTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCa                                                      >  2:1309939/1‑149 (MQ=255)
                                                                                                 tttGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCt                                                   <  2:1384071/149‑1 (MQ=255)
                                                                                                    gCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCCGCTCGa                                                >  2:3654921/1‑149 (MQ=255)
                                                                                                      gATGAACCTCGTCGAAAACGCGGCGTTGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATg                                              <  2:2089730/149‑1 (MQ=255)
                                                                                                      gATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATg                                              >  1:2870487/1‑149 (MQ=255)
                                                                                                       aTGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGa                                             >  2:3248724/1‑149 (MQ=255)
                                                                                                       aTGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGa                                             <  1:1006719/149‑1 (MQ=255)
                                                                                                        tGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGAc                                            >  1:1065006/1‑149 (MQ=255)
                                                                                                           aCCTCGTCGAAAACGCGGCTATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                         >  1:4101950/1‑149 (MQ=255)
                                                                                                           aCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                         >  2:1076913/1‑149 (MQ=255)
                                                                                                            ccTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAg                                        <  2:1298203/149‑1 (MQ=255)
                                                                                                             cTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGa                                       >  1:677410/1‑149 (MQ=255)
                                                                                                                gtcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTg                                    <  1:3020454/149‑1 (MQ=255)
                                                                                                                gtcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTg                                    >  2:2134965/1‑149 (MQ=255)
                                                                                                                 tcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCGGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTgg                                   <  2:4224330/149‑1 (MQ=255)
                                                                                                                    aaaaCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcg                                >  1:1871978/1‑149 (MQ=255)
                                                                                                                         gcggcgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCAt                           <  2:4595173/149‑1 (MQ=255)
                                                                                                                         gcggcgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCAt                           <  2:1767472/149‑1 (MQ=255)
                                                                                                                            gcgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                        >  2:3825722/1‑149 (MQ=255)
                                                                                                                              gATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCg                      <  2:1514271/149‑1 (MQ=255)
                                                                                                                              gATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCg                      >  1:956293/1‑149 (MQ=255)
                                                                                                                                     cAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCAt               >  2:4425422/1‑149 (MQ=255)
                                                                                                                                       ggATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATgg             <  2:4657689/149‑1 (MQ=255)
                                                                                                                                       ggATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATgg             <  1:3114953/149‑1 (MQ=255)
                                                                                                                                          tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa          >  2:520652/1‑149 (MQ=255)
                                                                                                                                          tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa          <  1:4089991/149‑1 (MQ=255)
                                                                                                                                           aaaCCGGTTCTTATCTTCTCGCTGGAGATACCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAt         <  1:731240/149‑1 (MQ=255)
                                                                                                                                            aaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt        >  1:2135343/1‑149 (MQ=255)
                                                                                                                                             aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt       >  2:4052867/1‑149 (MQ=255)
                                                                                                                                              ccGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtttt      >  2:1532050/1‑149 (MQ=255)
                                                                                                                                               cGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTg     <  2:169128/149‑1 (MQ=255)
                                                                                                                                               cGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACT‑AAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtttttc    >  2:3631789/1‑147 (MQ=255)
                                                                                                                                                ggTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGc    <  1:2216991/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCTC  >  CP009273/4254862‑4255156

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: ‑