Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,167,995 G→T 29.5% D109Y (GAT→TAT)  nusG → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,167,9950GT29.5% 646.8 / 235.0 302D109Y (GAT→TAT) nusGtranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (109/104);  new base T (40/49);  total (149/153)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 304 total aligned reads displayed.

AGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTG  >  CP009273/4167848‑4168144
                                                                                                                                                   |                                                                                                                                                     
  cGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGc                                                                                                                                                   >  2:3844505/1‑150 (MQ=255)
     aCGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAt                                                                                                                                                >  1:3873485/1‑150 (MQ=255)
      cGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGAtt                                                                                                                                               <  2:4616487/150‑1 (MQ=255)
       gTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTa                                                                                                                                              <  1:3205103/150‑1 (MQ=255)
       gTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTa                                                                                                                                              >  2:4279152/1‑150 (MQ=255)
        tAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTAt                                                                                                                                             >  2:929015/1‑150 (MQ=255)
        tAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTAt                                                                                                                                             >  2:3252010/1‑150 (MQ=255)
        tAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTAt                                                                                                                                             <  1:4498454/150‑1 (MQ=255)
                tcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGcc                                                                                                                                     >  1:1604716/1‑150 (MQ=255)
                tcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGcc                                                                                                                                     >  2:1097748/1‑150 (MQ=255)
                tcCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGcc                                                                                                                                     <  1:2254052/150‑1 (MQ=255)
                   cTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTgc                                                                                                                                  >  2:795019/1‑150 (MQ=255)
                    tGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTgca                                                                                                                                 >  2:1223341/1‑150 (MQ=255)
                     ggCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTgcag                                                                                                                                >  2:3813421/1‑150 (MQ=255)
                       cTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTgcagca                                                                                                                              <  2:3730294/150‑1 (MQ=255)
                        tACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGcagcag                                                                                                                             >  2:2628865/1‑150 (MQ=255)
                           gTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGtt                                                                                                                          >  2:4099037/1‑150 (MQ=255)
                           gTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGtt                                                                                                                          <  2:2718939/150‑1 (MQ=255)
                             ccTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTgg                                                                                                                        <  1:3421700/150‑1 (MQ=255)
                                cGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGa                                                                                                                     >  1:2891155/1‑150 (MQ=255)
                                cGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGa                                                                                                                     >  1:907423/1‑150 (MQ=255)
                                 gTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGAt                                                                                                                    >  1:882206/1‑150 (MQ=255)
                                 gTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGAt                                                                                                                    <  2:2529020/150‑1 (MQ=255)
                                   tCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATaa                                                                                                                  >  2:365497/1‑150 (MQ=255)
                                        tGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCcgc                                                                                                             >  1:1011080/1‑150 (MQ=255)
                                          gTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGt                                                                                                           >  2:1524262/1‑150 (MQ=255)
                                          gTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGt                                                                                                           >  2:871413/1‑150 (MQ=255)
                                           tgatgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTc                                                                                                          >  1:4638726/1‑150 (MQ=255)
                                           tgatgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTc                                                                                                          <  2:2097717/150‑1 (MQ=255)
                                             atgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCg                                                                                                        >  1:585094/1‑150 (MQ=255)
                                             atgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCg                                                                                                        >  2:2564298/1‑150 (MQ=255)
                                             atgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCg                                                                                                        <  1:2282539/150‑1 (MQ=255)
                                              tgaACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGa                                                                                                       >  2:1945003/1‑150 (MQ=255)
                                                aaCGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGaaa                                                                                                     >  1:1399179/1‑150 (MQ=255)
                                                  cgacgCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAAc                                                                                                   <  2:2865112/150‑1 (MQ=255)
                                                  cgacgCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAAc                                                                                                   >  1:115625/1‑150 (MQ=255)
                                                     cgcgAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCt                                                                                                <  1:3763434/150‑1 (MQ=255)
                                                       cgAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGt                                                                                              >  2:2270285/1‑150 (MQ=255)
                                                        gAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGtt                                                                                             >  1:1781382/1‑150 (MQ=255)
                                                         aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGttt                                                                                            >  1:2476158/1‑150 (MQ=255)
                                                         aGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGttt                                                                                            <  2:3944167/150‑1 (MQ=255)
                                                          gCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTg                                                                                           <  1:4202380/150‑1 (MQ=255)
                                                          gCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTg                                                                                           >  2:1773166/1‑150 (MQ=255)
                                                             ggCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAAc                                                                                        >  1:493523/1‑150 (MQ=255)
                                                                aCCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCgg                                                                                     >  2:659804/1‑150 (MQ=255)
                                                                aCCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCgg                                                                                     >  2:4435123/1‑150 (MQ=255)
                                                                aCCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCgg                                                                                     >  1:3996079/1‑150 (MQ=255)
                                                                   tGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTg                                                                                  >  2:3396237/1‑150 (MQ=255)
                                                                   tGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTg                                                                                  >  2:4353167/1‑150 (MQ=255)
                                                                    ggTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGa                                                                                 <  2:1750338/150‑1 (MQ=255)
                                                                     gTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGaa                                                                                >  2:129172/1‑150 (MQ=255)
                                                                       gcgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAAt                                                                              <  2:4382466/150‑1 (MQ=255)
                                                                       gcgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAAt                                                                              <  2:2052832/150‑1 (MQ=255)
                                                                        cgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                                                             <  1:2500336/150‑1 (MQ=255)
                                                                        cgcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATg                                                                             <  1:3668379/150‑1 (MQ=255)
                                                                         gcAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATgg                                                                            >  1:3350641/1‑150 (MQ=255)
                                                                            gCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTcc                                                                         <  1:3387650/150‑1 (MQ=255)
                                                                             cGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCg                                                                        >  1:4100893/1‑150 (MQ=255)
                                                                               tACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgtg                                                                      >  1:1546024/1‑150 (MQ=255)
                                                                                aCCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCgtgt                                                                     >  2:2415089/1‑150 (MQ=255)
                                                                                    cgTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAat                                                                 >  1:1339142/1‑150 (MQ=255)
                                                                                    cgTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAat                                                                 <  1:4525881/150‑1 (MQ=255)
                                                                                     gtgtGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAatg                                                                <  1:3954101/150‑1 (MQ=255)
                                                                                         gATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATgg                                                            >  2:221094/1‑150 (MQ=255)
                                                                                                  cATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGc                                                   >  2:107644/1‑150 (MQ=255)
                                                                                                   aTCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCt                                                  <  2:1378224/150‑1 (MQ=255)
                                                                                                    tCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTg                                                 <  2:2220894/150‑1 (MQ=255)
                                                                                                     cggcggTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGa                                                <  2:1935576/150‑1 (MQ=255)
                                                                                                      ggcggTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGAc                                               >  1:187658/1‑150 (MQ=255)
                                                                                                       gcggTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACt                                              >  2:3567388/1‑150 (MQ=255)
                                                                                                         ggTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTc                                            <  1:3844505/150‑1 (MQ=255)
                                                                                                          gTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCa                                           <  2:3092817/150‑1 (MQ=255)
                                                                                                           tACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCaa                                          <  1:2943882/150‑1 (MQ=255)
                                                                                                           tACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCTTGTTAATGATGGTCCGTTCGCTGACTTCaa                                          <  1:4289614/150‑1 (MQ=255)
                                                                                                             cTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACg                                        <  1:93085/150‑1 (MQ=255)
                                                                                                              ttCCGATCGTCCTGCGCCAATCAGCGATAANGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACgg                                       >  1:4281674/1‑150 (MQ=255)
                                                                                                              ttCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACgg                                       >  2:3646428/1‑150 (MQ=255)
                                                                                                               tCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGgt                                      >  1:302122/1‑150 (MQ=255)
                                                                                                                ccGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGgtg                                     >  2:4157740/1‑150 (MQ=255)
                                                                                                                    tCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGtgttgt                                 >  1:1928498/1‑150 (MQ=255)
                                                                                                                     cGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTgttgtt                                >  2:3894246/1‑150 (MQ=255)
                                                                                                                     cGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTgttgtt                                >  2:1702061/1‑150 (MQ=255)
                                                                                                                     cGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGAGGGTCCGTTCGCTGACTTCAACGGTgttgtt                                >  2:1009935/1‑150 (MQ=255)
                                                                                                                      gTCCTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGttgttg                               >  2:1169597/1‑150 (MQ=255)
                                                                                                                      gTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCAGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGATAATGATGGTCAGTTCTCTGACTTCAACGGTGttgttg                               >  1:1204632/1‑150 (MQ=255)
                                                                                                                        ccTGCGCCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaa                             >  2:3429058/1‑150 (MQ=255)
                                                                                                                        ccTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaa                             <  1:2792633/150‑1 (MQ=255)
                                                                                                                           gcgcCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaagaa                          <  2:3346801/150‑1 (MQ=255)
                                                                                                                           gcgcCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaagaa                          <  1:607303/150‑1 (MQ=255)
                                                                                                                           gcgcCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTgaagaa                          <  1:2124904/150‑1 (MQ=255)
                                                                                                                             gcCAATCAGCGATAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGt                        <  2:2374599/150‑1 (MQ=255)
                                                                                                                             gcCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGt                        <  2:182856/150‑1 (MQ=255)
                                                                                                                              ccAATCAGCGATACAGAAGTCGATGAGATTATGAAACGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTg                       <  1:1644181/150‑1 (MQ=255)
                                                                                                                                    aGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTAc                 <  2:894360/150‑1 (MQ=255)
                                                                                                                                        aTAAAGAAGTCTATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACgaga             <  1:2160871/150‑1 (MQ=255)
                                                                                                                                             gaagTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAAtct        <  2:3445924/150‑1 (MQ=255)
                                                                                                                                              aagTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAAtctc       <  1:4504016/150‑1 (MQ=255)
                                                                                                                                                 tCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTc    <  2:592025/150‑1 (MQ=255)
                                                                                                                                                 tCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTc    >  1:3855062/1‑150 (MQ=255)
                                                                                                                                                  cGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCt   <  2:1993655/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCGAACGTAAATTCTTCCCTGGCTACGTCCTCGTTCAGATGGTGATGAACGACGCGAGCTGGCACCTGGTGCGCAGCGTACCGCGTGTGATGGGCTTCATCGGCGGTACTTCCGATCGTCCTGCGCCAATCAGCGATAAAGAAGTCGATGCGATTATGAACCGCCTGCAGCAGGTTGGTGATAAGCCGCGTCCGAAAACGCTGTTTGAACCGGGTGAAATGGTCCGTGTTAATGATGGTCCGTTCGCTGACTTCAACGGTGTTGTTGAAGAAGTGGATTACGAGAAATCTCGTCTG  >  CP009273/4167848‑4168144

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: