Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,255,009 A→T 92.5% 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,0090AT92.5% 553.9 / 23.1 179V256V (GTA→GTTdnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base T (90/75);  total (95/84)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

ACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTG  >  CP009273/4254864‑4255153
                                                                                                                                                 |                                                                                                                                                
acacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCt                                                                                                                                                <  2:1744686/148‑1 (MQ=255)
   ccGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTAt                                                                                                                                             >  1:372988/1‑148 (MQ=255)
      gTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtctt                                                                                                                                          <  2:1199031/148‑1 (MQ=255)
       ttATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttc                                                                                                                                         <  2:1706678/148‑1 (MQ=255)
       ttATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttc                                                                                                                                         <  1:2017709/148‑1 (MQ=255)
          tGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCg                                                                                                                                      <  1:2475178/148‑1 (MQ=255)
          tGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCg                                                                                                                                      >  2:2219066/1‑148 (MQ=255)
          tGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCg                                                                                                                                      <  2:2675338/148‑1 (MQ=255)
          tGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCg                                                                                                                                      <  1:1365295/148‑1 (MQ=255)
            aCGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCt                                                                                                                                    >  2:685959/1‑148 (MQ=255)
              gATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTgg                                                                                                                                  >  2:1893221/1‑148 (MQ=255)
                tctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgag                                                                                                                                <  2:1503433/148‑1 (MQ=255)
                tctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGgag                                                                                                                                <  2:2213991/148‑1 (MQ=255)
                  tcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGAt                                                                                                                              >  1:788108/1‑148 (MQ=255)
                    aacaaaacaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGc                                                                                                                            <  2:3227696/148‑1 (MQ=255)
                     acaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGACAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGcc                                                                                                                           >  2:3031182/1‑148 (MQ=255)
                     acaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGcc                                                                                                                           >  2:2095493/1‑148 (MQ=255)
                     acaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGcc                                                                                                                           >  2:3032826/1‑148 (MQ=255)
                      caaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCa                                                                                                                          >  2:2757239/1‑148 (MQ=255)
                       aaaaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCAt                                                                                                                         >  2:3154042/1‑148 (MQ=255)
                       aaaaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCAt                                                                                                                         <  1:3081106/148‑1 (MQ=255)
                         aaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCa                                                                                                                       >  1:2756077/1‑148 (MQ=255)
                          aaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAg                                                                                                                      >  2:2973736/1‑148 (MQ=255)
                          aaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAg                                                                                                                      >  1:496995/1‑148 (MQ=255)
                          aaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAg                                                                                                                      >  2:2975656/1‑148 (MQ=255)
                          aaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAg                                                                                                                      >  2:272114/1‑148 (MQ=255)
                            aaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGaa                                                                                                                    <  2:1975563/148‑1 (MQ=255)
                            aaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGaa                                                                                                                    >  1:166502/1‑148 (MQ=255)
                             aCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAAc                                                                                                                   >  2:1292124/1‑148 (MQ=255)
                                gCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGa                                                                                                                <  1:2584953/148‑1 (MQ=255)
                                 cTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGAt                                                                                                               >  2:1310292/1‑148 (MQ=255)
                                 cTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGAt                                                                                                               >  1:2616295/1‑148 (MQ=255)
                                  tGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATc                                                                                                              <  2:1956568/148‑1 (MQ=255)
                                  tGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATc                                                                                                              <  2:1956526/148‑1 (MQ=255)
                                   ggCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCa                                                                                                             >  2:1875063/1‑148 (MQ=255)
                                    gCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCat                                                                                                            >  2:2270885/1‑148 (MQ=255)
                                          aGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCg                                                                                                      <  1:1182618/148‑1 (MQ=255)
                                            ccGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGtt                                                                                                    <  2:2325856/148‑1 (MQ=255)
                                               tCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtctc                                                                                                 <  2:592877/148‑1 (MQ=255)
                                               tCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtctc                                                                                                 >  1:2489798/1‑148 (MQ=255)
                                               tCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtatc                                                                                                 >  2:2725329/1‑148 (MQ=255)
                                                cGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTctct                                                                                                <  1:1635752/148‑1 (MQ=255)
                                                 ggATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTg                                                                                               <  1:2851123/148‑1 (MQ=255)
                                                 ggATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTg                                                                                               <  2:2248675/148‑1 (MQ=255)
                                                         tcatcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCt                                                                                       >  2:2460282/1‑148 (MQ=255)
                                                          catcGTCGCCGCGCGTCCCTCGATGGGTAAAACACCATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTg                                                                                      <  2:1951547/148‑1 (MQ=255)
                                                                          ccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGa                                                                      >  1:2060787/1‑148 (MQ=255)
                                                                          ccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGa                                                                      <  1:2253445/148‑1 (MQ=255)
                                                                             tcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTa                                                                   >  2:1658346/1‑148 (MQ=255)
                                                                              cGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaa                                                                  <  1:155329/148‑1 (MQ=255)
                                                                                aTGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTaaaa                                                                >  2:550451/1‑148 (MQ=255)
                                                                                 tGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAAt                                                               <  2:3109535/148‑1 (MQ=255)
                                                                                 tGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAAt                                                               <  2:2049502/148‑1 (MQ=255)
                                                                                 tGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAAt                                                               <  2:2049585/148‑1 (MQ=255)
                                                                                    gTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCg                                                            >  1:1236646/1‑148 (MQ=255)
                                                                                    gTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATACg                                                            >  1:1234361/1‑148 (MQ=255)
                                                                                        aacaacATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTAcc                                                        >  2:945541/1‑148 (MQ=255)
                                                                                            acaTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGc                                                    >  2:1788526/1‑148 (MQ=255)
                                                                                             caTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCa                                                   >  2:402372/1‑148 (MQ=255)
                                                                                                ttGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCt                                                >  1:1179956/1‑148 (MQ=255)
                                                                                                   cGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGa                                             <  2:615752/148‑1 (MQ=255)
                                                                                                     aTGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATg                                           <  1:958722/148‑1 (MQ=255)
                                                                                                     aTGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATg                                           <  1:958165/148‑1 (MQ=255)
                                                                                                     aTGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATg                                           >  1:1504721/1‑148 (MQ=255)
                                                                                                     aTGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATg                                           <  2:3081254/148‑1 (MQ=255)
                                                                                                     aTGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATg                                           <  2:1841588/148‑1 (MQ=255)
                                                                                                         aCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGa                                       >  1:2039076/1‑148 (MQ=255)
                                                                                                          ccTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                      <  1:227367/148‑1 (MQ=255)
                                                                                                          ccTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                      <  1:227237/148‑1 (MQ=255)
                                                                                                          ccTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                      <  2:2851739/148‑1 (MQ=255)
                                                                                                             cgtcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGAc                                   >  1:2766642/1‑148 (MQ=255)
                                                                                                                 gAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTggg                               <  1:516596/148‑1 (MQ=255)
                                                                                                                  aaaaCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgc                              >  1:1077262/1‑148 (MQ=255)
                                                                                                                    aaCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgc                            >  1:677362/1‑148 (MQ=255)
                                                                                                                    aaCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGATGACCAGACTAAAATCCGTACCTGGCAGCACGATGACGAATACTGGgcgc                            >  1:3254921/1‑148 (MQ=255)
                                                                                                                     aCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgcg                           >  1:2037815/1‑148 (MQ=255)
                                                                                                                      cgcgGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgcgc                          >  2:2152221/1‑148 (MQ=255)
                                                                                                                           cgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                     >  2:2679725/1‑148 (MQ=255)
                                                                                                                           cgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                     >  1:3197347/1‑148 (MQ=255)
                                                                                                                             aTGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGAGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTAACGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCg                   >  2:3224177/1‑148 (MQ=255)
                                                                                                                              tGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCgg                  >  1:2323846/1‑148 (MQ=255)
                                                                                                                              tGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCgg                  <  1:1506372/148‑1 (MQ=255)
                                                                                                                              tGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCgg                  <  1:1505920/148‑1 (MQ=255)
                                                                                                                               gTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGc                 >  1:2552915/1‑148 (MQ=255)
                                                                                                                                     ggATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATg           <  1:63023/148‑1 (MQ=255)
                                                                                                                                      gATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATgg          >  2:1762945/1‑148 (MQ=255)
                                                                                                                                        tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGt        <  1:2707375/148‑1 (MQ=255)
                                                                                                                                        tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGt        >  1:2331904/1‑148 (MQ=255)
                                                                                                                                        tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGt        <  1:2460282/148‑1 (MQ=255)
                                                                                                                                        tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGt        >  1:2331655/1‑148 (MQ=255)
                                                                                                                                         aaaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa       <  1:2150148/148‑1 (MQ=255)
                                                                                                                                          aaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAt      <  1:316612/148‑1 (MQ=255)
                                                                                                                                           aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCa             >  2:622631/1‑140 (MQ=255)
                                                                                                                                           aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt     <  1:1805906/148‑1 (MQ=255)
                                                                                                                                           aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt     <  1:2218735/148‑1 (MQ=255)
                                                                                                                                           aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt     <  1:2219066/148‑1 (MQ=255)
                                                                                                                                           aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt     >  1:1075307/1‑148 (MQ=255)
                                                                                                                                            ccGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTACAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt    <  2:2770663/148‑1 (MQ=255)
                                                                                                                                            ccGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt    >  2:1627498/1‑148 (MQ=255)
                                                                                                                                              ggTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTg  <  2:910545/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTG  >  CP009273/4254864‑4255153

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: