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,2470CA95.4% 532.5 / 15.6 187N404K (AAC→AAAyahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base C (7/2);  new base A (106/72);  total (113/74)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 187 total aligned reads displayed.

CAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTG  >  CP009273/335103‑335389
                                                                                                                                                |                                                                                                                                              
cAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTc                                                                                                                                             <  1:612095/148‑1 (MQ=255)
 aGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCa                                                                                                                                            <  1:3030499/148‑1 (MQ=255)
   cgcCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAcc                                                                                                                                          >  1:2943186/1‑148 (MQ=255)
      cATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAtt                                                                                                                                       >  1:2376301/1‑148 (MQ=255)
       aTCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTc                                                                                                                                      >  1:1715889/1‑148 (MQ=255)
         cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                    >  2:1949637/1‑148 (MQ=255)
         cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                    <  1:2239777/148‑1 (MQ=255)
          aCCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAACTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgc                                                                                                                                   >  2:2735503/1‑148 (MQ=255)
            cGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgctg                                                                                                                                 >  2:2467548/1‑148 (MQ=255)
            cGAAACCTACCTTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGACGCGCTGGAGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCAGTCAGATGCGCGAAATCACCCTCGGTGAAAACCACAAAGACACCATTACgctg                                                                                                                                 >  2:2467816/1‑148 (MQ=255)
                 ccTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCgg                                                                                                                            >  1:1000711/1‑148 (MQ=255)
                   tACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtt                                                                                                                          <  1:2849332/148‑1 (MQ=255)
                       gTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTAt                                                                                                                      <  2:1392046/148‑1 (MQ=255)
                       gTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTAt                                                                                                                      <  1:2026538/148‑1 (MQ=255)
                        tATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATg                                                                                                                     >  1:3412454/1‑148 (MQ=255)
                         aTTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTCCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATgg                                                                                                                    >  2:1041107/1‑148 (MQ=255)
                          ttGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATggg                                                                                                                   >  1:289950/1‑148 (MQ=255)
                           tGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGc                                                                                                                  >  1:3238865/1‑148 (MQ=255)
                           tGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGc                                                                                                                  >  1:2126540/1‑148 (MQ=255)
                            ggCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCg                                                                                                                 >  2:1265507/1‑148 (MQ=255)
                             gCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGt                                                                                                                >  2:863730/1‑148 (MQ=255)
                                gATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGcc                                                                                                             >  1:1654256/1‑148 (MQ=255)
                                gATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGcc                                                                                                             >  1:210844/1‑148 (MQ=255)
                                 agttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                            >  1:220739/3‑148 (MQ=255)
                                 aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                            >  1:220570/1‑148 (MQ=255)
                                 aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                            >  1:220589/1‑148 (MQ=255)
                                  tttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGt                                                                                                           <  2:341943/148‑1 (MQ=255)
                                  tttGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGt                                                                                                           <  2:344525/148‑1 (MQ=255)
                                    tGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCg                                                                                                         <  1:782574/148‑1 (MQ=255)
                                         tGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAAt                                                                                                    <  1:1140634/148‑1 (MQ=255)
                                              aCCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACa                                                                                               <  1:2430506/148‑1 (MQ=255)
                                               ccGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACAt                                                                                              <  1:295352/148‑1 (MQ=255)
                                               ccGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACAt                                                                                              >  2:2505191/1‑148 (MQ=255)
                                                cgcgCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                                             >  2:3395624/1‑148 (MQ=255)
                                                cgcgCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGAGTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                                             <  2:799531/148‑1 (MQ=255)
                                                      cGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgc                                                                                       >  2:2289865/1‑148 (MQ=255)
                                                        cTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                                                     <  2:2424210/148‑1 (MQ=255)
                                                          aTCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTgg                                                                                   >  1:3122641/1‑148 (MQ=255)
                                                          aTCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTgg                                                                                   >  1:3185271/1‑148 (MQ=255)
                                                              tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTag                                                                               >  2:1247361/1‑148 (MQ=255)
                                                               cgcgCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagc                                                                              >  1:2769080/1‑148 (MQ=255)
                                                                 cgcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcag                                                                            <  2:741518/148‑1 (MQ=255)
                                                                 cgcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcag                                                                            <  1:2386935/148‑1 (MQ=255)
                                                                               aCGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCgg                                                              <  2:3524599/148‑1 (MQ=255)
                                                                                cGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGt                                                             >  1:716217/1‑148 (MQ=255)
                                                                                   tGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt                                                          <  1:547972/148‑1 (MQ=255)
                                                                                    ggAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                                                         >  2:2029830/1‑148 (MQ=255)
                                                                                    ggAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                                                         >  2:2029836/1‑148 (MQ=255)
                                                                                     gaagaaGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCa                                                        <  1:3518565/148‑1 (MQ=255)
                                                                                      aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa                                                       <  1:1349958/148‑1 (MQ=255)
                                                                                       agacgCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAac                                                      >  2:3413549/1‑148 (MQ=255)
                                                                                         aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                    >  2:1605960/1‑148 (MQ=255)
                                                                                         aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                    <  1:928742/148‑1 (MQ=255)
                                                                                           gCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                                                  <  2:1584195/148‑1 (MQ=255)
                                                                                            cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                                                 <  2:3535695/148‑1 (MQ=255)
                                                                                               ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                                              <  2:3331932/148‑1 (MQ=255)
                                                                                               ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                                              <  2:1320473/148‑1 (MQ=255)
                                                                                               ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                                              <  2:521482/148‑1 (MQ=255)
                                                                                                  aTTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGc                                           >  1:1376991/1‑148 (MQ=255)
                                                                                                   tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGcc                                          >  1:1910744/1‑148 (MQ=255)
                                                                                                   tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGcc                                          <  2:1833825/148‑1 (MQ=255)
                                                                                                    ttCTCCCGTCAGATGCGCGAACTCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCTccc                                         <  2:2595747/148‑1 (MQ=255)
                                                                                                    ttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGccc                                         >  1:2086623/1‑148 (MQ=255)
                                                                                                    ttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGccc                                         <  2:2823566/148‑1 (MQ=255)
                                                                                                    ttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGccc                                         <  2:195244/148‑1 (MQ=255)
                                                                                                       tcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa                                      >  1:2156256/1‑148 (MQ=255)
                                                                                                       tcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa                                      >  1:367952/1‑148 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAACACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                   <  1:1889706/148‑1 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                   <  1:1117905/148‑1 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                   <  1:2616618/148‑1 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                   <  1:2617036/148‑1 (MQ=255)
                                                                                                          cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTACGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                   <  1:2836071/148‑1 (MQ=255)
                                                                                                           gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGa                                  <  2:2126492/148‑1 (MQ=255)
                                                                                                               gATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCg                              >  1:1851561/1‑148 (MQ=255)
                                                                                                                aTGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCgg                             >  1:1221327/1‑148 (MQ=255)
                                                                                                                aTGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCgg                             >  1:1678505/1‑148 (MQ=255)
                                                                                                                 tGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg                            >  2:1706869/1‑148 (MQ=255)
                                                                                                                  gcgcGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGc                           >  2:1406621/1‑148 (MQ=255)
                                                                                                                    gcgAAATCACCCTCGGTGAAAACCCCACAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGt                         <  1:1434134/148‑1 (MQ=255)
                                                                                                                    gcgAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTTCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGt                         <  2:147215/148‑1 (MQ=255)
                                                                                                                    gcgAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGt                         <  2:814206/148‑1 (MQ=255)
                                                                                                                       aaaTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCgg                      >  2:3293816/1‑148 (MQ=255)
                                                                                                                       aaaTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCgg                      >  2:3294134/1‑148 (MQ=255)
                                                                                                                         aTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCa                    <  1:451558/148‑1 (MQ=255)
                                                                                                                         aTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCa                    >  2:3241967/1‑148 (MQ=255)
                                                                                                                             cccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAt                >  2:346468/1‑148 (MQ=255)
                                                                                                                             cccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAt                >  1:1490216/1‑148 (MQ=255)
                                                                                                                              ccTCTGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt               >  1:1601148/1‑148 (MQ=255)
                                                                                                                              ccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt               <  1:1203094/148‑1 (MQ=255)
                                                                                                                              ccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt               <  1:990372/148‑1 (MQ=255)
                                                                                                                              ccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt               >  1:1317662/1‑148 (MQ=255)
                                                                                                                               cTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTg              <  1:1041107/148‑1 (MQ=255)
                                                                                                                                 cGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGgc            >  1:843648/1‑148 (MQ=255)
                                                                                                                                        aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa     >  2:2557874/1‑148 (MQ=255)
                                                                                                                                        aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa     >  1:660658/1‑148 (MQ=255)
                                                                                                                                        aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCAGGATTGGCGCGGGCa     >  2:2557958/1‑148 (MQ=255)
                                                                                                                                         aCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAt    <  2:2133543/148‑1 (MQ=255)
                                                                                                                                          ccccAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAtt   >  1:3146511/1‑148 (MQ=255)
                                                                                                                                          ccccAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAtt   >  1:3146267/1‑148 (MQ=255)
                                                                                                                                           cccAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATtg  <  1:2910905/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTG  >  CP009273/335103‑335389

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: