Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,956,165 C→A 45.1% E475D (GAG→GAT gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,956,1650CA45.1% 199.4 / inf 253E475D (GAG→GATgppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base C (81/58);  new base A (57/57);  total (138/115)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 258 total aligned reads displayed.

GCGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAG  >  CP009273/3956017‑3956311
                                                                                                                                                    |                                                                                                                                                  
    aTCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCttt                                                                                                                                               <  1:486210/150‑1 (MQ=255)
      cTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCACTAATCTCTTTAc                                                                                                                                             >  1:2125945/1‑150 (MQ=255)
        gAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTAccc                                                                                                                                           >  2:1993820/1‑150 (MQ=255)
        gAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTAccc                                                                                                                                           <  2:4127158/150‑1 (MQ=255)
          aaGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAg                                                                                                                                         >  1:1448409/1‑150 (MQ=255)
          aaGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAg                                                                                                                                         <  2:1736635/150‑1 (MQ=255)
            gCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCg                                                                                                                                       <  1:3632874/150‑1 (MQ=255)
             cGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCgg                                                                                                                                      >  1:1949264/1‑150 (MQ=255)
               atgatgGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGAt                                                                                                                                    >  2:2255683/1‑150 (MQ=255)
               atgatgGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGAt                                                                                                                                    >  1:2986669/1‑150 (MQ=255)
                tgatgGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATg                                                                                                                                   <  2:3900302/150‑1 (MQ=255)
                  atgGCGGCAAAACGAGGGAAATAATCAAGGCCGTAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCAATTCCAGCGGCCAGTGGACATATCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGtt                                                                                                                                 <  1:2507809/150‑1 (MQ=255)
                         cAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCCGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAg                                                                                                                          >  2:3870813/1‑150 (MQ=255)
                         cAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAg                                                                                                                          <  2:885596/150‑1 (MQ=255)
                             aCGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCaa                                                                                                                      >  1:124921/1‑150 (MQ=255)
                               gAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAAcc                                                                                                                    >  2:2485023/1‑150 (MQ=255)
                                aGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCt                                                                                                                   <  1:1723313/150‑1 (MQ=255)
                                aGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGGTGGGTTAGCCAACCt                                                                                                                   >  2:1943059/1‑150 (MQ=255)
                                      ataatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGa                                                                                                             >  2:678086/1‑150 (MQ=255)
                                        aatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAg                                                                                                           >  1:4108939/1‑150 (MQ=255)
                                        aatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAg                                                                                                           <  1:1357606/150‑1 (MQ=255)
                                        aatCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAg                                                                                                           >  1:38361/1‑150 (MQ=255)
                                        aatAAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTACTGTGCAATAAAAACTTAACCCAGCGGATGTTGGGTTAGCCAACCTTCCGGacg                                                                                                           >  2:2606309/1‑148 (MQ=255)
                                         atCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGc                                                                                                          <  1:1029705/150‑1 (MQ=255)
                                         atCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACGTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCATTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGc                                                                                                          <  2:1086902/150‑1 (MQ=255)
                                            aaGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCACCCTTGCGGAAGCGTc                                                                                                       <  1:2630765/150‑1 (MQ=255)
                                            aaGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTc                                                                                                       <  1:3260114/150‑1 (MQ=255)
                                             aGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCa                                                                                                      >  2:1299612/1‑150 (MQ=255)
                                               gCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAg                                                                                                    <  2:2284127/150‑1 (MQ=255)
                                                  ggAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTc                                                                                                 <  1:422263/150‑1 (MQ=255)
                                                  ggAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTc                                                                                                 <  2:911746/150‑1 (MQ=255)
                                                       aGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAGGCTTCAAGTTCAACAg                                                                                            >  1:2531959/1‑150 (MQ=255)
                                                         aTGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGAGGAAGCGTCAAGGTCAACAGtt                                                                                          <  2:1525719/150‑1 (MQ=255)
                                                          tGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTc                                                                                         <  2:3719030/150‑1 (MQ=255)
                                                           gCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCa                                                                                        <  2:3700311/150‑1 (MQ=255)
                                                            cGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCAt                                                                                       <  2:3314081/150‑1 (MQ=255)
                                                            cGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCAt                                                                                       <  1:1992314/150‑1 (MQ=255)
                                                            cGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCAt                                                                                       <  1:4040353/150‑1 (MQ=255)
                                                             gTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATg                                                                                      <  1:269303/150‑1 (MQ=255)
                                                             gTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATg                                                                                      >  1:3813509/1‑150 (MQ=255)
                                                                aGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGtt                                                                                   >  1:2999358/1‑150 (MQ=255)
                                                                aGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGtt                                                                                   >  1:669428/1‑150 (MQ=255)
                                                                aGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGtt                                                                                   <  1:106160/150‑1 (MQ=255)
                                                                aGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTGGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGtt                                                                                   >  2:354156/1‑150 (MQ=255)
                                                                   atcgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGc                                                                                >  1:4017336/1‑150 (MQ=255)
                                                                   atcgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGc                                                                                >  1:3445291/1‑150 (MQ=255)
                                                                   atcgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGc                                                                                >  1:4172219/1‑150 (MQ=255)
                                                                   atcgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGc                                                                                <  1:1612447/150‑1 (MQ=255)
                                                                   atcgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGc                                                                                >  2:11516/1‑150 (MQ=255)
                                                                     cgcatcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCt                                                                              <  1:152263/150‑1 (MQ=255)
                                                                       catcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGt                                                                            >  2:3532540/1‑150 (MQ=255)
                                                                        atcCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTa                                                                           <  1:533002/150‑1 (MQ=255)
                                                                           cGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATg                                                                        >  1:3784965/1‑150 (MQ=255)
                                                                           cGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATg                                                                        >  1:2753712/1‑150 (MQ=255)
                                                                           cGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATg                                                                        >  2:1058213/1‑150 (MQ=255)
                                                                              cACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCa                                                                     >  2:1200680/1‑150 (MQ=255)
                                                                               aCTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAt                                                                    >  2:3115884/1‑150 (MQ=255)
                                                                               aCTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAt                                                                    <  2:322355/150‑1 (MQ=255)
                                                                                 ttACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAtct                                                                  >  1:20659/1‑150 (MQ=255)
                                                                                 ttACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCAtct                                                                  <  2:1351230/150‑1 (MQ=255)
                                                                                    cTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTg                                                               >  2:496639/1‑150 (MQ=255)
                                                                                      cATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGc                                                             <  1:1258881/150‑1 (MQ=255)
                                                                                       aTTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCa                                                            <  2:829311/150‑1 (MQ=255)
                                                                                         tAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACg                                                          >  1:1078573/1‑150 (MQ=255)
                                                                                            tGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACgaga                                                       >  1:2636359/1‑150 (MQ=255)
                                                                                              cACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATc                                                     >  1:1862195/1‑150 (MQ=255)
                                                                                              cACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATc                                                     <  1:2369310/150‑1 (MQ=255)
                                                                                                 ttCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTc                                                  >  1:521554/1‑150 (MQ=255)
                                                                                                  tCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCa                                                 >  1:1175632/1‑150 (MQ=255)
                                                                                                   ccAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCAc                                                >  2:4171397/1‑150 (MQ=255)
                                                                                                   ccAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCAc                                                <  2:3574056/150‑1 (MQ=255)
                                                                                                           cAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACgg                                        >  1:1421725/1‑150 (MQ=255)
                                                                                                            aGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGc                                       >  2:2963453/1‑150 (MQ=255)
                                                                                                            aGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGc                                       <  2:2522982/150‑1 (MQ=255)
                                                                                                               ggACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTgg                                    >  1:3749913/1‑150 (MQ=255)
                                                                                                                gACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCAAGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGc                                   >  1:869825/1‑150 (MQ=255)
                                                                                                                  cATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCaa                                 <  2:3190795/150‑1 (MQ=255)
                                                                                                                  cATAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCaa                                 <  1:1769028/150‑1 (MQ=255)
                                                                                                                   aTAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCaaa                                >  1:3839053/1‑150 (MQ=255)
                                                                                                                    tAGCTCTGCCACTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCaaaa                               >  1:738210/1‑150 (MQ=255)
                                                                                                                      gCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCAAGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAat                             >  1:745952/1‑150 (MQ=255)
                                                                                                                          tGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAatgatg                         <  1:3754571/150‑1 (MQ=255)
                                                                                                                              aCTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCa                     >  2:2931162/1‑150 (MQ=255)
                                                                                                                              aCTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCa                     >  1:2190367/1‑150 (MQ=255)
                                                                                                                              aCTGGCTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCa                     >  1:319298/1‑150 (MQ=255)
                                                                                                                                 ggCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGc                  <  1:3187628/150‑1 (MQ=255)
                                                                                                                                  gCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCg                 >  2:2201583/1‑150 (MQ=255)
                                                                                                                                  gCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCg                 >  2:2056414/1‑150 (MQ=255)
                                                                                                                                   cTTTCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgt                <  2:65588/150‑1 (MQ=255)
                                                                                                                                     ttCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgtag              >  2:3744574/1‑150 (MQ=255)
                                                                                                                                     ttCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgtag              >  2:1774651/1‑150 (MQ=255)
                                                                                                                                      tCCTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCgtagt             <  2:2778067/150‑1 (MQ=255)
                                                                                                                                        cTGAGCAATAATATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTaa           <  2:1949264/150‑1 (MQ=255)
                                                                                                                                         tGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTaaa          >  1:1984439/1‑150 (MQ=255)
                                                                                                                                             cAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCCCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGc      >  2:3982763/1‑150 (MQ=255)
                                                                                                                                             cAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGc      >  1:2825056/1‑150 (MQ=255)
                                                                                                                                             cAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGc      <  1:3543935/150‑1 (MQ=255)
                                                                                                                                             cAATAATCTCTTTACCCAGCGGATGTTGGGTGAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGc      <  2:1288941/150‑1 (MQ=255)
                                                                                                                                              aataatATCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCa     <  2:3409368/150‑1 (MQ=255)
                                                                                                                                                taatCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCaga   >  2:1551173/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGAATCTGAAAGCGATGATGGCGGCAAAACGAGGGAAATAATCAAGGCCGGAAAAGATGCGTCAGCATCGCATCCGGCACTTACTCATTAATGCACTTCCAGCGGCCAGTGGACATAGCTCTGCCACTGGCTTTCCTGAGCAATAATCTCTTTACCCAGCGGATGTTGGGTTAGCCAACCTTGCGGAAGCGTCAAGGTCAACAGTTCATGGTTAGCCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAG  >  CP009273/3956017‑3956311

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: