Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 335,247 C→A 94.3% N404K (AAC→AAA)  yahG → DUF1116 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273335,2470CA94.3% 525.6 / 20.3 197N404K (AAC→AAA) yahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base C (8/4);  new base A (110/74);  total (118/79)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

AGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTG  >  CP009273/335104‑335391
                                                                                                                                               |                                                                                                                                                
aGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCa                                                                                                                                              >  2:1397516/1‑148 (MQ=255)
 gcgcCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAc                                                                                                                                             >  1:236164/1‑148 (MQ=255)
  cgcCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAcc                                                                                                                                            >  1:2215350/1‑148 (MQ=255)
      aTCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTc                                                                                                                                        >  2:1782630/1‑148 (MQ=255)
        cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                      >  1:1411628/1‑148 (MQ=255)
        cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                      >  2:2336503/1‑148 (MQ=255)
          ccGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgct                                                                                                                                    >  1:530151/1‑148 (MQ=255)
             aaaCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgctgct                                                                                                                                 <  1:2562467/148‑1 (MQ=255)
                 cTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGt                                                                                                                             <  1:595257/148‑1 (MQ=255)
                 cTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGt                                                                                                                             <  1:937356/148‑1 (MQ=255)
                  tACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtt                                                                                                                            >  2:14000/1‑148 (MQ=255)
                  tACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtt                                                                                                                            >  2:571546/1‑148 (MQ=255)
                   aCGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGttt                                                                                                                           >  2:2381938/1‑148 (MQ=255)
                    cGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtttt                                                                                                                          <  2:1729593/148‑1 (MQ=255)
                      gTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTAt                                                                                                                        >  2:2158754/1‑148 (MQ=255)
                       tATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATg                                                                                                                       >  1:3231538/1‑148 (MQ=255)
                       tATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATg                                                                                                                       >  1:3231554/1‑148 (MQ=255)
                       tATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATg                                                                                                                       >  1:2130559/1‑148 (MQ=255)
                          tGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGc                                                                                                                    <  2:1570975/148‑1 (MQ=255)
                           ggCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCTAAATCAACCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCg                                                                                                                   <  1:2455981/148‑1 (MQ=255)
                            gCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGt                                                                                                                  >  1:2787686/1‑148 (MQ=255)
                             cGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTg                                                                                                                 <  1:1700872/148‑1 (MQ=255)
                              ggATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGc                                                                                                                <  2:1503019/148‑1 (MQ=255)
                              ggATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGc                                                                                                                <  1:3404283/148‑1 (MQ=255)
                                aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                              <  2:1767040/148‑1 (MQ=255)
                                aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                              <  2:2559974/148‑1 (MQ=255)
                                  ttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                                                                            >  2:3115396/1‑148 (MQ=255)
                                   tGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCg                                                                                                           <  1:346662/148‑1 (MQ=255)
                                       aTGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCaa                                                                                                       >  2:1033889/1‑148 (MQ=255)
                                        tGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAAt                                                                                                      >  2:2187181/1‑148 (MQ=255)
                                         ggCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATc                                                                                                     >  2:2741115/1‑148 (MQ=255)
                                          gCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCg                                                                                                    <  2:2139663/148‑1 (MQ=255)
                                           cGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGa                                                                                                   <  1:2265741/148‑1 (MQ=255)
                                           cGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGa                                                                                                   >  2:2778646/1‑148 (MQ=255)
                                           cGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGa                                                                                                   >  1:145833/1‑148 (MQ=255)
                                               cgcgCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                                               <  2:3411901/148‑1 (MQ=255)
                                                gcgcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCa                                                                                              <  2:341474/148‑1 (MQ=255)
                                                      gCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgcg                                                                                        <  2:1706382/148‑1 (MQ=255)
                                                       cTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                                                       <  2:871062/148‑1 (MQ=255)
                                                           cgtcgCGCTGGTGGTCGGCACGGTGTAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                                                   <  2:3526939/148‑1 (MQ=255)
                                                           cgtcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                                                   >  2:446661/1‑148 (MQ=255)
                                                           cgtcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                                                   >  1:2499289/1‑148 (MQ=255)
                                                                 gcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAACTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                                                                             <  2:2579104/148‑1 (MQ=255)
                                                                 gcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                                                                             <  2:2710115/148‑1 (MQ=255)
                                                                   tggtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCACTCGACATCACCCGCGTGGGTAGCAGCgg                                                                           >  2:718285/1‑148 (MQ=255)
                                                                     gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                                                                         >  1:3261098/1‑148 (MQ=255)
                                                                     gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATAACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                                                                         <  1:1847558/148‑1 (MQ=255)
                                                                        ggcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTc                                                                      <  2:1999878/148‑1 (MQ=255)
                                                                          cggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTg                                                                    <  2:3353173/148‑1 (MQ=255)
                                                                            gcACGGTGGAAGAAGCTATTGATTTTTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                                                  >  2:759271/1‑148 (MQ=255)
                                                                            gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCTGTGAAAACCCCAAAGTCACCATTCAGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCAGCGTGGGTAGCAGCGGCATTCTGcc                                                                  >  1:425832/1‑148 (MQ=255)
                                                                            gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                                                  >  1:201979/1‑148 (MQ=255)
                                                                               cGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGt                                                               >  1:635407/1‑148 (MQ=255)
                                                                                 gTGGAAGAAGCTATTGATTTCTCCCGTCTTATGCGCGAAATCACCCTCGGTGAAAACCCCAATGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGa                                                             <  2:1422756/148‑1 (MQ=255)
                                                                                     aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa                                                         >  2:1309159/1‑148 (MQ=255)
                                                                                     aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa                                                         >  2:2970124/1‑148 (MQ=255)
                                                                                      agaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAac                                                        <  1:714247/148‑1 (MQ=255)
                                                                                      agaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAac                                                        <  1:714457/148‑1 (MQ=255)
                                                                                       gaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAaca                                                       >  1:266014/1‑148 (MQ=255)
                                                                                          gCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                                                    >  1:2827075/1‑148 (MQ=255)
                                                                                           cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                                                   >  2:1143945/1‑148 (MQ=255)
                                                                                           cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                                                   >  2:3037120/1‑148 (MQ=255)
                                                                                           cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                                                   >  1:3367465/1‑148 (MQ=255)
                                                                                             aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTCGCAGCGGCATTCTGCCGGTGATCAACACCGCCa                                                 >  2:2593150/1‑148 (MQ=255)
                                                                                             aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCa                                                 <  2:855873/148‑1 (MQ=255)
                                                                                             aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCa                                                 <  2:2263988/148‑1 (MQ=255)
                                                                                                      tcCCGTCAGATGCTCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa                                        >  2:1135751/1‑148 (MQ=255)
                                                                                                       cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaa                                       <  1:1679285/148‑1 (MQ=255)
                                                                                                        ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaaa                                      >  2:1547659/1‑148 (MQ=255)
                                                                                                        ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaaa                                      >  2:2197701/1‑148 (MQ=255)
                                                                                                        ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaaa                                      >  2:100135/1‑148 (MQ=255)
                                                                                                         cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                     >  1:3379551/1‑148 (MQ=255)
                                                                                                          gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGa                                    <  2:2423761/148‑1 (MQ=255)
                                                                                                            cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                  >  2:10258/1‑148 (MQ=255)
                                                                                                            cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                  <  2:324123/148‑1 (MQ=255)
                                                                                                             aGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGc                                 <  2:3489294/148‑1 (MQ=255)
                                                                                                               aTGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCgg                               >  1:1953958/1‑148 (MQ=255)
                                                                                                                 gcgcGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGc                             >  2:811164/1‑148 (MQ=255)
                                                                                                                  cgcgAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCg                            <  1:1198696/148‑1 (MQ=255)
                                                                                                                      aaaTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCgg                        >  1:1698888/1‑148 (MQ=255)
                                                                                                                       aaTCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGc                       >  1:293040/1‑148 (MQ=255)
                                                                                                                         tCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCAt                     <  1:1010321/148‑1 (MQ=255)
                                                                                                                          cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATg                    <  1:251745/148‑1 (MQ=255)
                                                                                                                          cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATg                    >  1:3464743/1‑148 (MQ=255)
                                                                                                                          cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATg                    <  1:3454652/148‑1 (MQ=255)
                                                                                                                          cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATg                    <  1:2965685/148‑1 (MQ=255)
                                                                                                                           aCCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGa                   >  2:1409591/1‑148 (MQ=255)
                                                                                                                             ccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt                 <  2:1056813/148‑1 (MQ=255)
                                                                                                                             ccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt                 <  2:2488644/148‑1 (MQ=255)
                                                                                                                               tCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTgg               >  2:2692563/1‑148 (MQ=255)
                                                                                                                                  gTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAAATGCGGGCGTCGGCATGATTGgcgc            >  2:2848325/1‑148 (MQ=255)
                                                                                                                                  gTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAAATGCGGGCGTCGGCATGATTGgcgc            >  2:2847910/1‑148 (MQ=255)
                                                                                                                                    gAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCgg          >  1:2910647/1‑148 (MQ=255)
                                                                                                                                     aaaaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCggg         >  2:938791/1‑148 (MQ=255)
                                                                                                                                      aaaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGc        <  1:2512063/148‑1 (MQ=255)
                                                                                                                                       aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa       <  1:831530/148‑1 (MQ=255)
                                                                                                                                       aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa       <  1:1397535/148‑1 (MQ=255)
                                                                                                                                          cccAAAGTCTCCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCTACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATtg    <  1:761308/148‑1 (MQ=255)
                                                                                                                                          cccAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATtg    <  1:1170674/148‑1 (MQ=255)
                                                                                                                                            cAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATtgtg  >  2:2916567/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTG  >  CP009273/335104‑335391

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: ‑