Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 608,575 T→A 9.8% L128Q (CTG→CAG)  fes → enterobactin/ferric enterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273608,5750TA 9.8% 652.1 / 52.5 225L128Q (CTG→CAG) fesenterobactin/ferric enterobactin esterase
Reads supporting (aligned to +/- strand):  ref base T (109/94);  new base A (11/11);  total (120/105)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

Only 100 of 227 total aligned reads displayed.

GACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTTGAAAAA  >  CP009273/608429‑608722
                                                                                                                                                  |                                                                                                                                                   
gacgacACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTg                                                                                                                                                    >  2:1246245/1‑148 (MQ=255)
 acgacACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGa                                                                                                                                                   >  1:890876/1‑148 (MQ=255)
    acacAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAAcc                                                                                                                                                <  1:3323302/148‑1 (MQ=255)
        aaCTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTCTTACCCCAGGCGATAGCCGATCCGCTGAACCTACa                                                                                                                                            >  2:2489305/1‑148 (MQ=255)
         aCTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACaa                                                                                                                                           <  1:1027372/148‑1 (MQ=255)
          cTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACaaa                                                                                                                                          <  1:472346/148‑1 (MQ=255)
           tCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAg                                                                                                                                         >  2:1010434/1‑148 (MQ=255)
           tCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAg                                                                                                                                         >  2:2798905/1‑148 (MQ=255)
            cAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGc                                                                                                                                        >  2:2604210/1‑148 (MQ=255)
             aaTGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCt                                                                                                                                       >  2:703201/1‑148 (MQ=255)
              aTGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTg                                                                                                                                      >  2:604424/1‑148 (MQ=255)
               tGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTgg                                                                                                                                     >  2:113348/1‑148 (MQ=255)
               tGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTgg                                                                                                                                     >  2:2107211/1‑148 (MQ=255)
               tGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTgg                                                                                                                                     >  2:314330/1‑148 (MQ=255)
                gCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGa                                                                                                                                    <  1:102626/148‑1 (MQ=255)
                 ccAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGaa                                                                                                                                   <  2:3297267/148‑1 (MQ=255)
                     cTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAggc                                                                                                                               >  1:1308001/1‑148 (MQ=255)
                     cTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAggc                                                                                                                               <  2:836221/148‑1 (MQ=255)
                     cTGGCGCGGCAGCTACTGCTTTAGTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAggc                                                                                                                               <  2:1456418/148‑1 (MQ=255)
                       ggCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAggcgg                                                                                                                             <  2:253431/148‑1 (MQ=255)
                        gcgcGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAggcggg                                                                                                                            <  1:1467506/148‑1 (MQ=255)
                         cgcgGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAggcgggc                                                                                                                           >  1:1711240/1‑148 (MQ=255)
                             gCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAgg                                                                                                                       >  1:774013/1‑148 (MQ=255)
                                gCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCa                                                                                                                    <  2:3376423/148‑1 (MQ=255)
                                gCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTAACATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAATACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCa                                                                                                                    <  2:355599/148‑1 (MQ=255)
                                 cTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCAc                                                                                                                   >  1:231860/1‑148 (MQ=255)
                                  tACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACg                                                                                                                  >  2:3409690/1‑148 (MQ=255)
                                  tACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATGACCCCAGGCGATAGCCGATCCGCTGAACCTACAACGCTGGAAAGGCGGGCGAGGGCACg                                                                                                                  <  2:2675539/148‑1 (MQ=255)
                                   aCTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGc                                                                                                                 >  1:2968083/1‑148 (MQ=255)
                                    cTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCt                                                                                                                >  1:362141/1‑148 (MQ=255)
                                    cTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCt                                                                                                                >  1:1650195/1‑148 (MQ=255)
                                     tGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTg                                                                                                               >  1:1113040/1‑148 (MQ=255)
                                     tGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTg                                                                                                               >  2:2086844/1‑148 (MQ=255)
                                             tCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCa                                                                                                       <  1:3501291/148‑1 (MQ=255)
                                                  ccGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGa                                                                                                  <  2:2507384/148‑1 (MQ=255)
                                                   cGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGaa                                                                                                 <  2:1711240/148‑1 (MQ=255)
                                                    gAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGaaa                                                                                                <  1:1071985/148‑1 (MQ=255)
                                                    gAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGaaa                                                                                                <  2:228259/148‑1 (MQ=255)
                                                      aCGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATg                                                                                              <  1:2369697/148‑1 (MQ=255)
                                                           aTGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCa                                                                                         <  1:1154528/148‑1 (MQ=255)
                                                           aTGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCa                                                                                         >  1:1914865/1‑148 (MQ=255)
                                                            tGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCaa                                                                                        >  1:3205348/1‑148 (MQ=255)
                                                            tGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCaa                                                                                        >  1:2148350/1‑148 (MQ=255)
                                                            tGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCaa                                                                                        >  1:3108165/1‑148 (MQ=255)
                                                            tGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCaa                                                                                        >  1:2842541/1‑148 (MQ=255)
                                                                   ttttCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCctc                                                                                 >  1:3343970/1‑148 (MQ=255)
                                                                   ttttCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCctc                                                                                 <  2:1445750/148‑1 (MQ=255)
                                                                      tCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGc                                                                              >  1:3196508/1‑148 (MQ=255)
                                                                        tGTACCAGCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCGGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCaa                                                                            <  2:1106950/148‑1 (MQ=255)
                                                                         gTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGATGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAAc                                                                           >  2:3220151/1‑148 (MQ=255)
                                                                            ccATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgg                                                                        >  1:3226260/1‑148 (MQ=255)
                                                                            ccATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgg                                                                        >  1:530828/1‑148 (MQ=255)
                                                                             cATCCCCCGATCGCCTCGATTTGCGCGAAGGCTGGCTAAAAATATTACCCCCGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCggg                                                                       <  1:1374688/148‑1 (MQ=255)
                                                                             cATCCCCCGATCGCCTCGAATTGCGCGATGGCTGGCGAAAACTATTACCCCTGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCTAGGGCACGCTGTTTCTGAACTCTAAATGCCGCAAGCGCCTCTGCAACCggg                                                                       <  1:2839448/148‑1 (MQ=255)
                                                                             cATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCggg                                                                       >  2:3210365/1‑148 (MQ=255)
                                                                              aTCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCggga                                                                      <  2:3463416/148‑1 (MQ=255)
                                                                                cccccGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCgggatg                                                                    >  1:2963852/1‑148 (MQ=255)
                                                                                       tCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAACGCTGGAACGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGt                                                             <  1:415418/148‑1 (MQ=255)
                                                                                       tCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGt                                                             >  2:78048/1‑148 (MQ=255)
                                                                                        cGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTc                                                            <  2:3217985/148‑1 (MQ=255)
                                                                                        cGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTc                                                            >  1:1214193/1‑148 (MQ=255)
                                                                                            tCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCa                                                        >  2:1866432/1‑148 (MQ=255)
                                                                                             cGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCaa                                                       >  2:847863/1‑148 (MQ=255)
                                                                                               aaTTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAgc                                                     >  2:1838841/1‑148 (MQ=255)
                                                                                                 ttGCGCTAAGGCTGGCGAAAACTATTCCCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGCGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAgcgc                                                   <  2:2912892/148‑1 (MQ=255)
                                                                                                 ttGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAgcgc                                                   >  1:2191502/1‑148 (MQ=255)
                                                                                                     gcgAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGa                                               <  1:2647646/148‑1 (MQ=255)
                                                                                                        aaGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAAt                                            <  2:1674260/148‑1 (MQ=255)
                                                                                                         aGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGCTGGGATTGTCCGCAAGCGCCAGAAATa                                           >  1:2769421/1‑148 (MQ=255)
                                                                                                         aGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATa                                           <  2:433644/148‑1 (MQ=255)
                                                                                                          ggCTGGCGAAAACTATTCCCCCCGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAACGGGGGGCGCGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAc                                          <  1:2160763/148‑1 (MQ=255)
                                                                                                          ggCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAc                                          <  2:1489381/148‑1 (MQ=255)
                                                                                                           gCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAcc                                         >  1:2921536/1‑148 (MQ=255)
                                                                                                           gCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCAGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATAcc                                         <  1:1743335/148‑1 (MQ=255)
                                                                                                              ggCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGc                                      >  2:1772960/1‑148 (MQ=255)
                                                                                                              ggCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGc                                      <  2:3046375/148‑1 (MQ=255)
                                                                                                                   aaaCTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAg                                 >  1:2749755/1‑148 (MQ=255)
                                                                                                                    aaCTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGa                                <  2:1360065/148‑1 (MQ=255)
                                                                                                                    aaCTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGa                                <  2:2148984/148‑1 (MQ=255)
                                                                                                                      ctatgcccccCGGCTATAGCCGATCCGCTGAACCTACAACGCTGGAAAGGCGGGCGTGGGCACGCTGTTTCTGCCCTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGaaa                              <  1:2096140/142‑1 (MQ=255)
                                                                                                                        aTTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAAtt                            <  2:421986/148‑1 (MQ=255)
                                                                                                                        aTTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAAtt                            >  1:1557985/1‑148 (MQ=255)
                                                                                                                          tACCCCAGGCGATAGCCGATCCGCTGAACCTACAACGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTAt                          <  2:1920148/148‑1 (MQ=255)
                                                                                                                          tACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTAt                          >  2:2997437/1‑148 (MQ=255)
                                                                                                                          tACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTAt                          <  2:1650195/148‑1 (MQ=255)
                                                                                                                           cccccCGGCGATAGCCTATCCGCTGAACCTACAACGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCTAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATc                         <  1:719964/147‑1 (MQ=255)
                                                                                                                                  gctaTAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGaaaa                  <  2:932841/145‑1 (MQ=255)
                                                                                                                                  gCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGaaaa                  >  2:805402/1‑148 (MQ=255)
                                                                                                                                    gATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGt                <  1:546125/148‑1 (MQ=255)
                                                                                                                                     aTAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTg               >  1:1229094/1‑148 (MQ=255)
                                                                                                                                      tAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGa              >  1:3414501/1‑148 (MQ=255)
                                                                                                                                       aGCCGATCCGCTGAACCTACAAAGCTGGAAATGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGaa             >  1:2784850/1‑148 (MQ=255)
                                                                                                                                       aGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGaa             <  2:1488038/148‑1 (MQ=255)
                                                                                                                                        gCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAAc            >  2:1201462/1‑148 (MQ=255)
                                                                                                                                           gATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGt         >  1:492049/1‑148 (MQ=255)
                                                                                                                                           gATCCGCAGAACCTACAACGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCTAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGt         <  2:2558992/148‑1 (MQ=255)
                                                                                                                                            aTCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGtt        >  1:2263114/1‑148 (MQ=255)
                                                                                                                                                 cTGAACCTACACAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTTGaaaa   <  1:2700074/148‑1 (MQ=255)
                                                                                                                                                 cTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTTGaaaa   <  2:286239/148‑1 (MQ=255)
                                                                                                                                                 cTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTTGaaaa   <  1:1053432/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACGACACAACTCAATGCCAACTGGCGCGGCAGCTACTGCTTTATTCCCACCGAACGCGATGACATTTTTTCTGTACCATCCCCCGATCGCCTCGAATTGCGCGAAGGCTGGCGAAAACTATTACCCCAGGCGATAGCCGATCCGCTGAACCTACAAAGCTGGAAAGGCGGGCGAGGGCACGCTGTTTCTGCACTCGAAATGCCGCAAGCGCCTCTGCAACCGGGATGGGATTGTCCGCAAGCGCCAGAAATACCTGCCAAAGAAATTATCTGGAAAAGTGAACGGTTGAAAAA  >  CP009273/608429‑608722

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