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,2470CA98.8% 1002.0 / 0.5 338N404K (AAC→AAAyahGDUF1116 family protein
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base A (190/140);  total (193/145)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 338 total aligned reads displayed.

CGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCA  >  CP009273/335100‑335393
                                                                                                                                                   |                                                                                                                                                  
cGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGTTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGt                                                                                                                                                  <  2:1286040/150‑1 (MQ=255)
cGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGt                                                                                                                                                  >  1:4768575/1‑150 (MQ=255)
   cAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTCTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAc                                                                                                                                               >  2:281944/1‑150 (MQ=255)
   cAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAc                                                                                                                                               >  2:673812/1‑150 (MQ=255)
   cAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAc                                                                                                                                               <  1:2232779/150‑1 (MQ=255)
   cAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAc                                                                                                                                               <  1:2138617/150‑1 (MQ=255)
    aGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCAcc                                                                                                                                              >  1:1179145/1‑150 (MQ=255)
     gcgcCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCa                                                                                                                                             <  1:5077417/150‑1 (MQ=255)
       gcCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCAtt                                                                                                                                           >  1:2007403/1‑150 (MQ=255)
        ccATCACCGAAACCGACGGTCTTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGGAGAAGCTATTGATTTCTCCCGTCAGATGCCCGAAATCACCCTCTCTGAAAACCCCAACGTCACCAtac                                                                                                                                          >  2:2633830/1‑148 (MQ=255)
         cATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTcc                                                                                                                                         <  1:29587/150‑1 (MQ=255)
          aTCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCg                                                                                                                                        >  1:4605705/1‑150 (MQ=255)
            cACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCgct                                                                                                                                      <  2:1829454/150‑1 (MQ=255)
                    ccTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtt                                                                                                                              >  1:3516407/1‑150 (MQ=255)
                      tACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGtttt                                                                                                                            <  2:3099383/150‑1 (MQ=255)
                        cGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTAt                                                                                                                          <  1:383329/150‑1 (MQ=255)
                          gTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATgg                                                                                                                        >  2:2156741/1‑150 (MQ=255)
                          gTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATgg                                                                                                                        >  2:3015444/1‑150 (MQ=255)
                             atGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCg                                                                                                                     >  2:1716368/2‑150 (MQ=255)
                                  ggATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                                >  2:2421473/1‑150 (MQ=255)
                                  ggATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                                                                                <  1:4117551/150‑1 (MQ=255)
                                   gATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGt                                                                                                               >  2:4340123/1‑150 (MQ=255)
                                    aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                                                                              >  1:947083/1‑150 (MQ=255)
                                    aTTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                                                                              >  1:761426/1‑150 (MQ=255)
                                         cTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCaa                                                                                                         <  1:5227547/150‑1 (MQ=255)
                                           aTGGCGACCGCGCCCGCTCTCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATc                                                                                                       <  1:2352492/150‑1 (MQ=255)
                                              gCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGAc                                                                                                    <  1:2435933/150‑1 (MQ=255)
                                                gACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACAt                                                                                                  >  2:4239704/1‑150 (MQ=255)
                                                 aCCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                                                 >  2:1943858/1‑150 (MQ=255)
                                                    gcgcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGAGTGCCGTCGGCAATCGACATCAcc                                                                                              >  2:2445712/1‑150 (MQ=255)
                                                     cgcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAccc                                                                                             <  1:149717/150‑1 (MQ=255)
                                                      gcCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCGCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                                                                            >  2:443803/1‑150 (MQ=255)
                                                       cccGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgc                                                                                           <  2:1757877/150‑1 (MQ=255)
                                                         cGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                                                         >  1:138219/1‑150 (MQ=255)
                                                           cTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTGTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTgg                                                                                       >  1:1701864/1‑150 (MQ=255)
                                                            tATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                                                      >  2:3417053/1‑150 (MQ=255)
                                                             aTCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                                                     >  1:3283521/1‑150 (MQ=255)
                                                               cgtcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTag                                                                                   >  2:1874178/1‑150 (MQ=255)
                                                                 tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagca                                                                                 <  1:2551986/150‑1 (MQ=255)
                                                                  cgcgCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcag                                                                                <  2:1598232/150‑1 (MQ=255)
                                                                  cgcgCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcag                                                                                >  1:4309088/1‑150 (MQ=255)
                                                                   gcgcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                                                                               >  2:1498480/1‑150 (MQ=255)
                                                                   gcgcTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                                                                               >  2:3950215/1‑150 (MQ=255)
                                                                      cTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGc                                                                            >  1:561532/1‑150 (MQ=255)
                                                                         gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAtt                                                                         <  2:1924418/150‑1 (MQ=255)
                                                                         gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAtt                                                                         >  1:1397253/1‑150 (MQ=255)
                                                                          tggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTc                                                                        <  1:512337/150‑1 (MQ=255)
                                                                          tggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTc                                                                        >  1:929648/1‑150 (MQ=255)
                                                                               ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                                                   <  2:766059/150‑1 (MQ=255)
                                                                               ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                                                   >  1:3045254/1‑150 (MQ=255)
                                                                                gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCgg                                                                  >  1:1264622/1‑150 (MQ=255)
                                                                                    ggTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt                                                              <  1:4998669/150‑1 (MQ=255)
                                                                                    ggTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt                                                              >  2:1439829/1‑150 (MQ=255)
                                                                                         aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAaca                                                         <  1:605653/150‑1 (MQ=255)
                                                                                          agaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                        >  2:4159117/1‑150 (MQ=255)
                                                                                          agaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                        <  1:3252067/150‑1 (MQ=255)
                                                                                          agaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                        >  2:3812992/1‑150 (MQ=255)
                                                                                          agaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAacac                                                        <  1:422969/150‑1 (MQ=255)
                                                                                            aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                                                      >  1:4469783/1‑150 (MQ=255)
                                                                                            aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                                                      >  2:2724969/1‑150 (MQ=255)
                                                                                            aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                                                      >  1:4750044/1‑150 (MQ=255)
                                                                                            aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                                                      <  2:4504680/150‑1 (MQ=255)
                                                                                              gCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGcc                                                    <  1:362087/150‑1 (MQ=255)
                                                                                                tATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCAt                                                  >  2:3912301/1‑150 (MQ=255)
                                                                                                 aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACCTCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATc                                                 >  2:849010/1‑150 (MQ=255)
                                                                                                 aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATc                                                 <  2:1371035/150‑1 (MQ=255)
                                                                                                  ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCg                                                <  2:2100209/150‑1 (MQ=255)
                                                                                                   tGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGc                                               >  1:778817/1‑150 (MQ=255)
                                                                                                     aTTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCAGCTCGGTATTATGGGCGTGCCGACAGCAAACTACATCACCCGCGTGGCTAGCAGCGGCATACTGCCGGAGATCAACACCGCCATCGccc                                             >  2:891650/1‑150 (MQ=255)
                                                                                                      tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa                                            <  2:3634938/150‑1 (MQ=255)
                                                                                                        tctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATa                                          <  2:5153010/150‑1 (MQ=255)
                                                                                                           cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg                                       <  1:2888749/150‑1 (MQ=255)
                                                                                                           cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGACGGTGATCAACACCGCCATCGCCCATAAAg                                       <  1:3392029/150‑1 (MQ=255)
                                                                                                            ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGa                                      >  2:1166161/1‑150 (MQ=255)
                                                                                                              gTCAGATGCGCGAAATCACCCTCGGTGAACACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                    <  2:3315550/150‑1 (MQ=255)
                                                                                                              gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg                                    >  2:1625622/1‑150 (MQ=255)
                                                                                                                 aGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCCTAAAGATGCgg                                 >  1:1590197/1‑150 (MQ=255)
                                                                                                                 aGATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCgg                                 >  2:2315918/1‑150 (MQ=255)
                                                                                                                  gATGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg                                >  2:3560993/1‑150 (MQ=255)
                                                                                                                   aTGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATTAAGATGCGGGc                               <  1:2735064/150‑1 (MQ=255)
                                                                                                                   aTGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGc                               >  1:3900368/1‑150 (MQ=255)
                                                                                                                    tGCGCGAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCg                              >  1:398923/1‑150 (MQ=255)
                                                                                                                      cgcgAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTc                            <  1:121354/150‑1 (MQ=255)
                                                                                                                        cgAAATCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCgg                          >  2:4873069/1‑150 (MQ=255)
                                                                                                                          aaaTCACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCa                        >  1:4699045/1‑150 (MQ=255)
                                                                                                                              cACCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAt                    <  2:3138945/150‑1 (MQ=255)
                                                                                                                              cACCCGCGGTGAAACCCCCCAAGACACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAGCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAt                    <  2:2004539/150‑1 (MQ=255)
                                                                                                                               aCCCTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGAtt                   <  2:584064/150‑1 (MQ=255)
                                                                                                                                cccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTg                  >  1:2670590/1‑150 (MQ=255)
                                                                                                                                cccTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTg                  <  1:3941177/150‑1 (MQ=255)
                                                                                                                                  cTCGGTGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGgc                <  1:3484992/150‑1 (MQ=255)
                                                                                                                                       tGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCggg           >  1:644826/1‑150 (MQ=255)
                                                                                                                                       tGAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCggg           >  2:5028488/1‑150 (MQ=255)
                                                                                                                                        gAAAACCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGc          >  1:1062811/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCa         >  2:4924372/1‑150 (MQ=255)
                                                                                                                                           aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCAtt       >  1:1054084/1‑150 (MQ=255)
                                                                                                                                           aaCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGACATGATTGGCGCGGGCAtt       >  1:795194/1‑150 (MQ=255)
                                                                                                                                            aCCCCAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGACGCGGGCGTCGGCATGATTGGCGCGGGCATtg      >  1:4837589/1‑150 (MQ=255)
                                                                                                                                               ccAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGc   >  2:5131257/1‑150 (MQ=255)
                                                                                                                                                cAAAGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCa  <  1:1817861/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGACAGCGCCATCACCGAAACCTACGGTATTGGCGGATTTGCTATGGCGACCGCGCCCGCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCA  >  CP009273/335100‑335393

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: