Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,673,380 T→A 9.9% S233T (TCA→ACA)  ydgH → DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,673,3800TA 9.9% 337.0 / 3.4 139S233T (TCA→ACA) ydgHDUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (55/67);  new base A (3/12);  total (58/81)
Fisher's exact test for biased strand distribution p-value = 9.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

AGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTG  >  CP009273/1673233‑1673527
                                                                                                                                                   |                                                                                                                                                   
agCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGtcat                                                                                                                                                     >  1:1049932/1‑148 (MQ=255)
agCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGtcat                                                                                                                                                     >  2:2269717/1‑148 (MQ=255)
agCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGtcat                                                                                                                                                     >  2:1642844/1‑148 (MQ=255)
agCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAAACCAGtcat                                                                                                                                                     >  1:1049944/1‑148 (MQ=255)
    cGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCaaa                                                                                                                                                 >  2:895496/1‑148 (MQ=255)
    cGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCaaa                                                                                                                                                 <  2:2690129/148‑1 (MQ=255)
      gATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCCAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAg                                                                                                                                               >  2:2919487/1‑148 (MQ=255)
      gATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAg                                                                                                                                               <  1:2885917/148‑1 (MQ=255)
       aTGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAgg                                                                                                                                              <  1:749854/148‑1 (MQ=255)
       aTGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAgg                                                                                                                                              <  1:2266383/148‑1 (MQ=255)
       aTGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAgg                                                                                                                                              <  1:2266369/148‑1 (MQ=255)
        tgtgATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAggc                                                                                                                                             <  2:3221683/148‑1 (MQ=255)
         gtgATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAggcg                                                                                                                                            <  1:2619831/148‑1 (MQ=255)
          tgATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCGTCTTTGTAACCCACTCATCAAAAtgcgg                                                                                                                                           >  1:620507/1‑148 (MQ=255)
             tCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAtgcgggcg                                                                                                                                        >  1:2835603/1‑148 (MQ=255)
             tCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAggcgggcg                                                                                                                                        <  2:514712/148‑1 (MQ=255)
              cccGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTGTTTGTAACCCCGTCATCAAAAGGCGGGCGt                                                                                                                                       >  1:1500485/1‑148 (MQ=255)
              cccGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGt                                                                                                                                       <  2:2293511/148‑1 (MQ=255)
              cccGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGt                                                                                                                                       <  2:982558/148‑1 (MQ=255)
                 ggCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTac                                                                                                                                    <  2:2382351/148‑1 (MQ=255)
                 ggCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTac                                                                                                                                    <  2:2382786/148‑1 (MQ=255)
                  gCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAAACCACTCATCAAAAGGCGGGCGTTaaa                                                                                                                                   >  1:678626/1‑146 (MQ=255)
                   cAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTacac                                                                                                                                  >  2:2035137/1‑148 (MQ=255)
                     gATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCg                                                                                                                                <  2:1327402/148‑1 (MQ=255)
                      aTTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACGGt                                                                                                                               <  1:1374735/148‑1 (MQ=255)
                      aTTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACGGt                                                                                                                               <  1:2140449/148‑1 (MQ=255)
                      aTTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGt                                                                                                                               <  2:104218/148‑1 (MQ=255)
                        tCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCa                                                                                                                             <  1:956044/148‑1 (MQ=255)
                          cGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACg                                                                                                                           <  2:2936740/148‑1 (MQ=255)
                             aGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTc                                                                                                                        <  2:1538884/148‑1 (MQ=255)
                              gCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTcc                                                                                                                       <  1:423205/148‑1 (MQ=255)
                              gCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTAAACCGTCACGCTcc                                                                                                                       >  1:1749032/1‑148 (MQ=255)
                                 ggACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCgg                                                                                                                    <  2:2574220/148‑1 (MQ=255)
                                 ggACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCgg                                                                                                                    <  2:2573780/148‑1 (MQ=255)
                                 ggACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGGTACACCGTCACGCTCCCgg                                                                                                                    >  2:1999392/1‑148 (MQ=255)
                                  gACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGa                                                                                                                   >  2:289969/1‑148 (MQ=255)
                                  gACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGa                                                                                                                   >  1:111925/1‑148 (MQ=255)
                                    cGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATg                                                                                                                 <  1:2043507/148‑1 (MQ=255)
                                       gCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCa                                                                                                              <  2:3282696/148‑1 (MQ=255)
                                           ctctGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTGTTTGAAACCCAGTCATCAAAAGGCGGGCGGTACACCGTCACGCTCCCGGATGGCACTaa                                                                                                          >  1:2359636/1‑148 (MQ=255)
                                           ctctGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGGTACACCGTCACGCTCCCGGATGGcccgaa                                                                                                          >  1:2361312/1‑143 (MQ=255)
                                           ctctGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATAAAAAGGCGGGCGTTACACCGTCACGCTCCCCGATGGCACTaa                                                                                                          >  1:734436/1‑148 (MQ=255)
                                                     cGGTGTCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGAGACTACCGCATCACCTAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCgaagaa                                                                                                <  1:3114035/148‑1 (MQ=255)
                                                         ggCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGa                                                                                            <  1:2522227/148‑1 (MQ=255)
                                                            gAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACa                                                                                         <  1:97824/148‑1 (MQ=255)
                                                             aaGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATAAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTTAACaa                                                                                        >  2:1915450/1‑148 (MQ=255)
                                                              aGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACcaa                                                                                       >  2:1357039/1‑148 (MQ=255)
                                                               gCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAg                                                                                      <  2:1235866/148‑1 (MQ=255)
                                                               gCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAg                                                                                      <  2:1170365/148‑1 (MQ=255)
                                                                     aagaaAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTATTTTAAACCCAGTCAACAAAAGGCGGGCGTTAAAACCTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAACGAccc                                                                                >  1:659022/1‑147 (MQ=255)
                                                                      agaaAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAATCGACCGc                                                                               >  2:190312/1‑148 (MQ=255)
                                                                          aGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCGTCTTTTAAACCCAGTCATCAAAAGGCGGGCGGTACACCCTCACGCTCCCCGATGGCACTAAAGTCGAAGAACAGAACAAAACCACCGCAGCg                                                                           >  1:81496/1‑148 (MQ=255)
                                                                            ttGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACGAAAGTCGAAGAACTGAACAAAGCGACCGCAGCgct                                                                         >  2:632443/1‑146 (MQ=255)
                                                                             tGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTGGGTCGCTTCTTTGAACCCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCgatg                                                                        <  2:2979323/148‑1 (MQ=255)
                                                                              tataTCCCGGTTGTTGCGATTACCGCATCACAAATTTCTCAAGTCGGTCGATTATTTTAAACCCATTCAAAAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGAACTAAATTCGAAGAACTGAAAAAAGCGACCGCAGCgatga                                                                       <  2:558675/145‑1 (MQ=255)
                                                                              gagaTCCGGGGTTTTGCGACTACCGCATCACCAAGTTCTGAAGTCGTTCGCTTCTTTGAAACCCAGTCATAAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGAACTAAAGTCGAAGAACTGAACAACGCGACCGCAGCgatga                                                                       <  2:719699/148‑1 (MQ=255)
                                                                               agaTCCGGGGTGTTGCGACTCCCGCATCACCATGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAACGGGGGGCGTTACACCGTCACGCTCCGGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCgatgat                                                                      <  2:1590136/148‑1 (MQ=255)
                                                                               agaTCCGGGGTGTTGCGACTACCGAATCACCATGTTCTTATGTCGGTCGCTTATTTGAAACCCAGTCATAAAAAGGGGGGCGTTACACCGTCACGCTCCGGTATGGCACTAAAGTCGAATAACTCAACAAAGCGACCGCAGCgatgat                                                                      <  2:2835603/148‑1 (MQ=255)
                                                                                 aTCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAATCGACCGCAGCGATGATgg                                                                    >  2:2608616/1‑148 (MQ=255)
                                                                                    ccGGTTGTTGCGACTACCGCATCACAAATTTATGAAGTCGGTCGATAATTTTAAACCCAGACAAAAAAAGGCGGGCGTTACACCGACACGCACCCGGATGGCACTAAAGTCCAAGAACTGAAAAAAGCGACCGCAGCGATGATGGTcc                                                                 <  2:723721/148‑1 (MQ=255)
                                                                                    ccGGGTGTTGCCAATCCCGCATCACCAAGTTCTGACGTCGGTCGCTTCTTTGAAACCCCGTCATCAAACGGCGGGGGTTACACCGTCCCGCTCCCGGATGGCACTAACGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTcc                                                                 <  2:2657183/148‑1 (MQ=255)
                                                                                        gtggtgCGACTACCCAATCACCAAGTGCTGAAGGCGGTCGCTTCTTTGAAACCCAGCCATCAAACGGCGGGCGTAACACCGTCACGCACCCGGATGGCACTAAAGTCGAAGAACTGACCAAAGCGACCGCAGCGATGATGGTCCCGtt                                                             <  2:2005924/148‑1 (MQ=255)
                                                                                                aCTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAATAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCa                                                     >  2:1839019/1‑148 (MQ=255)
                                                                                                aCTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCa                                                     >  2:1839042/1‑148 (MQ=255)
                                                                                                 ctcccGCATCACCTTTTTGTGTTGTGGGTCGTTTTTTTTAAACCACTTAATCAAACGGGGGGGTTTACTCGTTCACGCTCCGTTTTGCCACTAATGTCGATGAACTTAAAAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCAt                                                    <  2:303344/145‑1 (MQ=255)
                                                                                                   aCCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCa                                                  >  2:821451/1‑148 (MQ=255)
                                                                                                   aCCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATAAAAAGGCGGGCGTTAAACCGTCACGCTCCCGGATGGCACAAAAGTCGAAGAACAGAACAAAGCGACCGCAGCGATGATTGTCCCGGTCGACAGCAACa                                                  >  1:507516/1‑148 (MQ=255)
                                                                                                     cGCATCACAAAGTACTGACGTGGGTCGCTTCTTTGAAACCCAGTCAAAAAACGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGACGAACTGAACAACGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCaaa                                                <  1:817846/148‑1 (MQ=255)
                                                                                                      gCAACACCAAGTTCTGAAGTCGTTCGCTTATTTCAAACCCAGACAAAAAACGGCGGGCGTTACACCGTCACGCTCCCGGAGGCCACAAACGCCCAACAACTCAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAAt                                               <  1:3271182/148‑1 (MQ=255)
                                                                                                       cATCACCAAGTTCTGAAGTCGGTCGCGTGTTTTTAACCCAGTCATCAAAAGGCGGGCGGTACAACCTCCCCCGCCCCGCTGGCACTAAAGTCGAAGAACAGAACAAAGCGACCGCAGCCATGATGGTCCCCTTCGACACCCTCAAAtt                                              >  1:957713/1‑148 (MQ=255)
                                                                                                         tCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGGTCGACAGCATCAAATtct                                            >  2:3203719/1‑148 (MQ=255)
                                                                                                         tCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACGAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATtct                                            >  2:3204166/1‑148 (MQ=255)
                                                                                                         tCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACGAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATtct                                            >  2:1929578/1‑148 (MQ=255)
                                                                                                          cacaatgTTCTGATGTCGGTCGGTTATTTGAAACCCAGTCATAAAAAGGCGGGCGTTACACCGACACGCTCCCGGATGGAACTAAAGTCGAACAACTGAAAAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATtctc                                           <  1:494216/142‑1 (MQ=255)
                                                                                                          cACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATtcta                                           >  2:2479176/1‑147 (MQ=255)
                                                                                                            ccAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTg                                         >  2:1187139/1‑148 (MQ=255)
                                                                                                                  tCTGAAGTGGGTCGGTTCTTTTAAACCCAGTAATAAAAAGGCGTGCTTTACACCGTCACGCTCCCGTATGGCACTAACGGCGACGTATTTAAAAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACt                                   <  2:692288/148‑1 (MQ=255)
                                                                                                                  tCTGAAGTCGGTCGCTTCTTTTAAACCCAGTCATAAAAAGGCGGGCTTTACACCGGCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAACGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACt                                   <  2:2359636/148‑1 (MQ=255)
                                                                                                                   cTGACGTCGGTCGCTTCTTTGAACCCCCGTCATCAACAGGCGGGCGTAACCCCGTCACGCGCCCGGATGGCACTAAAGTCGAAGAACTGAAAAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTa                                  <  2:654392/148‑1 (MQ=255)
                                                                                                                      aaGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACgg                               >  2:3006196/1‑148 (MQ=255)
                                                                                                                        gtcggtcgCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCa                             <  2:206020/148‑1 (MQ=255)
                                                                                                                           ggtcgCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACa                          <  2:2182522/148‑1 (MQ=255)
                                                                                                                             tcgatacttTGAAACCCAGTCATAAAATGGCGGGCGTTAACCCGACACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAACGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAAGTCTCTGGCAACTACGGCAACATg                        <  2:3158562/142‑1 (MQ=255)
                                                                                                                              cgtttcttTCAAACCCAGGCAACAAAGGGCGGGGGTTACCCCGCCCCGCTCCCGGTTGGCACAAACGGCCACCAACTTAACAACGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGa                       <  1:2641534/145‑1 (MQ=255)
                                                                                                                                cttcttTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGAcc                     >  1:341113/1‑148 (MQ=255)
                                                                                                                                  tcttTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGa                   >  1:3040086/1‑148 (MQ=255)
                                                                                                                                  tattttAAACCCACACAAAAAATGGGGGGGGTAACCCCGACACGCTCCGGTTTGCCAATAATGACCAATAAATCAAAAACGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGa                   <  1:360191/146‑1 (MQ=255)
                                                                                                                                    tttGAACCCCAGTCACCAAACGGGGGGGGTTACCCCGTCCCGCACCGGGATGGCACTAACGTCGAAGAACTGACCAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAg                 <  1:3133264/148‑1 (MQ=255)
                                                                                                                                         aaCCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTcc            >  2:81574/1‑148 (MQ=255)
                                                                                                                                         aaCCCAGTCAACAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTcc            <  2:116744/148‑1 (MQ=255)
                                                                                                                                          aCCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCt           >  1:2458270/1‑148 (MQ=255)
                                                                                                                                          aCCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCt           >  2:2539087/1‑148 (MQ=255)
                                                                                                                                              aGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCa       >  2:1903764/1‑148 (MQ=255)
                                                                                                                                               gTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAg      >  1:3301078/1‑148 (MQ=255)
                                                                                                                                               gTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAg      >  1:1318858/1‑148 (MQ=255)
                                                                                                                                                tcatcaAATGGCGGGCGTTCCACCGTCACGCTCCCGGATGGCACAAAAGTCGAACAACTGAACCAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAgg     <  1:1753315/148‑1 (MQ=255)
                                                                                                                                                tcatcaAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAgg     >  2:1435588/1‑148 (MQ=255)
                                                                                                                                                tcatcaAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAgg     >  1:1176629/1‑148 (MQ=255)
                                                                                                                                                tcataaAACGGGGGGCTTTACACCGACACGCACCGGGTTGGCAATAAAGTCCAACAACTCAAAAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAgg     <  1:1572090/148‑1 (MQ=255)
                                                                                                                                                   tcaAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTg  >  2:1521328/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCCCGGATGTGATCCCGGCAGATTCCGAAGCAGGACGTGCAGCTCTGGCTGCCGGTGGCGAAGCCGCGAAGAAAGTTGAGATCCCGGGTGTTGCGACTACCGCATCACCAAGTTCTGAAGTCGGTCGCTTCTTTGAAACCCAGTCATCAAAAGGCGGGCGTTACACCGTCACGCTCCCGGATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTG  >  CP009273/1673233‑1673527

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: