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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273335,2470CA91.2% 519.4 / 45.4 206N404K (AAC→AAA) yahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base C (11/7);  new base A (108/80);  total (119/87)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

Only 100 of 206 total aligned reads displayed.

ACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTG  >  CP009273/335102‑335391
                                                                                                                                                 |                                                                                                                                                
aCAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGt                                                                                                                                                >  2:1709354/1‑148 (MQ=255)
aCAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGt                                                                                                                                                >  2:1301269/1‑148 (MQ=255)
aCAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGt                                                                                                                                                >  2:2605856/1‑148 (MQ=255)
 cAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTc                                                                                                                                               >  1:2454945/1‑148 (MQ=255)
   gcgcCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAc                                                                                                                                             <  2:2361046/148‑1 (MQ=255)
      ccATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAt                                                                                                                                          >  1:1785210/1‑148 (MQ=255)
      ccATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAt                                                                                                                                          <  2:2257387/148‑1 (MQ=255)
       cATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAtt                                                                                                                                         <  1:294402/148‑1 (MQ=255)
         tCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTcc                                                                                                                                       >  1:2221968/1‑148 (MQ=255)
          cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                      <  2:1763623/148‑1 (MQ=255)
            ccGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgct                                                                                                                                    >  1:2446228/1‑148 (MQ=255)
            ccGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgct                                                                                                                                    >  1:4857/1‑148 (MQ=255)
             cGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGGCACCATTCCgctg                                                                                                                                   <  2:1154827/148‑1 (MQ=255)
                aaCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTc                                                                                                                                >  1:1345651/1‑148 (MQ=255)
                aaCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTc                                                                                                                                >  2:465564/1‑148 (MQ=255)
                aaCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTc                                                                                                                                >  1:2198278/1‑148 (MQ=255)
                aaCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTc                                                                                                                                >  1:214005/1‑148 (MQ=255)
                 aCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCg                                                                                                                               >  1:2920363/1‑148 (MQ=255)
                  ccTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCgg                                                                                                                              <  2:2018933/148‑1 (MQ=255)
                   cTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGt                                                                                                                             >  1:1370306/1‑148 (MQ=255)
                      cGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtttt                                                                                                                          <  2:1020928/148‑1 (MQ=255)
                       ggTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTa                                                                                                                         <  2:1544199/148‑1 (MQ=255)
                       ggTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTa                                                                                                                         <  2:696518/148‑1 (MQ=255)
                          aTTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATgg                                                                                                                      >  1:2813583/1‑148 (MQ=255)
                          aTTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATgg                                                                                                                      >  1:761110/1‑148 (MQ=255)
                             ggCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCg                                                                                                                   <  1:10104/148‑1 (MQ=255)
                                  aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                              >  2:112222/1‑148 (MQ=255)
                                  aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                              >  1:2717660/1‑148 (MQ=255)
                                       cTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGc                                                                                                         >  2:1476735/1‑148 (MQ=255)
                                        tATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCa                                                                                                        >  1:1325002/1‑148 (MQ=255)
                                            gCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCg                                                                                                    <  1:2686606/148‑1 (MQ=255)
                                            gCGACCGCGCCCGCTATCGTCGCGCGGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCg                                                                                                    >  2:1768070/1‑148 (MQ=255)
                                               aCCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACa                                                                                                 <  2:2689869/148‑1 (MQ=255)
                                                   cgcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAc                                                                                             <  1:1086548/148‑1 (MQ=255)
                                                    gcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAcc                                                                                            <  1:881728/148‑1 (MQ=255)
                                                    gcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAcc                                                                                            <  1:2535901/148‑1 (MQ=255)
                                                      ccGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                                                                          <  2:2919472/148‑1 (MQ=255)
                                                        gCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgcg                                                                                        <  2:2472058/148‑1 (MQ=255)
                                                         cTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                                                       <  1:2814563/148‑1 (MQ=255)
                                                             cgtcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                                                   >  1:229924/1‑148 (MQ=255)
                                                             cgtcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                                                   >  1:230189/1‑148 (MQ=255)
                                                                cgcgCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagc                                                                                >  2:1324901/1‑148 (MQ=255)
                                                                cgcgCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagc                                                                                >  2:2707387/1‑148 (MQ=255)
                                                                       gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                                                                         >  1:3162551/1‑148 (MQ=255)
                                                                       gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                                                                         <  2:639852/148‑1 (MQ=255)
                                                                              gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                                                  >  1:2583987/1‑148 (MQ=255)
                                                                               cACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                                                 >  2:2127250/1‑148 (MQ=255)
                                                                                   gTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGa                                                             >  1:643156/1‑148 (MQ=255)
                                                                                   gTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGa                                                             <  2:747362/148‑1 (MQ=255)
                                                                                    tGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt                                                            <  1:1431246/148‑1 (MQ=255)
                                                                                     ggAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                                                           >  2:402410/1‑148 (MQ=255)
                                                                                       aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa                                                         >  1:2784647/1‑148 (MQ=255)
                                                                                         gaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAaca                                                       >  2:2137363/1‑148 (MQ=255)
                                                                                         gaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAaca                                                       >  1:1736630/1‑148 (MQ=255)
                                                                                             cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                                                   <  1:1436127/148‑1 (MQ=255)
                                                                                             cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                                                   >  2:2782786/1‑148 (MQ=255)
                                                                                             cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGAAACGGCATTCTGCCGGTGATCAACACCGc                                                   <  2:573577/148‑1 (MQ=255)
                                                                                              tATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGcc                                                  >  2:2870131/1‑148 (MQ=255)
                                                                                                ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                                                >  2:1083203/1‑148 (MQ=255)
                                                                                                ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                                                <  2:214005/148‑1 (MQ=255)
                                                                                                 tGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATc                                               >  1:1197868/1‑148 (MQ=255)
                                                                                                 tGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATc                                               >  2:727378/1‑148 (MQ=255)
                                                                                                     ttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGccc                                           <  2:2221968/148‑1 (MQ=255)
                                                                                                      tctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa                                          >  2:1093626/1‑148 (MQ=255)
                                                                                                        tcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa                                        <  1:336352/148‑1 (MQ=255)
                                                                                                        tcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa                                        <  1:336407/148‑1 (MQ=255)
                                                                                                           cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                     <  1:1917434/148‑1 (MQ=255)
                                                                                                           cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                     <  1:387246/148‑1 (MQ=255)
                                                                                                            gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGa                                    <  1:465564/148‑1 (MQ=255)
                                                                                                            gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGa                                    <  2:2851564/148‑1 (MQ=255)
                                                                                                             tCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGAt                                   <  1:1217853/148‑1 (MQ=255)
                                                                                                              cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                  >  1:2779204/1‑148 (MQ=255)
                                                                                                              cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                  >  2:1353045/1‑148 (MQ=255)
                                                                                                                gATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCg                                >  1:29824/1‑148 (MQ=255)
                                                                                                                       gAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCg                         >  2:3007493/1‑148 (MQ=255)
                                                                                                                           tCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCAt                     <  1:1768040/148‑1 (MQ=255)
                                                                                                                            cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATg                    <  2:1223239/148‑1 (MQ=255)
                                                                                                                            cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATg                    <  1:166070/148‑1 (MQ=255)
                                                                                                                             aCCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGa                   <  2:2198278/148‑1 (MQ=255)
                                                                                                                             aCCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGa                   <  2:1345651/148‑1 (MQ=255)
                                                                                                                             aCCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGa                   >  2:2771580/1‑148 (MQ=255)
                                                                                                                                cTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTg                <  2:1405864/148‑1 (MQ=255)
                                                                                                                                   ggTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGgcg             <  1:2060438/148‑1 (MQ=255)
                                                                                                                                    gTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAGTGgcgc            >  2:1763902/1‑148 (MQ=255)
                                                                                                                                     tGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGcgcg           >  1:2107626/1‑148 (MQ=255)
                                                                                                                                     tGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGcgcg           >  1:2419397/1‑148 (MQ=255)
                                                                                                                                     tGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGcgcg           >  1:443109/1‑148 (MQ=255)
                                                                                                                                      gAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCgg          >  2:1036919/1‑148 (MQ=255)
                                                                                                                                      gAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCgg          >  2:2390825/1‑148 (MQ=255)
                                                                                                                                      gAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCgg          >  1:932156/1‑148 (MQ=255)
                                                                                                                                        aaaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGc        >  1:2252340/1‑148 (MQ=255)
                                                                                                                                        aaaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGc        >  1:419054/1‑148 (MQ=255)
                                                                                                                                        aaaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGc        >  1:419003/1‑148 (MQ=255)
                                                                                                                                         aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa       >  2:1456920/1‑148 (MQ=255)
                                                                                                                                         aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa       >  2:1760787/1‑148 (MQ=255)
                                                                                                                                         aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa       <  2:590653/148‑1 (MQ=255)
                                                                                                                                          aCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAt      >  1:767694/1‑148 (MQ=255)
                                                                                                                                            cccAAAGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCGGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATtg    <  2:2168266/148‑1 (MQ=255)
                                                                                                                                             ccAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATtgt   <  1:1652223/148‑1 (MQ=255)
                                                                                                                                              cAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATtgtg  <  1:3069761/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTG  >  CP009273/335102‑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: ‑