Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,024,141 G→T 7.9% R13S (CGC→AGC)  hspQ ← heat shock protein involved in degradation of mutant DnaA; hemimethylated oriC DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,024,1410GT 7.9% 638.2 / 19.0 175R13S (CGC→AGC) hspQheat shock protein involved in degradation of mutant DnaA; hemimethylated oriC DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (73/88);  new base T (5/9);  total (78/97)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 176 total aligned reads displayed.

GGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCTTCCGGATAGGT  >  CP009273/1023994‑1024285
                                                                                                                                                   |                                                                                                                                                
   cGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGAc                                                                                                                                               <  1:2009820/148‑1 (MQ=255)
     ttATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGCATGGCGGACct                                                                                                                                             >  2:149283/1‑148 (MQ=255)
      tATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACctg                                                                                                                                            >  2:1332833/1‑148 (MQ=255)
       aTCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCTGACctgc                                                                                                                                           >  1:192597/1‑148 (MQ=255)
       aTCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACctgc                                                                                                                                           >  1:1701577/1‑148 (MQ=255)
        tCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACctgct                                                                                                                                          >  1:1035400/1‑148 (MQ=255)
        tCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACctgct                                                                                                                                          >  2:2000867/1‑148 (MQ=255)
           tcctccATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGc                                                                                                                                       >  2:1394511/1‑148 (MQ=255)
            cctccATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGACTGGCGGACCTGCTGGcc                                                                                                                                      >  1:2575133/1‑148 (MQ=255)
            cctccATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACAGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGcc                                                                                                                                      >  1:3221844/1‑148 (MQ=255)
                cATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATa                                                                                                                                  <  1:410741/148‑1 (MQ=255)
                  tCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTCACAGGGAATGGCGGACCTGCTGGCCGATAcc                                                                                                                                >  1:13019/1‑148 (MQ=255)
                      tACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAAt                                                                                                                            <  2:1973236/148‑1 (MQ=255)
                         gTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTg                                                                                                                         <  2:1265434/148‑1 (MQ=255)
                           gATACCACTGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCt                                                                                                                       <  2:120725/148‑1 (MQ=255)
                           gATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCt                                                                                                                       <  2:805073/148‑1 (MQ=255)
                            aTACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTg                                                                                                                      <  2:1107283/148‑1 (MQ=255)
                            aTACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTg                                                                                                                      >  1:568395/1‑148 (MQ=255)
                             tACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTgg                                                                                                                     >  2:860871/1‑148 (MQ=255)
                               ccACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCa                                                                                                                   <  2:2896904/148‑1 (MQ=255)
                                cACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGGTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCTGGCaa                                                                                                                  <  2:1948460/148‑1 (MQ=255)
                                 aCGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAAt                                                                                                                 >  1:1391437/1‑148 (MQ=255)
                                 aCGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAAt                                                                                                                 >  1:1390584/1‑148 (MQ=255)
                                 aCGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAAt                                                                                                                 <  2:2138601/148‑1 (MQ=255)
                                    gAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGAGATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCAt                                                                                                              <  2:2127220/148‑1 (MQ=255)
                                       cggcgCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGt                                                                                                           <  2:105030/148‑1 (MQ=255)
                                       cggcgCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGt                                                                                                           <  2:1077018/148‑1 (MQ=255)
                                          cgcgAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCAc                                                                                                        <  2:259658/148‑1 (MQ=255)
                                          cgcgAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCGGGCCGATACCGAATTTGCTGGCAATCATAGTCAc                                                                                                        <  2:970316/148‑1 (MQ=255)
                                              aaGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTcc                                                                                                    <  2:1790176/148‑1 (MQ=255)
                                                 cTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCgg                                                                                                 >  1:3036985/1‑148 (MQ=255)
                                                  tCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCggg                                                                                                <  1:2774496/148‑1 (MQ=255)
                                                     tCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGaaa                                                                                             <  1:2618702/148‑1 (MQ=255)
                                                      cATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAAt                                                                                            <  1:1617825/148‑1 (MQ=255)
                                                      cATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAAt                                                                                            <  1:1617819/148‑1 (MQ=255)
                                                         tCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCtt                                                                                         <  2:2258769/148‑1 (MQ=255)
                                                         tCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCtt                                                                                         >  2:2548076/1‑148 (MQ=255)
                                                         tCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCtt                                                                                         <  2:1782014/148‑1 (MQ=255)
                                                            ccGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTAc                                                                                      <  2:808000/148‑1 (MQ=255)
                                                              gCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCt                                                                                    <  2:1628917/148‑1 (MQ=255)
                                                              gCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCt                                                                                    <  2:3266294/148‑1 (MQ=255)
                                                                  gCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGa                                                                                >  2:775439/1‑148 (MQ=255)
                                                                   cTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGat                                                                               <  1:1942060/148‑1 (MQ=255)
                                                                   cTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGat                                                                               >  1:386746/1‑148 (MQ=255)
                                                                       tcaGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATaa                                                                           <  2:155268/148‑1 (MQ=255)
                                                                       tcaGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATaa                                                                           >  2:1100582/1‑148 (MQ=255)
                                                                          ggCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAggg                                                                        <  1:935230/148‑1 (MQ=255)
                                                                              aCGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTaaa                                                                    <  2:1629678/148‑1 (MQ=255)
                                                                               cGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTaaaa                                                                   >  1:931397/1‑148 (MQ=255)
                                                                               cGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTaaaa                                                                   <  1:2961890/148‑1 (MQ=255)
                                                                                 gTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAAtt                                                                 <  2:1544084/148‑1 (MQ=255)
                                                                                  ttCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTc                                                                <  2:2949908/148‑1 (MQ=255)
                                                                                       agagCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACAt                                                           >  1:2194834/1‑148 (MQ=255)
                                                                                       agagCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACAt                                                           <  2:1081913/148‑1 (MQ=255)
                                                                                         agCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACAttt                                                         >  2:2758436/1‑148 (MQ=255)
                                                                                         agCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACAttt                                                         <  1:791127/148‑1 (MQ=255)
                                                                                           cGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGAAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCa                                                       >  1:931070/1‑148 (MQ=255)
                                                                                           cGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGAAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCa                                                       >  1:2585121/1‑148 (MQ=255)
                                                                                              aTAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGt                                                    <  2:945209/148‑1 (MQ=255)
                                                                                               tAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGtt                                                   <  1:1180740/148‑1 (MQ=255)
                                                                                                aaaCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTa                                                  >  2:2131355/1‑148 (MQ=255)
                                                                                                aaaCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTa                                                  >  1:1527399/1‑148 (MQ=255)
                                                                                                 aaCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTAc                                                 >  1:399042/1‑148 (MQ=255)
                                                                                                 aaCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTAc                                                 >  1:2045954/1‑148 (MQ=255)
                                                                                                  aCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACa                                                >  1:2905529/1‑148 (MQ=255)
                                                                                                  aCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACa                                                <  2:5556/148‑1 (MQ=255)
                                                                                                  aCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACa                                                <  2:1701577/148‑1 (MQ=255)
                                                                                                  aCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACa                                                >  1:2905508/1‑148 (MQ=255)
                                                                                                  aCCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCCGTCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACa                                                <  2:1701765/148‑1 (MQ=255)
                                                                                                   ccGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACAt                                               >  2:1027360/1‑148 (MQ=255)
                                                                                                    cGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATg                                              <  2:1885655/148‑1 (MQ=255)
                                                                                                    cGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATg                                              <  1:852800/148‑1 (MQ=255)
                                                                                                           atatCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGAc                                       >  1:3085719/1‑148 (MQ=255)
                                                                                                           atatCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGAc                                       <  1:3230283/148‑1 (MQ=255)
                                                                                                            tatCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACa                                      <  1:2617449/148‑1 (MQ=255)
                                                                                                               ccACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGc                                   >  2:1820394/1‑148 (MQ=255)
                                                                                                                 aCGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCaa                                 >  1:1004107/1‑148 (MQ=255)
                                                                                                                  cGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCaaa                                <  1:1850377/148‑1 (MQ=255)
                                                                                                                     ccACTCCGAGATAACCTAACAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAAccc                             <  2:501792/148‑1 (MQ=255)
                                                                                                                     ccACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAAccc                             >  1:663021/1‑148 (MQ=255)
                                                                                                                     ccACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAAccc                             >  2:432291/1‑148 (MQ=255)
                                                                                                                      cACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCt                            >  2:2195437/1‑148 (MQ=255)
                                                                                                                      cACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCt                            <  2:687936/148‑1 (MQ=255)
                                                                                                                       aCTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCtt                           <  1:2390972/148‑1 (MQ=255)
                                                                                                                        cTCCGAGATAACCGAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCttt                          <  2:3260830/148‑1 (MQ=255)
                                                                                                                           cGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAg                       >  2:2465072/1‑148 (MQ=255)
                                                                                                                            gagaTAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGa                      >  1:393252/1‑148 (MQ=255)
                                                                                                                               aTAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGatat                   <  1:2559995/148‑1 (MQ=255)
                                                                                                                               aTAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGatat                   <  1:1731166/148‑1 (MQ=255)
                                                                                                                                  aCCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGc                >  2:2746513/1‑148 (MQ=255)
                                                                                                                                  aCCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGc                >  1:1801170/1‑148 (MQ=255)
                                                                                                                                  aCCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGc                >  1:2402298/1‑148 (MQ=255)
                                                                                                                                    cTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGccc              <  1:2851252/148‑1 (MQ=255)
                                                                                                                                    cTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGccc              <  1:2004345/148‑1 (MQ=255)
                                                                                                                                      aaCAGGGAATGGCTGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCtt            <  2:399042/148‑1 (MQ=255)
                                                                                                                                        cAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCTTcc          >  1:522688/1‑148 (MQ=255)
                                                                                                                                        cAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCTTcc          >  1:522692/1‑148 (MQ=255)
                                                                                                                                           ggAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCTTCCGGa       <  1:149283/148‑1 (MQ=255)
                                                                                                                                            gAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCTTCCGGAt      <  2:113753/148‑1 (MQ=255)
                                                                                                                                             aaTGGCGGACCTACTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCTTCCGGATa     >  1:3102522/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTGGCCGATACCGAATTTGCTGGCAATCATAGTCACCTCCCGGGAAATCTTTACCTTTGATATAAGGGTAAAATTCCACATTTCAAGTTACATGACGCGACAGGCAAACCCTTTCAGATATAGCCCTTCCGGATAGGT  >  CP009273/1023994‑1024285

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: