Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,255,009 A→T 100% 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,0090AT95.0% 630.7 / 17.0 200V256V (GTA→GTT) dnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (4/6);  new base T (102/88);  total (106/94)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 200 total aligned reads displayed.

ACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCT  >  CP009273/4254864‑4255155
                                                                                                                                                 |                                                                                                                                                  
acacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCt                                                                                                                                                  >  2:2276306/1‑148 (MQ=255)
acacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCt                                                                                                                                                  >  2:1390042/1‑148 (MQ=255)
 cacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCtt                                                                                                                                                 <  2:1917539/148‑1 (MQ=255)
    cGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtc                                                                                                                                              >  1:2319232/1‑148 (MQ=255)
     ggTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtct                                                                                                                                             >  2:1431935/1‑148 (MQ=255)
        tATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttct                                                                                                                                          <  2:565204/148‑1 (MQ=255)
        tATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttct                                                                                                                                          <  2:2357192/148‑1 (MQ=255)
             cGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTg                                                                                                                                     >  1:2644007/1‑148 (MQ=255)
             cGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATATTCTCGCTg                                                                                                                                     <  1:3542932/148‑1 (MQ=255)
              gATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCTAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTgg                                                                                                                                    <  1:2883092/148‑1 (MQ=255)
               aTCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGga                                                                                                                                   >  1:2551645/1‑148 (MQ=255)
                tctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgag                                                                                                                                  <  1:1093223/148‑1 (MQ=255)
                tctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgag                                                                                                                                  >  2:1606664/1‑148 (MQ=255)
                 ctcAACAAACAAACCGCGGGCTTGCAGCCGTCGGATTTGATCCTCGTCGCCGCGCGTCCGTCGATGGGTAAAACTACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgaga                                                                                                                                 <  1:2066555/148‑1 (MQ=255)
                 ctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATACACCGGTTCTTATCTTCTCGCTGgaga                                                                                                                                 >  1:2169712/1‑148 (MQ=255)
                 ctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgaga                                                                                                                                 >  1:2169867/1‑148 (MQ=255)
                 ctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgaga                                                                                                                                 >  1:637478/1‑148 (MQ=255)
                  tcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGAt                                                                                                                                >  2:2011239/1‑148 (MQ=255)
                   caacaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATg                                                                                                                               <  2:791125/148‑1 (MQ=255)
                    aacaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGc                                                                                                                              >  2:1761757/1‑148 (MQ=255)
                      caaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCa                                                                                                                            >  2:1792257/1‑148 (MQ=255)
                       aaaaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCAt                                                                                                                           <  1:1706469/148‑1 (MQ=255)
                         aaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCa                                                                                                                         >  1:195897/1‑148 (MQ=255)
                         aaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCa                                                                                                                         >  1:767677/1‑148 (MQ=255)
                           aaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGa                                                                                                                       <  1:1667233/148‑1 (MQ=255)
                            aaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGaa                                                                                                                      <  1:921772/148‑1 (MQ=255)
                               cGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAg                                                                                                                   <  1:1768194/148‑1 (MQ=255)
                                     cTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatg                                                                                                             <  1:3043623/148‑1 (MQ=255)
                                      ttGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatga                                                                                                            <  1:606713/148‑1 (MQ=255)
                                       tGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatgat                                                                                                           >  2:2003186/1‑148 (MQ=255)
                                        gCAGCCGTAGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatgatg                                                                                                          <  1:3239250/148‑1 (MQ=255)
                                            ccGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGtt                                                                                                      >  1:2256880/1‑148 (MQ=255)
                                             cGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtc                                                                                                     <  2:3132787/148‑1 (MQ=255)
                                               tCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtctc                                                                                                   >  1:762569/1‑148 (MQ=255)
                                               tCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtctc                                                                                                   >  2:1117267/1‑148 (MQ=255)
                                                cGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTctct                                                                                                  >  2:2505341/1‑148 (MQ=255)
                                                cGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTctct                                                                                                  <  2:2845817/148‑1 (MQ=255)
                                                cGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTctct                                                                                                  >  2:3053037/1‑148 (MQ=255)
                                                cGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTctct                                                                                                  <  1:113132/148‑1 (MQ=255)
                                                 ggATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTg                                                                                                 <  2:3524239/148‑1 (MQ=255)
                                                    tttGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCg                                                                                              >  1:2404546/1‑148 (MQ=255)
                                                       gATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCg                                                                                           <  2:767677/148‑1 (MQ=255)
                                                          catcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTg                                                                                        <  2:2086306/148‑1 (MQ=255)
                                                             cgtcgCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTcg                                                                                     >  1:1928128/1‑148 (MQ=255)
                                                                 gccgcGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTcgcgcg                                                                                 <  2:3167217/148‑1 (MQ=255)
                                                                      gcgTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGAc                                                                            <  1:2650910/148‑1 (MQ=255)
                                                                       cgtccgtcGATGGGTAAAACAACATTTGCTATGAAACTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTAGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGAcc                                                                           >  2:2251109/1‑148 (MQ=255)
                                                                         tccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAg                                                                         >  2:1362601/1‑148 (MQ=255)
                                                                              cGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaa                                                                    >  1:1402556/1‑148 (MQ=255)
                                                                               gATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaaa                                                                   <  2:1833411/148‑1 (MQ=255)
                                                                               gATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaaa                                                                   <  2:2404546/148‑1 (MQ=255)
                                                                                aTGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaaaa                                                                  >  2:685157/1‑148 (MQ=255)
                                                                                aTGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaaaa                                                                  >  2:92961/1‑148 (MQ=255)
                                                                                aTGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaaaa                                                                  <  1:1502003/148‑1 (MQ=255)
                                                                                aTGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaaaa                                                                  >  2:93515/1‑148 (MQ=255)
                                                                                    gTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCg                                                              >  1:397416/1‑148 (MQ=255)
                                                                                    gTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCg                                                              >  1:397360/1‑148 (MQ=255)
                                                                                    gTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCg                                                              >  1:2784671/1‑148 (MQ=255)
                                                                                              aTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAg                                                    <  1:2374745/148‑1 (MQ=255)
                                                                                                ttGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCt                                                  <  1:2956816/148‑1 (MQ=255)
                                                                                                   cGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATTCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGa                                               <  1:3544890/148‑1 (MQ=255)
                                                                                                   cGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGa                                               <  1:1042183/148‑1 (MQ=255)
                                                                                                   cGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGa                                               <  1:2017429/148‑1 (MQ=255)
                                                                                                       gAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGAc                                           >  1:2777510/1‑148 (MQ=255)
                                                                                                       gAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGAc                                           >  1:2861707/1‑148 (MQ=255)
                                                                                                       gAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGAc                                           >  2:2125119/1‑148 (MQ=255)
                                                                                                        aaCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACg                                          >  1:2882908/1‑148 (MQ=255)
                                                                                                        aaCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACg                                          >  1:895315/1‑148 (MQ=255)
                                                                                                        aaCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACg                                          >  1:2051949/1‑148 (MQ=255)
                                                                                                         aCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGa                                         >  1:1437498/1‑148 (MQ=255)
                                                                                                          ccTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                        >  2:3338195/1‑148 (MQ=255)
                                                                                                           cTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAg                                       >  2:2528547/1‑148 (MQ=255)
                                                                                                            tcgtcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGa                                      <  1:417537/148‑1 (MQ=255)
                                                                                                             cgtcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGAc                                     >  2:2720281/1‑148 (MQ=255)
                                                                                                              gtcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACt                                    <  2:1189099/148‑1 (MQ=255)
                                                                                                               tcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTg                                   >  1:828400/1‑148 (MQ=255)
                                                                                                                cgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTgg                                  >  1:609843/1‑148 (MQ=255)
                                                                                                                cgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTgg                                  >  1:1522720/1‑148 (MQ=255)
                                                                                                                cgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTgg                                  <  1:1571090/148‑1 (MQ=255)
                                                                                                                  aaaaCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgc                                >  2:1373662/1‑148 (MQ=255)
                                                                                                                        cggcgATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCAt                          <  2:2487200/148‑1 (MQ=255)
                                                                                                                          gcgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCAttt                        <  1:1329743/148‑1 (MQ=255)
                                                                                                                           cgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                       <  1:557719/148‑1 (MQ=255)
                                                                                                                           cgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                       >  2:1758455/1‑148 (MQ=255)
                                                                                                                           cgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                       >  1:2352866/1‑148 (MQ=255)
                                                                                                                             aTGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCg                     >  2:2815280/1‑148 (MQ=255)
                                                                                                                              tGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCgg                    >  2:1235684/1‑148 (MQ=255)
                                                                                                                              tGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCgg                    >  2:2614159/1‑148 (MQ=255)
                                                                                                                                      gATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATgg            <  2:2578640/148‑1 (MQ=255)
                                                                                                                                        tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGt          >  2:369275/1‑148 (MQ=255)
                                                                                                                                        tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGt          <  1:519873/148‑1 (MQ=255)
                                                                                                                                         aaaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa         >  1:1636707/1‑148 (MQ=255)
                                                                                                                                         aaaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa         >  2:2339541/1‑148 (MQ=255)
                                                                                                                                           aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt       <  1:685157/148‑1 (MQ=255)
                                                                                                                                           aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt       >  1:2096647/1‑148 (MQ=255)
                                                                                                                                           aCCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt       <  2:262974/148‑1 (MQ=255)
                                                                                                                                            ccGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt      >  2:2743079/1‑148 (MQ=255)
                                                                                                                                            ccGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGCTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt      <  2:2644007/148‑1 (MQ=255)
                                                                                                                                             cGGTTCTTATCTTCTCTCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtttt     >  1:339867/1‑148 (MQ=255)
                                                                                                                                                tACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCt  <  1:2298248/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGCT  >  CP009273/4254864‑4255155

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