Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,335,356 A→T 100% H181L (CAT→CTT)  yciM → envelope integrity maintenance protein; EnvC‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,335,3560AT97.5% 594.7 / 2.4 183H181L (CAT→CTT) yciMenvelope integrity maintenance protein; EnvC‑interacting protein
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base T (90/88);  total (94/89)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 192 total aligned reads displayed.

CTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAGGAGA  >  CP009273/1335211‑1335500
                                                                                                                                                 |                                                                                                                                                
cTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGACATTGCCCttt                                                                                                                                                >  2:3379979/1‑148 (MQ=255)
cTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCttt                                                                                                                                                >  1:2028218/1‑148 (MQ=255)
cTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCttt                                                                                                                                                <  2:2018440/148‑1 (MQ=255)
 tGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCtttt                                                                                                                                               <  1:2166220/148‑1 (MQ=255)
  gACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTc                                                                                                                                              <  1:2960357/148‑1 (MQ=255)
   aCCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCACTTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCt                                                                                                                                             >  2:3203493/1‑148 (MQ=255)
   aCCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCt                                                                                                                                             >  2:781211/1‑148 (MQ=255)
    ccGATGAACCTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTa                                                                                                                                            <  2:3254177/148‑1 (MQ=255)
     cGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTTAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTAc                                                                                                                                           <  2:627263/148‑1 (MQ=255)
      gATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACt                                                                                                                                          <  1:1039/148‑1 (MQ=255)
       aTGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACtg                                                                                                                                         <  1:2668891/148‑1 (MQ=255)
         gAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACtgtg                                                                                                                                       <  1:3295522/148‑1 (MQ=255)
         gAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACtgtg                                                                                                                                       <  1:309146/148‑1 (MQ=255)
          aaaCTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGa                                                                                                                                      <  1:2762827/148‑1 (MQ=255)
          aaaCTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGa                                                                                                                                      <  1:141779/148‑1 (MQ=255)
           aaCTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAg                                                                                                                                     <  2:3715/148‑1 (MQ=255)
           aaCTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAg                                                                                                                                     <  1:1977583/148‑1 (MQ=255)
            aCTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGt                                                                                                                                    >  1:3388597/1‑148 (MQ=255)
               gACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAg                                                                                                                                 <  1:2297224/148‑1 (MQ=255)
                aCTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGc                                                                                                                                <  2:3279795/148‑1 (MQ=255)
                aCTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGc                                                                                                                                >  2:1042738/1‑148 (MQ=255)
                aCTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGc                                                                                                                                >  2:1044211/1‑148 (MQ=255)
                aCTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGc                                                                                                                                <  2:3279823/148‑1 (MQ=255)
                   tCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCt                                                                                                                             >  1:3043213/1‑148 (MQ=255)
                          tGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCat                                                                                                                      >  1:3102784/1‑148 (MQ=255)
                            gcgcgcTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCatat                                                                                                                    <  1:2678337/148‑1 (MQ=255)
                            gcgcgcTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCatat                                                                                                                    <  2:2294791/148‑1 (MQ=255)
                            gcgcgcTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCatat                                                                                                                    <  2:3081326/148‑1 (MQ=255)
                                 cTGCAAGTGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCa                                                                                                               <  1:1854874/148‑1 (MQ=255)
                                 cTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCa                                                                                                               <  1:1854280/148‑1 (MQ=255)
                                 cTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCa                                                                                                               >  1:1560545/1‑148 (MQ=255)
                                     aaCAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGcga                                                                                                           <  2:1104042/148‑1 (MQ=255)
                                      aCAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGcgac                                                                                                          >  2:681618/1‑148 (MQ=255)
                                      aCAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGcgac                                                                                                          >  1:3117969/1‑148 (MQ=255)
                                       cAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGcgacg                                                                                                         >  1:1935832/1‑148 (MQ=255)
                                       cAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGcgacg                                                                                                         >  1:810481/1‑148 (MQ=255)
                                          ttGCTAAAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGATTTAGCCCTGCAGCATATGGCCAGCGACGAtc                                                                                                      <  1:1286622/148‑1 (MQ=255)
                                           tGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGAtct                                                                                                     >  2:2178967/1‑148 (MQ=255)
                                             cTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCg                                                                                                   <  2:3059134/148‑1 (MQ=255)
                                             cTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCg                                                                                                   <  2:2696155/148‑1 (MQ=255)
                                                cAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATc                                                                                                >  2:1794785/1‑148 (MQ=255)
                                                  aaTCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGTCAGCGACGATCTCGATCGt                                                                                              >  2:34894/1‑148 (MQ=255)
                                                   aTCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTg                                                                                             <  1:3420821/148‑1 (MQ=255)
                                                   aTCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTg                                                                                             <  1:611193/148‑1 (MQ=255)
                                                    tCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGc                                                                                            >  2:855096/1‑148 (MQ=255)
                                                     cTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGcc                                                                                           <  1:723597/148‑1 (MQ=255)
                                                     cTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGcc                                                                                           >  2:324499/1‑148 (MQ=255)
                                                        ccAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATg                                                                                        <  1:107029/148‑1 (MQ=255)
                                                          aGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGAc                                                                                      <  1:509894/148‑1 (MQ=255)
                                                            gCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCt                                                                                    >  2:1385359/1‑148 (MQ=255)
                                                            gCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCt                                                                                    <  1:289147/148‑1 (MQ=255)
                                                               aCCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGc                                                                                 >  1:1503819/1‑148 (MQ=255)
                                                               aCCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGc                                                                                 <  2:1006709/148‑1 (MQ=255)
                                                                ccAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCt                                                                                >  1:620315/1‑148 (MQ=255)
                                                                 cAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTa                                                                               >  1:2360168/1‑148 (MQ=255)
                                                                  agtgagtgGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTaa                                                                              <  2:2852471/148‑1 (MQ=255)
                                                                  agtgagtgGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTaa                                                                              >  2:2201124/1‑148 (MQ=255)
                                                                     gagtgGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTaaaaa                                                                           >  1:3311624/1‑148 (MQ=255)
                                                                     gagtgGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTaaaaa                                                                           >  1:3311630/1‑148 (MQ=255)
                                                                          gCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAggg                                                                      <  1:2717904/148‑1 (MQ=255)
                                                                          gCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAggg                                                                      <  2:2444842/148‑1 (MQ=255)
                                                                             gAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGggcg                                                                   >  2:2539188/1‑148 (MQ=255)
                                                                             gAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGggcg                                                                   >  2:2800545/1‑148 (MQ=255)
                                                                             gAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGggcg                                                                   <  1:44448/148‑1 (MQ=255)
                                                                                 gCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGggcggcgg                                                               >  2:186918/1‑148 (MQ=255)
                                                                                   aaTTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCa                                                             >  1:750098/1‑148 (MQ=255)
                                                                                        aTGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATaa                                                        <  1:2173325/148‑1 (MQ=255)
                                                                                          gTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATaaaa                                                      <  1:2042316/148‑1 (MQ=255)
                                                                                                gAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAgcg                                                <  1:2453724/148‑1 (MQ=255)
                                                                                                  aCGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGcc                                              >  2:693838/1‑148 (MQ=255)
                                                                                                  aCGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGcc                                              <  2:2609878/148‑1 (MQ=255)
                                                                                                         gTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTAt                                       <  2:1935832/148‑1 (MQ=255)
                                                                                                          tGAAGCTGGTTAAATATAAACAGCGCGTCTAAATTGCCATTTTATACTGTGAGTTAGCCCTGCAGCATAGGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATc                                      <  2:2159748/148‑1 (MQ=255)
                                                                                                            aaGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCa                                    >  2:1018093/1‑148 (MQ=255)
                                                                                                              gCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATa                                  >  1:1152365/1‑148 (MQ=255)
                                                                                                                 ggttAAATATAAACAGCGCGTCTAAATTGCCCTTTTCTACTGTTTGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAatg                               <  1:3520700/145‑1 (MQ=255)
                                                                                                                    tAAAGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAatgatg                            >  2:1987752/1‑148 (MQ=255)
                                                                                                                      aaGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATggg                          >  2:1729789/1‑148 (MQ=255)
                                                                                                                       ataTAAACAGCGCGTCTAATTTGCCTTTTTCTACTGTGTGTTAGCCCTGCAGAATAGGCCCCGCGACGATCGCGATGGTGCCATGACGTTGCAAAAAAAAGGGGGGGCGGCATATAAAAACCGCGCCCGCGTATCCATAATGATGGGa                         <  1:1042756/146‑1 (MQ=255)
                                                                                                                       aGATAAACCGCGCGTCTAATTTGCCCTTTTCTACTGTGTTTTAGCCCTGCAGCATATGGCCCCCGACGATCTCGATCGTGCAATGACCTTGATAAAAAATGGGGCGGGGGCATATAAAAACAGCGCCCGCGTATCCATAATGATGGGa                         <  1:1042279/148‑1 (MQ=255)
                                                                                                                       aGATAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGa                         >  1:1669513/1‑148 (MQ=255)
                                                                                                                          tAAACAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGAcgc                      >  1:992380/1‑148 (MQ=255)
                                                                                                                             aCAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCgtg                   <  2:1118074/148‑1 (MQ=255)
                                                                                                                             aCAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCgtg                   >  2:1058151/1‑148 (MQ=255)
                                                                                                                             aCAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCgtg                   >  2:3332917/1‑148 (MQ=255)
                                                                                                                              cAGCGCGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCgtgt                  >  1:2812088/1‑148 (MQ=255)
                                                                                                                                gcgcGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGttt                >  1:1344948/1‑148 (MQ=255)
                                                                                                                                gcgcGTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGttt                >  1:2467830/1‑148 (MQ=255)
                                                                                                                                 cgcgTCTAAATTCCCCTTTTCTACTGTTAGTTAGCCCTGCAGCATTGGGCCAGCGACTATCGCGATCGTGCCATGACCTTGATAAAAAACGGGGGGGGGGCATATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTa               <  1:415607/148‑1 (MQ=255)
                                                                                                                                 cgcgTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTa               <  2:241481/148‑1 (MQ=255)
                                                                                                                                 cgcgTCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTa               <  2:595443/148‑1 (MQ=255)
                                                                                                                                     tCGAAATTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGc           >  2:1606640/1‑148 (MQ=255)
                                                                                                                                         aaTTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGaaa       >  2:477968/1‑148 (MQ=255)
                                                                                                                                          aTTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAg      >  2:1897026/1‑148 (MQ=255)
                                                                                                                                          aTTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAg      >  2:2371078/1‑148 (MQ=255)
                                                                                                                                          aTTGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAg      >  2:2775785/1‑148 (MQ=255)
                                                                                                                                           ttGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAgg     >  2:2149859/1‑148 (MQ=255)
                                                                                                                                            tGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTANAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAGga    >  1:1808872/1‑148 (MQ=255)
                                                                                                                                            tGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAGga    >  1:943830/1‑148 (MQ=255)
                                                                                                                                            tGCCCTTTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAGga    >  1:2479848/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGCGGCAGATAAAAACAGCGCCCGCGTATCCATAATGATGGGACGCGTGTTTATGGCGAAAGGAGA  >  CP009273/1335211‑1335500

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: