Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,890,960 T→A 89.9% F32Y (TTC→TAC)  cbrB → PRK09823 family inner membrane protein, creBC regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,890,9600TA89.9% 883.3 / 74.6 298F32Y (TTC→TAC) cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base T (23/7);  new base A (121/146);  total (144/154)
Fisher's exact test for biased strand distribution p-value = 1.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

Only 100 of 299 total aligned reads displayed.

CCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATC  >  CP009273/3890812‑3891103
                                                                                                                                                    |                                                                                                                                               
  ttCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACg                                                                                                                                               >  2:3276457/1‑149 (MQ=255)
     aCTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCaa                                                                                                                                            >  1:2737560/1‑149 (MQ=255)
         aCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGa                                                                                                                                        >  1:4718203/1‑149 (MQ=255)
         aCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGa                                                                                                                                        >  2:2201068/1‑149 (MQ=255)
           cccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAAc                                                                                                                                      >  1:4018935/1‑149 (MQ=255)
           cccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAAc                                                                                                                                      <  1:2854974/149‑1 (MQ=255)
           cccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAAc                                                                                                                                      >  2:1181326/1‑149 (MQ=255)
            ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAAcc                                                                                                                                     <  2:48329/149‑1 (MQ=255)
              gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGc                                                                                                                                   <  1:3180388/149‑1 (MQ=255)
              gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGc                                                                                                                                   <  1:2328147/149‑1 (MQ=255)
               cTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTGTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCt                                                                                                                                  >  2:2748641/1‑149 (MQ=255)
                  aaaTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGt                                                                                                                               <  1:1054318/149‑1 (MQ=255)
                         cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTgg                                                                                                                        <  1:2366526/149‑1 (MQ=255)
                          gcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTggg                                                                                                                       <  2:15129/149‑1 (MQ=255)
                          gcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTggg                                                                                                                       <  1:771218/149‑1 (MQ=255)
                             cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGa                                                                                                                    >  2:3225799/1‑149 (MQ=255)
                             cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCGGGTTCTTTGGGTGa                                                                                                                    >  2:2670865/1‑149 (MQ=255)
                               aaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGata                                                                                                                  <  2:854569/149‑1 (MQ=255)
                                   taagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATAc                                                                                                              <  1:162531/149‑1 (MQ=255)
                                      gtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                           <  2:1375285/149‑1 (MQ=255)
                                      gtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                           <  1:4723791/149‑1 (MQ=255)
                                       taaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATcc                                                                                                          >  1:4559892/1‑149 (MQ=255)
                                          aGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAt                                                                                                       >  2:1044966/1‑149 (MQ=255)
                                          aGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAt                                                                                                       <  2:1825886/149‑1 (MQ=255)
                                          aGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAt                                                                                                       <  2:843708/149‑1 (MQ=255)
                                            ggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAttt                                                                                                     <  1:3304642/149‑1 (MQ=255)
                                               gTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGAAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTtc                                                                                                  >  1:728440/1‑149 (MQ=255)
                                                  tGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTctct                                                                                               <  1:957376/149‑1 (MQ=255)
                                                   gATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCtt                                                                                              <  2:3358327/149‑1 (MQ=255)
                                                    atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                             <  2:1916121/149‑1 (MQ=255)
                                                    atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                             <  1:568925/149‑1 (MQ=255)
                                                     tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAt                                                                                            <  1:4590172/149‑1 (MQ=255)
                                                     tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAt                                                                                            >  2:588274/1‑149 (MQ=255)
                                                     tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAt                                                                                            <  1:2219266/149‑1 (MQ=255)
                                                      atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAtt                                                                                           <  1:4227042/149‑1 (MQ=255)
                                                      atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACACCCAATTTTTCTCTTAtt                                                                                           <  2:3673530/149‑1 (MQ=255)
                                                       tGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTg                                                                                          <  2:2503258/149‑1 (MQ=255)
                                                        gTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGt                                                                                         >  2:918623/1‑149 (MQ=255)
                                                           tGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGa                                                                                      <  1:2482000/149‑1 (MQ=255)
                                                            gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGAt                                                                                     <  1:133302/149‑1 (MQ=255)
                                                              ttCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATaa                                                                                   <  1:2684783/149‑1 (MQ=255)
                                                                 aCGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTa                                                                                <  2:230265/149‑1 (MQ=255)
                                                                   gTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                              <  2:727051/149‑1 (MQ=255)
                                                                   gTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                              >  2:2792804/1‑149 (MQ=255)
                                                                   gTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                              >  2:400290/1‑149 (MQ=255)
                                                                    tCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                             <  2:3476244/149‑1 (MQ=255)
                                                                      gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGa                                                                           >  1:2614902/1‑149 (MQ=255)
                                                                       ggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                          <  2:4288967/149‑1 (MQ=255)
                                                                       ggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                          >  2:4775906/1‑149 (MQ=255)
                                                                       ggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                          >  2:2978075/1‑149 (MQ=255)
                                                                        gTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGc                                                                         <  1:1814824/149‑1 (MQ=255)
                                                                         tAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCt                                                                        >  2:4718993/1‑149 (MQ=255)
                                                                         tAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCt                                                                        >  1:1322624/1‑149 (MQ=255)
                                                                          aaTACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTa                                                                       >  2:4677848/1‑149 (MQ=255)
                                                                             aCATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTc                                                                    >  2:4436015/1‑149 (MQ=255)
                                                                              cATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTcc                                                                   >  1:3630726/1‑149 (MQ=255)
                                                                                   cGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                              <  1:4635364/149‑1 (MQ=255)
                                                                                   cGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgt                                                              <  2:2399774/149‑1 (MQ=255)
                                                                                    ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                             <  1:450418/149‑1 (MQ=255)
                                                                                     gACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttg                                                            >  1:640052/1‑149 (MQ=255)
                                                                                       cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                          <  1:3611385/149‑1 (MQ=255)
                                                                                       cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                          <  2:2179045/149‑1 (MQ=255)
                                                                                       cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCATGCgttgtt                                                          >  2:2021572/1‑149 (MQ=255)
                                                                                        tttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTTTCTATAACTACTGGAGCTATTCATGAGTTGTTa                                                         >  2:3726104/1‑149 (MQ=255)
                                                                                         ttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTaa                                                        >  1:2740540/1‑149 (MQ=255)
                                                                                           aTTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAcc                                                      <  2:4225576/149‑1 (MQ=255)
                                                                                           aTTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCGACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAcc                                                      <  1:4046659/149‑1 (MQ=255)
                                                                                                gCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtgt                                                 >  1:2451066/1‑149 (MQ=255)
                                                                                                gCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtgt                                                 >  2:2848031/1‑149 (MQ=255)
                                                                                                     tttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATgg                                            <  1:33009/149‑1 (MQ=255)
                                                                                                      ttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGt                                           >  1:1118210/1‑149 (MQ=255)
                                                                                                        aGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTg                                         >  1:380453/1‑149 (MQ=255)
                                                                                                         ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                        >  2:3145362/1‑149 (MQ=255)
                                                                                                         ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGGTGTTAACCGGTGTAATGGTTGc                                        >  1:4404178/1‑149 (MQ=255)
                                                                                                          gACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                       <  1:1789316/149‑1 (MQ=255)
                                                                                                           aCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                      <  1:1960499/149‑1 (MQ=255)
                                                                                                           aCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                      <  1:2184893/149‑1 (MQ=255)
                                                                                                            ccGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTg                                     <  1:3479171/149‑1 (MQ=255)
                                                                                                             cGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGt                                    >  1:1446260/1‑149 (MQ=255)
                                                                                                               ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCt                                  >  2:1703770/1‑149 (MQ=255)
                                                                                                                 aaTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGc                                <  1:797578/149‑1 (MQ=255)
                                                                                                                  aTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                               >  1:3558438/1‑149 (MQ=255)
                                                                                                                   ttGGTGTTCTGTTTCTTGTCCTCTACATATTCTACGCAAAATAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcct                              >  2:2350732/1‑148 (MQ=255)
                                                                                                                      gtgtTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGACGGaa                           >  2:3906861/1‑149 (MQ=255)
                                                                                                                       tgtTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaa                          <  2:2882491/149‑1 (MQ=255)
                                                                                                                       tgtTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaa                          <  2:1325655/149‑1 (MQ=255)
                                                                                                                        gtTCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaaa                         >  1:4121590/1‑149 (MQ=255)
                                                                                                                         ttCTGTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAg                        >  1:1196990/1‑149 (MQ=255)
                                                                                                                             gTTTCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCg                    <  2:2128662/149‑1 (MQ=255)
                                                                                                                              tttCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                   >  2:2602450/1‑149 (MQ=255)
                                                                                                                              tttCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                   >  1:3057663/1‑149 (MQ=255)
                                                                                                                               ttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCggg                  >  1:526602/1‑149 (MQ=255)
                                                                                                                               ttCTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCggg                  <  1:102096/149‑1 (MQ=255)
                                                                                                                                 cTTGTCCTCTACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTc                <  2:1610660/149‑1 (MQ=255)
                                                                                                                                       ctctACATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa          >  1:1423952/1‑149 (MQ=255)
                                                                                                                                            cATATTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGtt     <  1:2978075/149‑1 (MQ=255)
                                                                                                                                             atatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa    >  1:1088385/1‑149 (MQ=255)
                                                                                                                                              tatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTAt   >  1:2787915/1‑149 (MQ=255)
                                                                                                                                              tatTCTACGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTAt   >  1:1234115/1‑149 (MQ=255)
                                                                                                                                               atTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATc  >  2:3050033/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATC  >  CP009273/3890812‑3891103

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: