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,0090AT98.8% 839.2 / ‑0.6 252V256V (GTA→GTTdnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (143/106);  total (144/108)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 252 total aligned reads displayed.

AACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGC  >  CP009273/4254863‑4255154
                                                                                                                                                  |                                                                                                                                                 
aaCACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCt                                                                                                                                                 >  2:18388/1‑149 (MQ=255)
 acacCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACtt                                                                                                                                                >  1:2115144/1‑149 (MQ=255)
    ccGTTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtc                                                                                                                                             <  1:2196041/149‑1 (MQ=255)
    ccGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAAAAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtc                                                                                                                                             >  1:2985549/1‑149 (MQ=255)
     cGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtct                                                                                                                                            <  2:3674627/149‑1 (MQ=255)
       gTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttc                                                                                                                                          <  1:3409391/149‑1 (MQ=255)
       gTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttc                                                                                                                                          >  2:159307/1‑149 (MQ=255)
       gTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttc                                                                                                                                          <  1:3001386/149‑1 (MQ=255)
        ttATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTAtcttct                                                                                                                                         <  1:626047/149‑1 (MQ=255)
          aTGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCg                                                                                                                                       >  1:1591049/1‑149 (MQ=255)
           tGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGc                                                                                                                                      <  1:240883/149‑1 (MQ=255)
             aCGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTg                                                                                                                                    >  2:2568039/1‑149 (MQ=255)
             aCGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTg                                                                                                                                    >  1:1413008/1‑149 (MQ=255)
              cGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTgg                                                                                                                                   >  2:4022998/1‑149 (MQ=255)
              cGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTAGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTgg                                                                                                                                   >  2:3283421/1‑149 (MQ=255)
                  ctcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGAt                                                                                                                               <  2:1232648/149‑1 (MQ=255)
                   tcAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATg                                                                                                                              >  1:1037137/1‑149 (MQ=255)
                    caacaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGc                                                                                                                             <  2:2436227/149‑1 (MQ=255)
                     aacaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGcc                                                                                                                            >  1:3919894/1‑149 (MQ=255)
                     aacaaAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGcc                                                                                                                            <  1:845961/149‑1 (MQ=255)
                        aaaaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATc                                                                                                                         >  2:1710418/1‑149 (MQ=255)
                         aaaaaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCa                                                                                                                        <  1:3315147/149‑1 (MQ=255)
                            aaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGaa                                                                                                                     <  2:3535475/149‑1 (MQ=255)
                            aaaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGaa                                                                                                                     <  2:3535332/149‑1 (MQ=255)
                             aaCCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAAc                                                                                                                    >  2:3382890/1‑149 (MQ=255)
                                 gCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGAt                                                                                                                <  2:58057/149‑1 (MQ=255)
                                 gCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGAt                                                                                                                <  2:2221585/149‑1 (MQ=255)
                                 gCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGAt                                                                                                                >  2:666774/1‑149 (MQ=255)
                                  cTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACGTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATc                                                                                                               >  1:1681516/1‑149 (MQ=255)
                                  cTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCTCTGGAGATGCCATCAGAACAGAta                                                                                                               >  1:2110176/1‑148 (MQ=255)
                                    ggCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCat                                                                                                             >  2:1456017/1‑149 (MQ=255)
                                    ggCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCat                                                                                                             >  2:1986096/1‑149 (MQ=255)
                                     gCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatg                                                                                                            >  2:3701528/1‑149 (MQ=255)
                                     gCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatg                                                                                                            >  2:1753093/1‑149 (MQ=255)
                                     gCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatg                                                                                                            >  2:3795617/1‑149 (MQ=255)
                                      cTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatga                                                                                                           <  2:1937043/149‑1 (MQ=255)
                                        tGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatgatg                                                                                                         <  1:1944723/149‑1 (MQ=255)
                                        tGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCatgatg                                                                                                         <  1:764280/149‑1 (MQ=255)
                                         gCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGc                                                                                                        <  1:3827321/149‑1 (MQ=255)
                                             ccGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtc                                                                                                    <  1:1291721/149‑1 (MQ=255)
                                              cGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTtct                                                                                                   <  1:3785641/149‑1 (MQ=255)
                                                     tttGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGt                                                                                            >  2:1717919/1‑149 (MQ=255)
                                                     tttGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGt                                                                                            >  1:3237138/1‑149 (MQ=255)
                                                      ttGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTc                                                                                           >  1:222232/1‑149 (MQ=255)
                                                      ttGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTc                                                                                           >  1:2592564/1‑149 (MQ=255)
                                                        gATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGc                                                                                         >  1:580414/1‑149 (MQ=255)
                                                         atcatcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCt                                                                                        >  2:3011563/1‑149 (MQ=255)
                                                           catcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGt                                                                                      <  1:1747973/149‑1 (MQ=255)
                                                            atcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTc                                                                                     <  1:1178391/149‑1 (MQ=255)
                                                            atcGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTc                                                                                     >  1:1380853/1‑149 (MQ=255)
                                                                  gccgcGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGt                                                                               <  1:3775889/149‑1 (MQ=255)
                                                                       gcgTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGAcc                                                                          >  2:468712/1‑149 (MQ=255)
                                                                        cgtccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGAGATGTTTCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCa                                                                         <  2:3094875/149‑1 (MQ=255)
                                                                         gtccgtcGATGGTTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAg                                                                        >  1:3717476/1‑149 (MQ=255)
                                                                         gtccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAg                                                                        >  1:3719201/1‑149 (MQ=255)
                                                                         gtccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAg                                                                        >  2:1589227/1‑149 (MQ=255)
                                                                          tccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGa                                                                       <  2:1662081/149‑1 (MQ=255)
                                                                           ccgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGAc                                                                      <  1:162141/149‑1 (MQ=255)
                                                                            cgtcGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACt                                                                     <  1:3961060/149‑1 (MQ=255)
                                                                                gATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGATCGTTGACCAGACTaaaa                                                                 >  2:1793852/1‑149 (MQ=255)
                                                                                   gggTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATcc                                                              >  1:787155/1‑149 (MQ=255)
                                                                                   gggTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATcc                                                              >  2:100190/1‑149 (MQ=255)
                                                                                     gTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGt                                                            >  1:3587698/1‑149 (MQ=255)
                                                                                       aaaaCAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTAc                                                          <  1:887341/149‑1 (MQ=255)
                                                                                       aaaaCAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTAc                                                          <  1:3194800/149‑1 (MQ=255)
                                                                                          acaacaTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCgg                                                       >  1:447200/1‑149 (MQ=255)
                                                                                              caTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAg                                                   >  2:360537/1‑149 (MQ=255)
                                                                                               aTTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGc                                                  >  2:2619873/1‑149 (MQ=255)
                                                                                                tttGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCt                                                 >  2:2292038/1‑149 (MQ=255)
                                                                                                tttGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCt                                                 >  2:3417442/1‑149 (MQ=255)
                                                                                                  tGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCg                                               <  2:633281/149‑1 (MQ=255)
                                                                                                         aaCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGa                                        >  1:461277/1‑149 (MQ=255)
                                                                                                         aaCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGa                                        >  1:2037323/1‑149 (MQ=255)
                                                                                                          aCCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGaa                                       >  2:1314582/1‑149 (MQ=255)
                                                                                                             tcgtcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGAc                                    <  1:3806759/149‑1 (MQ=255)
                                                                                                                tcgAAAACGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTgg                                 >  2:2792400/1‑149 (MQ=255)
                                                                                                                   aaaaCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcg                              >  1:3023451/1‑149 (MQ=255)
                                                                                                                     aaCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgcg                            >  2:531766/1‑149 (MQ=255)
                                                                                                                      aCGCGGCGATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgcgc                           >  2:798093/1‑149 (MQ=255)
                                                                                                                      aCGCGGCGATGTTGCAGGATAAACCGGATCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGgcgcgc                           >  2:122949/1‑149 (MQ=255)
                                                                                                                           gcgATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTc                      <  1:2547222/149‑1 (MQ=255)
                                                                                                                             gATGTTGCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCg                    >  2:2063727/1‑149 (MQ=255)
                                                                                                                                   gCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCa              <  1:180706/149‑1 (MQ=255)
                                                                                                                                   gCAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCa              <  2:802830/149‑1 (MQ=255)
                                                                                                                                    cAGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCAt             <  2:1858459/149‑1 (MQ=255)
                                                                                                                                     aGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATg            <  1:4090934/149‑1 (MQ=255)
                                                                                                                                     aGGATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATg            >  1:1565096/1‑149 (MQ=255)
                                                                                                                                       gATAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATggg          >  1:4078673/1‑149 (MQ=255)
                                                                                                                                       gATAAACCGGTTCTTATATTCTCGCTGGAGATGCCATCAGAACAGATCATGATTCGTTCTCTGGCGTCGCTTTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACAATggg          >  2:3973062/1‑149 (MQ=255)
                                                                                                                                         tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCTCATTTCCGGCACCATGGGTa        <  2:4089567/149‑1 (MQ=255)
                                                                                                                                         tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa        >  2:2149589/1‑149 (MQ=255)
                                                                                                                                         tAAACCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTa        >  2:2312115/1‑149 (MQ=255)
                                                                                                                                          aaaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAt       >  2:3216730/1‑149 (MQ=255)
                                                                                                                                           aaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATTATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCAGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt      >  2:1999129/1‑149 (MQ=255)
                                                                                                                                           aaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt      >  2:1336332/1‑149 (MQ=255)
                                                                                                                                           aaCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAtt      >  1:1299268/1‑149 (MQ=255)
                                                                                                                                            aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt     >  2:896858/1‑149 (MQ=255)
                                                                                                                                            aCCGGTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTAttt     <  1:358302/149‑1 (MQ=255)
                                                                                                                                               ggTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGc  <  1:2697984/149‑1 (MQ=255)
                                                                                                                                               ggTTCTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGc  <  2:2738201/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACACCGGTTATGACGATCTCAACAAAAAAACCGCTGGCTTGCAGCCGTCGGATTTGATCATCGTCGCCGCGCGTCCGTCGATGGGTAAAACAACATTTGCGATGAACCTCGTCGAAAACGCGGCGATGTTGCAGGATAAACCGGTACTTATCTTCTCGCTGGAGATGCCATCAGAACAGATCATGATGCGTTCTCTGGCGTCGCTGTCGCGCGTTGACCAGACTAAAATCCGTACCGGGCAGCTCGATGACGAAGACTGGGCGCGCATTTCCGGCACCATGGGTATTTTGC  >  CP009273/4254863‑4255154

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: