Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,285 C→T 6.4% 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.4% 1063.4 / 38.6 298A133A (GCG→GCA) gatYD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (128/150);  new base T (9/10);  total (138/160)
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 301 total aligned reads displayed.

ACATCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGT  >  CP009273/2170137‑2170433
                                                                                                                                                    |                                                                                                                                                    
    cccATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCt                                                                                                                                                  >  2:3553661/1‑149 (MQ=255)
       aTGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGa                                                                                                                                               >  2:3589244/1‑149 (MQ=255)
       aTGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGa                                                                                                                                               <  1:1484407/149‑1 (MQ=255)
        tGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGAc                                                                                                                                              >  2:3414557/1‑149 (MQ=255)
                gCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACAt                                                                                                                                      <  2:2584688/149‑1 (MQ=255)
                 ccGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATc                                                                                                                                     <  1:1816920/149‑1 (MQ=255)
                     tGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAg                                                                                                                                 <  1:813498/149‑1 (MQ=255)
                     tGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAg                                                                                                                                 >  1:1388231/1‑149 (MQ=255)
                          aCCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATg                                                                                                                            <  2:2816106/149‑1 (MQ=255)
                             gCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCa                                                                                                                         <  2:4586875/149‑1 (MQ=255)
                              ccAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaa                                                                                                                        <  2:1598992/149‑1 (MQ=255)
                              ccAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaa                                                                                                                        <  2:124502/149‑1 (MQ=255)
                                 gggAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCaaaaa                                                                                                                     <  1:51786/149‑1 (MQ=255)
                                 gggAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCaaaaa                                                                                                                     <  1:2392762/149‑1 (MQ=255)
                                 gggAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaaaaa                                                                                                                     <  2:1450880/149‑1 (MQ=255)
                                 gggAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaaaaa                                                                                                                     <  1:2094237/149‑1 (MQ=255)
                                 gggAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaaaaa                                                                                                                     <  1:2144099/149‑1 (MQ=255)
                                 gggAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCaaaaa                                                                                                                     >  1:2630521/1‑149 (MQ=255)
                                    aaTCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATcc                                                                                                                  <  2:1825030/149‑1 (MQ=255)
                                     aTCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATcca                                                                                                                 <  2:1942762/149‑1 (MQ=255)
                                     aTCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATcca                                                                                                                 >  2:4104888/1‑149 (MQ=255)
                                     aTCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATcca                                                                                                                 <  2:897469/149‑1 (MQ=255)
                                      tCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATccac                                                                                                                >  1:1969712/1‑149 (MQ=255)
                                      tCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCAGCCTCGACGCTGACATCAAAGCGATGGCAAAAATccac                                                                                                                >  1:4055464/1‑149 (MQ=255)
                                       cAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATccacc                                                                                                               <  2:2144810/149‑1 (MQ=255)
                                        aaTTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATccacca                                                                                                              >  2:233879/1‑149 (MQ=255)
                                            ccGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACctc                                                                                                          <  1:1673748/149‑1 (MQ=255)
                                             aggTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACctct                                                                                                         >  1:2248881/2‑149 (MQ=255)
                                             aggTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCACAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACctct                                                                                                         >  1:4589065/2‑149 (MQ=255)
                                               gTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCttt                                                                                                       >  2:2973122/1‑149 (MQ=255)
                                                 tGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGa                                                                                                     <  1:1592062/149‑1 (MQ=255)
                                                  gCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGAc                                                                                                    >  1:4140034/1‑149 (MQ=255)
                                                   ccTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGAcc                                                                                                   <  2:4453229/149‑1 (MQ=255)
                                                     tCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCg                                                                                                 >  1:3951833/1‑149 (MQ=255)
                                                      cGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGt                                                                                                <  2:3863149/149‑1 (MQ=255)
                                                       ggCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTg                                                                                               <  1:3416971/149‑1 (MQ=255)
                                                       ggCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTg                                                                                               <  2:2281580/149‑1 (MQ=255)
                                                        gCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGa                                                                                              >  1:4256319/1‑149 (MQ=255)
                                                          aaaTTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGaaa                                                                                            <  2:2212927/149‑1 (MQ=255)
                                                            aTTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAata                                                                                          <  1:3489329/149‑1 (MQ=255)
                                                             ttCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAatat                                                                                         <  2:4503200/149‑1 (MQ=255)
                                                             ttCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAatat                                                                                         <  2:2054006/149‑1 (MQ=255)
                                                                  gcgcCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTgc                                                                                    >  1:1138245/1‑149 (MQ=255)
                                                                  gcgcCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTgc                                                                                    >  1:3473696/1‑149 (MQ=255)
                                                                      cTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCaa                                                                                <  1:4369085/149‑1 (MQ=255)
                                                                        gAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCaaaa                                                                              >  2:4376503/1‑149 (MQ=255)
                                                                        gAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCaaaa                                                                              <  1:2287095/149‑1 (MQ=255)
                                                                         aGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTGTGCCTCGACGCTGACATCAAAGCGTTGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAg                                                                             <  1:4611023/149‑1 (MQ=255)
                                                                         aGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAg                                                                             <  1:868327/149‑1 (MQ=255)
                                                                           cGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGc                                                                           <  1:307867/149‑1 (MQ=255)
                                                                            ggggTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCa                                                                          <  1:1255786/149‑1 (MQ=255)
                                                                               gTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAAt                                                                       <  1:569619/149‑1 (MQ=255)
                                                                               gTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAAt                                                                       >  1:323566/1‑149 (MQ=255)
                                                                               gTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAAt                                                                       >  2:2106179/1‑149 (MQ=255)
                                                                                 tGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGc                                                                     >  1:1882395/1‑149 (MQ=255)
                                                                                  ggTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCg                                                                    <  2:3939994/149‑1 (MQ=255)
                                                                                   gtgtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGa                                                                   <  1:3724951/149‑1 (MQ=255)
                                                                                   gtgtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGa                                                                   <  2:1522869/149‑1 (MQ=255)
                                                                                    tgtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAg                                                                  <  2:4614040/149‑1 (MQ=255)
                                                                                    tgtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAg                                                                  >  1:3896820/1‑149 (MQ=255)
                                                                                    tgtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAg                                                                  >  2:1936189/1‑149 (MQ=255)
                                                                                     gtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAgg                                                                 <  1:4635837/149‑1 (MQ=255)
                                                                                     gtACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAgg                                                                 >  1:1022960/1‑149 (MQ=255)
                                                                                        cAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGt                                                              <  1:107371/149‑1 (MQ=255)
                                                                                              aTCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCa                                                        <  2:1744239/149‑1 (MQ=255)
                                                                                               tCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCAt                                                       >  2:4240900/1‑149 (MQ=255)
                                                                                                  gCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGAc                                                    >  2:3164853/1‑149 (MQ=255)
                                                                                                  gCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGAc                                                    <  1:4104888/149‑1 (MQ=255)
                                                                                                      cATTGACTTGCACATCATCTTCCTGGACGCCAAGTTGCCCCAGCTCCGCCTAGACGCTGACATCAAAGCTATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGat                                                <  2:573906/149‑2 (MQ=255)
                                                                                                       aTTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAgc                                               <  2:363889/149‑1 (MQ=255)
                                                                                                       aTTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAgc                                               >  2:674203/1‑149 (MQ=255)
                                                                                                          gACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCa                                            <  1:3317073/149‑1 (MQ=255)
                                                                                                               gCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCa                                       >  1:1214006/1‑149 (MQ=255)
                                                                                                                cacaTCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAg                                      >  1:757888/1‑149 (MQ=255)
                                                                                                                  catcatCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaa                                    <  1:721821/149‑1 (MQ=255)
                                                                                                                  catcatCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaa                                    <  2:3709925/149‑1 (MQ=255)
                                                                                                                   atcatcTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaac                                   <  1:4516100/149‑1 (MQ=255)
                                                                                                                   atcatcTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaac                                   <  1:2181394/149‑1 (MQ=255)
                                                                                                                    tcatcTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaacg                                  >  1:941496/1‑149 (MQ=255)
                                                                                                                      atcTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAgaacgaa                                <  1:4449021/149‑1 (MQ=255)
                                                                                                                         ttCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCt                             <  1:1155581/149‑1 (MQ=255)
                                                                                                                         gtCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCt                             <  1:4039477/148‑1 (MQ=255)
                                                                                                                          tCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCtt                            <  2:812881/149‑1 (MQ=255)
                                                                                                                           ccTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTc                           <  2:4676193/149‑1 (MQ=255)
                                                                                                                           ccTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTc                           <  2:311580/149‑1 (MQ=255)
                                                                                                                            cTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCt                          <  1:1919416/149‑1 (MQ=255)
                                                                                                                             tGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTg                         >  1:1708599/1‑149 (MQ=255)
                                                                                                                             tGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTg                         >  1:2811277/1‑149 (MQ=255)
                                                                                                                               gccgccAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAg                       <  2:3474764/149‑1 (MQ=255)
                                                                                                                                    cAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGata                  <  2:1764991/149‑1 (MQ=255)
                                                                                                                                     aaGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGatat                 >  1:631207/1‑149 (MQ=255)
                                                                                                                                     aaGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGatat                 <  1:3400465/149‑1 (MQ=255)
                                                                                                                                       gTTGCCCCAGCTCTGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCg               >  2:3797824/1‑149 (MQ=255)
                                                                                                                                        ttGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGt              <  1:3587346/149‑1 (MQ=255)
                                                                                                                                          gCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCa            >  1:249777/1‑149 (MQ=255)
                                                                                                                                              cAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTGTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAAtt        <  1:2951253/149‑1 (MQ=255)
                                                                                                                                                gCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTc      >  2:502387/1‑149 (MQ=255)
                                                                                                                                                gCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTc      >  2:1048912/1‑149 (MQ=255)
                                                                                                                                                gCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTc      >  1:3763181/1‑149 (MQ=255)
                                                                                                                                                 cTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCg     <  2:1672164/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACATCCCATGAGCCGTGCCGATGGCGACCGCCAGGGAATCAATTCCGGTTGCCTCGGCAAATTCACGCGCCTGAGCGGGGTTGGTGTACAACGCATCGGCTTCATTGACTTGCACATCATCTTCCTGGCCGCCAAGTTGCCCCAGCTCCGCCTCGACGCTGACATCAAAGCGATGGCAAAAATCCACCACCTCTTTGACCCGTGAAATATTTTGCGCAAAAGGCAAATGCGAGGCGTCAATCATGACTGAGCGCACGCCAGAACGAACCTTCTGAGCGATATCGTCAAATTTCGTGT  >  CP009273/2170137‑2170433

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: ‑