Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,953,876 C→A 100% P76P (CCG→CCT yecM ← putative metal‑binding enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,953,8760CA96.5% 612.0 / 12.0 203P76P (CCG→CCTyecMputative metal‑binding enzyme
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (105/91);  total (108/95)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 203 total aligned reads displayed.

GCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGG  >  CP009273/1953731‑1954017
                                                                                                                                                 |                                                                                                                                             
gCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCgg                                                                                                                                             <  1:3371632/148‑1 (MQ=255)
gCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAgg                                                                                                                                             <  1:3615287/148‑1 (MQ=255)
gCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAgg                                                                                                                                             >  1:911985/1‑148 (MQ=255)
gCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAgg                                                                                                                                             >  1:1466306/1‑148 (MQ=255)
gCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCAAAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAgg                                                                                                                                             <  2:3296379/148‑1 (MQ=255)
 cGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGt                                                                                                                                            >  2:2197366/1‑148 (MQ=255)
 cGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGt                                                                                                                                            <  1:1979329/148‑1 (MQ=255)
 cGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGt                                                                                                                                            <  1:1147230/148‑1 (MQ=255)
  gTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTc                                                                                                                                           <  1:1349164/148‑1 (MQ=255)
   tttCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCt                                                                                                                                          >  1:2028255/1‑148 (MQ=255)
   tttCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCt                                                                                                                                          >  1:3098533/1‑148 (MQ=255)
   tttCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCt                                                                                                                                          >  2:2846556/1‑148 (MQ=255)
    ttCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCTCGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTg                                                                                                                                         <  2:2731285/148‑1 (MQ=255)
    ttCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTg                                                                                                                                         <  1:3113391/148‑1 (MQ=255)
    ttCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTg                                                                                                                                         >  2:3432738/1‑148 (MQ=255)
      cTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGcc                                                                                                                                       <  2:1949248/148‑1 (MQ=255)
       tGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCa                                                                                                                                      <  2:2915229/148‑1 (MQ=255)
         gATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCAtt                                                                                                                                    >  2:3658427/1‑148 (MQ=255)
          aTCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACCAATAGGTCTGCCATTg                                                                                                                                   >  2:2219882/1‑148 (MQ=255)
            cACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGAt                                                                                                                                 <  1:3463514/148‑1 (MQ=255)
            cACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGAt                                                                                                                                 <  2:3653767/148‑1 (MQ=255)
             aCCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATc                                                                                                                                <  1:1207126/148‑1 (MQ=255)
             aCCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTATGCCATTGATc                                                                                                                                <  1:2362644/148‑1 (MQ=255)
              cccGGCAGAACGATTTCAATATGATCCCAACCTTCGTGTGGGTAACGTTTTTCCCCCGGCCACGGTAGCACCACAATAGAAAACTGCCAGTGCTCAACCAGTACCGGTTCAAGCACTTTAAACAGACAAATCGGTCTGACATTGATCa                                                                                                                               >  1:2478567/1‑148 (MQ=255)
               ccGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCat                                                                                                                              >  2:2696624/1‑148 (MQ=255)
               ccGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCat                                                                                                                              >  1:1138775/1‑148 (MQ=255)
               ccGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCat                                                                                                                              >  2:587223/1‑148 (MQ=255)
                     gAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTc                                                                                                                        >  1:989464/1‑148 (MQ=255)
                     gAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTc                                                                                                                        >  1:905060/1‑148 (MQ=255)
                      aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCt                                                                                                                       >  1:3301052/1‑148 (MQ=255)
                      aaCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCt                                                                                                                       >  2:734983/1‑148 (MQ=255)
                       aCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTg                                                                                                                      <  2:3385240/148‑1 (MQ=255)
                       aCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTg                                                                                                                      <  1:1175314/148‑1 (MQ=255)
                       aCGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTg                                                                                                                      >  1:1706324/1‑148 (MQ=255)
                        cGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGa                                                                                                                     >  1:2914764/1‑148 (MQ=255)
                        cGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGa                                                                                                                     >  2:3282065/1‑148 (MQ=255)
                         gATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGAc                                                                                                                    >  2:367751/1‑148 (MQ=255)
                         gATTTCAATATGTTCCCAACATTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGAc                                                                                                                    >  2:1551761/1‑148 (MQ=255)
                           tttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaa                                                                                                                  >  1:3045315/1‑148 (MQ=255)
                           tttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaa                                                                                                                  >  1:3710745/1‑148 (MQ=255)
                            ttCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACaaa                                                                                                                 >  2:408647/1‑148 (MQ=255)
                             tCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACaaaa                                                                                                                >  2:1573918/1‑148 (MQ=255)
                               aaTATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGc                                                                                                              <  2:1713860/148‑1 (MQ=255)
                                 tatGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTc                                                                                                            >  2:1682141/1‑148 (MQ=255)
                                 tatGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTc                                                                                                            >  2:1682320/1‑148 (MQ=255)
                                  atGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCg                                                                                                           >  2:2954684/1‑148 (MQ=255)
                                       cccAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCcaca                                                                                                      >  1:11827/1‑148 (MQ=255)
                                        ccAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCacac                                                                                                     >  1:1243618/1‑148 (MQ=255)
                                         cAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACt                                                                                                    >  1:3613485/1‑148 (MQ=255)
                                           aCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGt                                                                                                  >  1:2155575/1‑148 (MQ=255)
                                           aCCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGt                                                                                                  >  2:3696698/1‑148 (MQ=255)
                                             cTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTc                                                                                                >  1:3777310/1‑148 (MQ=255)
                                             cTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTc                                                                                                >  1:2355718/1‑148 (MQ=255)
                                              ttCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCa                                                                                               >  1:3817747/1‑148 (MQ=255)
                                               tCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCaa                                                                                              <  2:869133/148‑1 (MQ=255)
                                               tCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGAAAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACAGTTCaa                                                                                              >  2:2242949/1‑148 (MQ=255)
                                                  tgtgGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAAcc                                                                                           >  2:1327359/1‑148 (MQ=255)
                                                         aaCGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGAGCGCAACCTGTACCGGTTCATGCAATTTAAACATACAAATAGGTCTGCCATTGAACAAATTTTCTGACAAAAGCTAGCCACACTGTTCAAACCCGCGACg                                                                                    >  2:3208243/1‑148 (MQ=255)
                                                              ttttCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGc                                                                               >  1:714692/1‑148 (MQ=255)
                                                                 tCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGtt                                                                            <  1:766723/148‑1 (MQ=255)
                                                                 tCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGtt                                                                            <  1:3294911/148‑1 (MQ=255)
                                                                     cGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTcagc                                                                        <  1:2731044/148‑1 (MQ=255)
                                                                            cGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCg                                                                 >  1:543085/1‑148 (MQ=255)
                                                                            cGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCg                                                                 >  1:3758459/1‑148 (MQ=255)
                                                                              gTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGtt                                                               <  1:3658427/148‑1 (MQ=255)
                                                                               tAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGttt                                                              <  1:2197366/148‑1 (MQ=255)
                                                                                 gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTg                                                            >  2:1149059/1‑148 (MQ=255)
                                                                                 gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTg                                                            >  1:1626877/1‑148 (MQ=255)
                                                                                 gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTg                                                            <  2:948920/148‑1 (MQ=255)
                                                                                 gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTg                                                            <  1:3403096/148‑1 (MQ=255)
                                                                                 gCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTg                                                            >  2:1883991/1‑148 (MQ=255)
                                                                                  cTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGa                                                           >  2:1935012/1‑148 (MQ=255)
                                                                                    ccACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATg                                                         <  1:3050692/148‑1 (MQ=255)
                                                                                     cacaATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATgg                                                        <  1:3429677/148‑1 (MQ=255)
                                                                                        aaTAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCaa                                                     <  1:1906601/148‑1 (MQ=255)
                                                                                          tAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACg                                                   >  1:1047642/1‑148 (MQ=255)
                                                                                             aaaaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCaa                                                <  1:2428228/148‑1 (MQ=255)
                                                                                               aaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAg                                              >  2:2390222/1‑148 (MQ=255)
                                                                                               aaCTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAg                                              >  1:3029041/1‑148 (MQ=255)
                                                                                                 cTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGaa                                            <  2:3314787/148‑1 (MQ=255)
                                                                                                    ccAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAAtg                                         <  1:168359/148‑1 (MQ=255)
                                                                                                               ccTGTACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGt                              <  2:3211929/148‑1 (MQ=255)
                                                                                                                   tACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGc                          <  1:1660140/148‑1 (MQ=255)
                                                                                                                   tACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGc                          >  1:2162207/1‑148 (MQ=255)
                                                                                                                   tACCGGTTCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAAAGTGATCGGCAGTCAGc                          >  1:3511924/1‑148 (MQ=255)
                                                                                                                         ttCATGCAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTg                    <  2:2581833/148‑1 (MQ=255)
                                                                                                                               cAATTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTc              >  2:2089862/1‑148 (MQ=255)
                                                                                                                                 aTTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCaa            <  2:2028255/148‑1 (MQ=255)
                                                                                                                                 aTTTAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCaa            <  2:2771170/148‑1 (MQ=255)
                                                                                                                                  tttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAc           >  2:119005/1‑148 (MQ=255)
                                                                                                                                  tttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAc           >  2:10994/1‑148 (MQ=255)
                                                                                                                                  tttAAACAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAAc           <  1:2119894/148‑1 (MQ=255)
                                                                                                                                      aaCAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCa       >  2:3561417/1‑148 (MQ=255)
                                                                                                                                       aCAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAg      <  1:2846556/148‑1 (MQ=255)
                                                                                                                                       aCAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAg      <  1:1682164/148‑1 (MQ=255)
                                                                                                                                        cAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGa     <  2:2261584/148‑1 (MQ=255)
                                                                                                                                        cAGACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGa     <  1:3069402/148‑1 (MQ=255)
                                                                                                                                          gACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGAcg   >  2:962014/1‑148 (MQ=255)
                                                                                                                                          gACAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGAcg   <  2:3137431/148‑1 (MQ=255)
                                                                                                                                           aCAAATAGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGAcgg  >  2:3262367/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GCGTTTCTGGATCACCCGGCAGAACGATTTCAATATGTTCCCAACCTTCGTGTGGGTAACGTTTTTCCCCGGGCCACGGTAGCTCCACAATAGAAAACTGCCAGTGCGCAACCTGTACCGGTTCATGCAATTTAAACAGACAAATCGGTCTGCCATTGATCATATTTTCTGACAAAAGCTCGCCACACTGTTCAAACCCGCGACGCCAGCGTTCAGCAGTGGCGTTTTGATGGCAACGCAAAGAAATGTGATCGGCAGTCAGCGGAGTGATATTCAACCCCAGACGG  >  CP009273/1953731‑1954017

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: