Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,285 C→T 6.7% A133A (GCG→GCA gatY ← D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,170,2850CT 6.7% 1074.0 / 42.8 299A133A (GCG→GCAgatYD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (132/147);  new base T (9/11);  total (141/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

ACATCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCG  >  CP009273/2170137‑2170430
                                                                                                                                                    |                                                                                                                                                 
  aTCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGc                                                                                                                                                 >  2:2171835/1‑149 (MQ=255)
   tCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGcc                                                                                                                                                >  2:5118400/1‑149 (MQ=255)
   tCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGcc                                                                                                                                                <  2:4536798/149‑1 (MQ=255)
   tCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGcc                                                                                                                                                >  2:2516392/1‑149 (MQ=255)
    ccccTGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCt                                                                                                                                               <  2:3163079/149‑1 (MQ=255)
       aTGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGa                                                                                                                                            <  1:4143109/149‑1 (MQ=255)
        tGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGAc                                                                                                                                           >  1:958394/1‑149 (MQ=255)
             cGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGa                                                                                                                                      <  2:911164/149‑1 (MQ=255)
              gTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGAc                                                                                                                                     <  1:2285747/149‑1 (MQ=255)
               tGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACa                                                                                                                                    >  2:670044/1‑149 (MQ=255)
                   gATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCaa                                                                                                                                >  2:1966262/1‑149 (MQ=255)
                   gATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCaa                                                                                                                                <  2:3956699/149‑1 (MQ=255)
                     tGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAg                                                                                                                              >  1:1184267/1‑149 (MQ=255)
                     tGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAg                                                                                                                              <  2:2743204/149‑1 (MQ=255)
                      ggCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGc                                                                                                                             <  1:4805692/149‑1 (MQ=255)
                         gACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGAt                                                                                                                          >  2:2650837/1‑149 (MQ=255)
                         gACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGAt                                                                                                                          >  1:1228245/1‑149 (MQ=255)
                           ccGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATgg                                                                                                                        <  2:4863676/149‑1 (MQ=255)
                              ccAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaa                                                                                                                     <  2:2033599/149‑1 (MQ=255)
                               cAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaaa                                                                                                                    >  1:1949943/1‑149 (MQ=255)
                                  ggAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAAt                                                                                                                 <  2:1068030/149‑1 (MQ=255)
                                    aaTCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATcc                                                                                                               >  1:2453902/1‑149 (MQ=255)
                                    aaTCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATcc                                                                                                               >  2:2414412/1‑149 (MQ=255)
                                     aTCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATcca                                                                                                              <  2:2361122/149‑1 (MQ=255)
                                         aTTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGCACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCcaccac                                                                                                          >  2:2938899/1‑149 (MQ=255)
                                            ccGGGTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACctc                                                                                                       <  2:3395317/149‑1 (MQ=255)
                                              ggTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCtt                                                                                                     <  2:3985751/149‑1 (MQ=255)
                                                ttGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTg                                                                                                   <  1:1530190/149‑1 (MQ=255)
                                                ttGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTg                                                                                                   >  2:4610409/1‑149 (MQ=255)
                                                ttGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTg                                                                                                   >  2:3031757/1‑149 (MQ=255)
                                                  gCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGAc                                                                                                 >  1:3617783/1‑149 (MQ=255)
                                                  gCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGAc                                                                                                 >  2:1744815/1‑149 (MQ=255)
                                                    cTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGAccc                                                                                               <  2:4052092/149‑1 (MQ=255)
                                                    cTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGAccc                                                                                               <  2:136594/149‑1 (MQ=255)
                                                        gCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGa                                                                                           <  2:4382989/149‑1 (MQ=255)
                                                          aaaTTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGaaa                                                                                         >  1:248946/1‑149 (MQ=255)
                                                             ttCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAatat                                                                                      <  2:2862510/149‑1 (MQ=255)
                                                                    gcCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGAGCCGTGAAATATTTTgcgc                                                                               <  1:1521966/149‑1 (MQ=255)
                                                                     ccTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCa                                                                              >  1:2863923/1‑149 (MQ=255)
                                                                     ccTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCa                                                                              >  2:1637442/1‑149 (MQ=255)
                                                                     ccTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCa                                                                              >  2:4196671/1‑149 (MQ=255)
                                                                        gAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCaaaa                                                                           <  1:5070683/149‑1 (MQ=255)
                                                                         aGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAg                                                                          <  2:85455/149‑1 (MQ=255)
                                                                           cGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGc                                                                        <  2:1169228/149‑1 (MQ=255)
                                                                              ggTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCaaa                                                                     >  2:3685152/1‑149 (MQ=255)
                                                                              ggTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCaaa                                                                     <  2:3614065/149‑1 (MQ=255)
                                                                                 tGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGc                                                                  >  2:166720/1‑149 (MQ=255)
                                                                                    tgtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAg                                                               >  1:1773702/1‑149 (MQ=255)
                                                                                       aCAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCg                                                            <  1:1612817/149‑1 (MQ=255)
                                                                                       aCAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCg                                                            >  1:293662/1‑149 (MQ=255)
                                                                                        cAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGt                                                           >  2:4334613/1‑149 (MQ=255)
                                                                                          aCGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCa                                                         <  1:3141566/149‑1 (MQ=255)
                                                                                             cATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATc                                                      <  2:3119860/149‑1 (MQ=255)
                                                                                             cATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATc                                                      <  2:3457335/149‑1 (MQ=255)
                                                                                               tCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCAt                                                    >  1:733763/1‑149 (MQ=255)
                                                                                                cGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATg                                                   <  2:1294232/149‑1 (MQ=255)
                                                                                                cGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATCGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATg                                                   <  2:3889887/149‑1 (MQ=255)
                                                                                                 ggCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGa                                                  <  2:1694366/149‑1 (MQ=255)
                                                                                                 ggCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGa                                                  >  1:563907/1‑149 (MQ=255)
                                                                                                   cTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACt                                                >  1:928450/1‑149 (MQ=255)
                                                                                                     tCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGa                                              <  2:4079298/149‑1 (MQ=255)
                                                                                                      cATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAg                                             >  2:2404411/1‑149 (MQ=255)
                                                                                                      cATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAg                                             <  2:5177715/149‑1 (MQ=255)
                                                                                                       aTTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAgc                                            <  1:4981341/149‑1 (MQ=255)
                                                                                                        ttGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAgcg                                           >  1:620845/1‑149 (MQ=255)
                                                                                                              tGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGcc                                     >  1:3946823/1‑149 (MQ=255)
                                                                                                              tGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGcc                                     <  1:4089548/149‑1 (MQ=255)
                                                                                                                cacaTCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAg                                   <  1:2708761/149‑1 (MQ=255)
                                                                                                                  catcatCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaa                                 <  1:2175584/149‑1 (MQ=255)
                                                                                                                  catcatCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaa                                 <  1:2666066/149‑1 (MQ=255)
                                                                                                                     catcTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaacga                              >  1:4236056/1‑149 (MQ=255)
                                                                                                                      atcTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaacgaa                             >  2:1766853/1‑149 (MQ=255)
                                                                                                                        cTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAAcc                           <  1:4161082/149‑1 (MQ=255)
                                                                                                                          tCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCtt                         <  1:4938839/149‑1 (MQ=255)
                                                                                                                            cTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCt                       <  2:3865178/149‑1 (MQ=255)
                                                                                                                            cTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCt                       <  2:919997/149‑1 (MQ=255)
                                                                                                                             tGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTg                      <  2:4973593/149‑1 (MQ=255)
                                                                                                                             tGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTg                      <  2:1436387/149‑1 (MQ=255)
                                                                                                                              ggCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGa                     <  1:1273319/149‑1 (MQ=255)
                                                                                                                                 cgccAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTACGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCg                  >  1:3904150/1‑149 (MQ=255)
                                                                                                                                  gccAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGa                 >  2:3505789/1‑149 (MQ=255)
                                                                                                                                  gccAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGa                 <  1:3573520/149‑1 (MQ=255)
                                                                                                                                   ccAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGat                <  2:3243038/149‑1 (MQ=255)
                                                                                                                                   ccAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGat                >  2:4517685/1‑149 (MQ=255)
                                                                                                                                    cAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGata               >  2:4764317/1‑149 (MQ=255)
                                                                                                                                     aaGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGatat              <  2:1949903/149‑1 (MQ=255)
                                                                                                                                     aaGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGatat              <  1:1431904/149‑1 (MQ=255)
                                                                                                                                        ttGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGt           >  2:1515937/1‑149 (MQ=255)
                                                                                                                                         tGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTc          <  1:2777833/149‑1 (MQ=255)
                                                                                                                                         tGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTc          <  1:4459953/149‑1 (MQ=255)
                                                                                                                                          gCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCa         <  1:3707442/149‑1 (MQ=255)
                                                                                                                                          gCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCa         <  1:3529586/149‑1 (MQ=255)
                                                                                                                                          gCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCa         <  1:2649578/149‑1 (MQ=255)
                                                                                                                                          gCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCa         <  1:46814/149‑1 (MQ=255)
                                                                                                                                           ccccAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCaa        >  1:13433/1‑149 (MQ=255)
                                                                                                                                           ccccAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCaa        >  1:13258/1‑149 (MQ=255)
                                                                                                                                            cccAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCaaa       >  1:1245324/1‑149 (MQ=255)
                                                                                                                                             ccAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAAt      >  1:1913885/1‑149 (MQ=255)
                                                                                                                                             ccAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAAt      >  1:3261097/1‑149 (MQ=255)
                                                                                                                                                 cTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCg  <  1:5149128/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACATCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCG  >  CP009273/2170137‑2170430

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: