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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273335,2470CA97.8% 578.1 / 3.5 202N404K (AAC→AAAyahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (115/82);  total (117/85)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 202 total aligned reads displayed.

CAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGC  >  CP009273/335103‑335392
                                                                                                                                                |                                                                                                                                                 
cAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTc                                                                                                                                                <  1:7924/148‑1 (MQ=255)
 aGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCa                                                                                                                                               <  2:200805/148‑1 (MQ=255)
      cATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAtt                                                                                                                                          >  2:1001976/1‑148 (MQ=255)
       aTCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTc                                                                                                                                         <  1:656732/148‑1 (MQ=255)
       aTCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTc                                                                                                                                         >  1:342812/1‑148 (MQ=255)
        tCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTcc                                                                                                                                        >  1:1582705/1‑148 (MQ=255)
        tCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTcc                                                                                                                                        >  2:2694326/1‑148 (MQ=255)
        tCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTcc                                                                                                                                        >  2:2775701/1‑148 (MQ=255)
         cACCGAAACCTACGGTATTTGCGGATTTGCTATGGCGACCGCTCCCGCTATCGTAGCGCTGGTGGGCGGCACTGTGGAAGAAGCTATTGAATTCTCCCGTCAGATGCGCGAAATCACCCTCGGAGAAAACCACAAAGTCACCATTCCg                                                                                                                                       >  1:3217699/1‑148 (MQ=255)
         cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                       <  2:2197240/148‑1 (MQ=255)
         cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                       >  1:340998/1‑148 (MQ=255)
            cGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgctg                                                                                                                                    >  2:2862594/1‑148 (MQ=255)
             gAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgctgc                                                                                                                                   <  1:2943336/148‑1 (MQ=255)
             gAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgctgc                                                                                                                                   <  1:2942618/148‑1 (MQ=255)
               aaCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTc                                                                                                                                 >  2:1678581/1‑148 (MQ=255)
                aCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCg                                                                                                                                >  2:411802/1‑148 (MQ=255)
                aCCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCg                                                                                                                                >  1:1204020/1‑148 (MQ=255)
                   tACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtt                                                                                                                             <  1:2516049/148‑1 (MQ=255)
                   tACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtt                                                                                                                             >  2:2966577/1‑148 (MQ=255)
                       gTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTAt                                                                                                                         <  2:1289371/148‑1 (MQ=255)
                        tATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATg                                                                                                                        >  2:1446496/1‑148 (MQ=255)
                        tATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATg                                                                                                                        >  1:480340/1‑148 (MQ=255)
                             gCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGt                                                                                                                   >  2:407496/1‑148 (MQ=255)
                                gATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGcc                                                                                                                >  2:1483054/1‑148 (MQ=255)
                                  tttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGt                                                                                                              >  1:1813286/1‑148 (MQ=255)
                                   ttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                                                                             <  2:439047/148‑1 (MQ=255)
                                   ttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                                                                             >  1:1273176/1‑148 (MQ=255)
                                   ttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                                                                             >  1:2981082/1‑148 (MQ=255)
                                      cTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGc                                                                                                          >  1:405581/1‑148 (MQ=255)
                                        aTGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCaa                                                                                                        >  1:2509782/1‑148 (MQ=255)
                                             gACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGAc                                                                                                   <  2:300555/148‑1 (MQ=255)
                                                cgcgCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                                                >  2:548873/1‑148 (MQ=255)
                                                  cgcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAc                                                                                              >  1:2737581/1‑148 (MQ=255)
                                                   gcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAcc                                                                                             <  1:3319951/148‑1 (MQ=255)
                                                     ccGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                                                                           >  2:2962597/1‑148 (MQ=255)
                                                     ccGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                                                                           <  2:1329377/148‑1 (MQ=255)
                                                      cGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgc                                                                                          >  2:1430660/1‑148 (MQ=255)
                                                       gCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgcg                                                                                         <  1:3091011/148‑1 (MQ=255)
                                                         tATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTg                                                                                       >  2:531885/1‑148 (MQ=255)
                                                            cgtcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                                                    <  1:411522/148‑1 (MQ=255)
                                                              tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCACAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTag                                                                                  <  2:2347068/148‑1 (MQ=255)
                                                              tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTag                                                                                  >  2:115168/1‑148 (MQ=255)
                                                                 cgcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcag                                                                               <  2:2830028/148‑1 (MQ=255)
                                                                  gcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                                                                              <  2:902132/148‑1 (MQ=255)
                                                                     ggtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGc                                                                           >  2:744214/1‑148 (MQ=255)
                                                                      gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                                                                          <  1:1032897/148‑1 (MQ=255)
                                                                        gggCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAtt                                                                        <  2:2217623/148‑1 (MQ=255)
                                                                        gggCGGCACGGTGGAAGAAGCTAGTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAtt                                                                        <  1:879040/148‑1 (MQ=255)
                                                                             gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                                                   >  1:2219193/1‑148 (MQ=255)
                                                                              cACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                                                  <  2:826009/148‑1 (MQ=255)
                                                                               aCGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCgg                                                                 <  1:3135144/148‑1 (MQ=255)
                                                                                cGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGt                                                                <  2:1885459/148‑1 (MQ=255)
                                                                                   tGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt                                                             <  1:2224164/148‑1 (MQ=255)
                                                                                     gaagaaGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCa                                                           >  1:922490/1‑148 (MQ=255)
                                                                                     gaagaaGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCa                                                           <  1:3060073/148‑1 (MQ=255)
                                                                                        gaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAaca                                                        >  2:492875/1‑148 (MQ=255)
                                                                                         aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                       >  2:2747945/1‑148 (MQ=255)
                                                                                         aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                       >  2:989432/1‑148 (MQ=255)
                                                                                          agCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACAcc                                                      >  2:1875526/1‑148 (MQ=255)
                                                                                            cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                                                    >  2:3375717/1‑148 (MQ=255)
                                                                                             tATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGcc                                                   >  1:1776954/1‑148 (MQ=255)
                                                                                               ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                                                 >  1:3268853/1‑148 (MQ=255)
                                                                                                tGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATc                                                >  1:983331/1‑148 (MQ=255)
                                                                                                 gATGTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCACAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGGAGCCGCGGCACTCTGCCGGCGAACAACACCGCCATCg                                               <  2:3349426/148‑1 (MQ=255)
                                                                                                  aTTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGc                                              >  2:577165/1‑148 (MQ=255)
                                                                                                  aTTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGc                                              <  2:2497991/148‑1 (MQ=255)
                                                                                                   tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGcc                                             <  2:2153908/148‑1 (MQ=255)
                                                                                                   tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGcc                                             <  2:2072797/148‑1 (MQ=255)
                                                                                                       tcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa                                         <  2:3270693/148‑1 (MQ=255)
                                                                                                        cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaa                                        >  1:1885279/1‑148 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                      >  1:2311371/1‑148 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                      >  1:3423343/1‑148 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                      >  2:1327363/1‑148 (MQ=255)
                                                                                                           gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGa                                     <  2:3344881/148‑1 (MQ=255)
                                                                                                            tCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGAt                                    >  2:1084642/1‑148 (MQ=255)
                                                                                                             cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                   >  1:2001267/1‑148 (MQ=255)
                                                                                                             cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                   <  1:2510544/148‑1 (MQ=255)
                                                                                                              aGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGc                                  >  1:3378411/1‑148 (MQ=255)
                                                                                                              aGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGc                                  <  1:3365131/148‑1 (MQ=255)
                                                                                                              aGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGc                                  <  2:757189/148‑1 (MQ=255)
                                                                                                               gATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCg                                 >  2:727875/1‑148 (MQ=255)
                                                                                                                 tGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCGGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg                               <  2:249479/148‑1 (MQ=255)
                                                                                                                    gcgAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGt                            <  2:2472984/148‑1 (MQ=255)
                                                                                                                    gcgAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGt                            <  2:1794622/148‑1 (MQ=255)
                                                                                                                       aaaTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCgg                         >  1:2204777/1‑148 (MQ=255)
                                                                                                                         aTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGAGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCTCCATCGCCCATAAAGATGCTGGCGTCGGCa                       >  1:1073874/1‑148 (MQ=255)
                                                                                                                           cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATg                     <  1:270105/148‑1 (MQ=255)
                                                                                                                            aCCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTTATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGa                    <  1:2649940/148‑1 (MQ=255)
                                                                                                                            aCCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGa                    <  2:1582705/148‑1 (MQ=255)
                                                                                                                             cccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAt                   >  2:3436282/1‑148 (MQ=255)
                                                                                                                              ccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt                  >  2:1655391/1‑148 (MQ=255)
                                                                                                                                tCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTgg                >  1:936161/1‑148 (MQ=255)
                                                                                                                                tCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTgg                >  1:1374598/1‑148 (MQ=255)
                                                                                                                                    tGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGcgcg            >  1:109135/1‑148 (MQ=255)
                                                                                                                                      aaaaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCggg          >  2:2093373/1‑148 (MQ=255)
                                                                                                                                        aaCCCCAAAGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa        <  1:3211156/148‑1 (MQ=255)
                                                                                                                                         aCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAt       >  1:2026365/1‑148 (MQ=255)
                                                                                                                                         aCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAt       >  1:3292151/1‑148 (MQ=255)
                                                                                                                                          ccccAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAtt      >  1:2667051/1‑148 (MQ=255)
                                                                                                                                          ccccAAAGTCACCAGTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAtt      <  2:61031/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGC  >  CP009273/335103‑335392

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: