Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,022,182 G→T 7.4% P97Q (CCG→CAG)  mgsA ← methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,022,1820GT 7.4% 1226.0 / 43.4 320P97Q (CCG→CAG) mgsAmethylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base G (125/170);  new base T (12/13);  total (137/183)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

GATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTA  >  CP009273/1022033‑1022331
                                                                                                                                                     |                                                                                                                                                     
 aTAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGcgg                                                                                                                                                      <  2:1875532/150‑1 (MQ=255)
    aCGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGcggcac                                                                                                                                                   >  2:3003665/1‑150 (MQ=255)
    aCGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGcggcac                                                                                                                                                   >  2:2661999/1‑150 (MQ=255)
      gCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGcggcacgg                                                                                                                                                 >  2:3486139/1‑150 (MQ=255)
       cTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGcggcacggc                                                                                                                                                >  1:1831986/1‑150 (MQ=255)
       cTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGcggcacggc                                                                                                                                                >  1:4405840/1‑150 (MQ=255)
           tAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCAttt                                                                                                                                            <  2:3108036/150‑1 (MQ=255)
           tAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCAttt                                                                                                                                            >  1:4229103/1‑150 (MQ=255)
           tAATCGGGGATCAGAATATCGACCGCGTCGATGAAAAGCGGCGACTGGATTATGAAGACTGCCGTTGCCACCTTGGTGGCGACCGGAATGTTACATACCGTCGCCATACGCAGCAAGGCTTTCACGGCAGGATCGTGCGGCACGGCAttt                                                                                                                                            >  1:3599889/1‑150 (MQ=255)
            aaTCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCTGCACGGCATTTa                                                                                                                                           >  1:836231/1‑150 (MQ=255)
            aaTCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTa                                                                                                                                           <  2:3202113/150‑1 (MQ=255)
              tCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGt                                                                                                                                         <  2:2476406/150‑1 (MQ=255)
              tCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGt                                                                                                                                         >  1:2832033/1‑150 (MQ=255)
                ggggATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCTGCACGGCAGTGAGTgg                                                                                                                                       >  2:5111369/1‑150 (MQ=255)
                ggggATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTgg                                                                                                                                       >  1:4131509/1‑150 (MQ=255)
                 gggATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGa                                                                                                                                      <  1:814869/150‑1 (MQ=255)
                   gATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATc                                                                                                                                    <  1:4449002/150‑1 (MQ=255)
                    aTCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCTGCACGGCATTTAGTGGATcc                                                                                                                                   >  1:1633922/1‑150 (MQ=255)
                    aTCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATcc                                                                                                                                   <  2:1297875/150‑1 (MQ=255)
                     tCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATccc                                                                                                                                  <  2:1815668/150‑1 (MQ=255)
                      cAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCa                                                                                                                                 >  2:2513434/1‑150 (MQ=255)
                      cAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCa                                                                                                                                 <  2:1329603/150‑1 (MQ=255)
                        gAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCaga                                                                                                                               <  1:503281/150‑1 (MQ=255)
                         aaTATCGACCGCGTCGTTGAAAGGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCAGACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCagaa                                                                                                                              <  2:3601655/150‑1 (MQ=255)
                           tatCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCagaaga                                                                                                                            >  2:4343055/1‑150 (MQ=255)
                              cGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGaaaa                                                                                                                         <  2:3164701/150‑1 (MQ=255)
                               gACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAAt                                                                                                                        >  2:41471/1‑150 (MQ=255)
                               gACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAAt                                                                                                                        <  2:1557168/150‑1 (MQ=255)
                               gACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAAt                                                                                                                        >  1:2880088/1‑150 (MQ=255)
                                aCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCTGCACGGCATTTAGTGGATCCCAGAAGAAAATc                                                                                                                       <  2:2163994/150‑1 (MQ=255)
                                aCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCTGCACGGCATTTAGTGGATCCCAGAAGAAAATc                                                                                                                       >  1:4521837/1‑150 (MQ=255)
                                aCCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATc                                                                                                                       >  1:2744598/1‑150 (MQ=255)
                                  cgcgTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCTGCACGGCATTTAGTGGATCCCAGAAGAAAATCaa                                                                                                                     >  1:5076487/1‑150 (MQ=255)
                                  cgcgTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCaa                                                                                                                     >  2:461737/1‑150 (MQ=255)
                                      tcgtTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACa                                                                                                                 <  1:2039675/150‑1 (MQ=255)
                                        gtTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATc                                                                                                               >  1:3221237/1‑150 (MQ=255)
                                        gtTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATc                                                                                                               >  1:15204/1‑150 (MQ=255)
                                        gtTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATc                                                                                                               >  1:1740553/1‑150 (MQ=255)
                                        gtTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATc                                                                                                               <  1:1925672/150‑1 (MQ=255)
                                        gtTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATc                                                                                                               <  1:4106403/150‑1 (MQ=255)
                                          tGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCaa                                                                                                             >  2:3516060/1‑150 (MQ=255)
                                          tGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCaa                                                                                                             >  2:4064088/1‑150 (MQ=255)
                                            aaaTGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAAtt                                                                                                           <  2:4504666/150‑1 (MQ=255)
                                                   gcgACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCtt                                                                                                    >  2:917871/1‑150 (MQ=255)
                                                    cgACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTc                                                                                                   <  1:4452470/150‑1 (MQ=255)
                                                     gACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGTCGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCt                                                                                                  <  2:4699524/150‑1 (MQ=255)
                                                     gACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCt                                                                                                  <  1:3356957/150‑1 (MQ=255)
                                                       cTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTga                                                                                                <  2:2466544/150‑1 (MQ=255)
                                                        tGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCTGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTgag                                                                                               >  2:1124785/1‑150 (MQ=255)
                                                        tGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTgag                                                                                               >  1:2818785/1‑150 (MQ=255)
                                                          gATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGAt                                                                                             <  2:820326/150‑1 (MQ=255)
                                                           aTTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGACTGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATc                                                                                            <  2:2634636/150‑1 (MQ=255)
                                                               tGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATg                                                                                        >  2:3723138/1‑150 (MQ=255)
                                                               tGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATg                                                                                        >  1:461926/1‑150 (MQ=255)
                                                                 aaGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATgcg                                                                                      <  1:706029/150‑1 (MQ=255)
                                                                    tCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCa                                                                                   >  2:1696436/1‑150 (MQ=255)
                                                                      tGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAAc                                                                                 >  1:1278722/1‑150 (MQ=255)
                                                                        ccGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACct                                                                               <  2:4236683/150‑1 (MQ=255)
                                                                               cACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTc                                                                        >  1:3941593/1‑150 (MQ=255)
                                                                                   ttGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAccc                                                                    <  2:1441204/150‑1 (MQ=255)
                                                                                   ttGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAccc                                                                    <  1:1056560/150‑1 (MQ=255)
                                                                                    tggtggCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCAcccc                                                                   >  1:5127226/1‑150 (MQ=255)
                                                                                         gCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCAtt                                                              >  1:294367/1‑150 (MQ=255)
                                                                                           gACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTgg                                                            <  1:2375611/150‑1 (MQ=255)
                                                                                            aCCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTggg                                                           >  1:4694053/1‑150 (MQ=255)
                                                                                                 aaTGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACt                                                      >  1:4738274/1‑150 (MQ=255)
                                                                                                   tGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCa                                                    <  2:679621/150‑1 (MQ=255)
                                                                                                    gTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCaa                                                   >  2:1434648/1‑150 (MQ=255)
                                                                                                      tCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACa                                                 >  2:5137618/1‑150 (MQ=255)
                                                                                                      tCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACa                                                 <  1:2441264/150‑1 (MQ=255)
                                                                                                        cATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATc                                               <  1:256562/150‑1 (MQ=255)
                                                                                                        cATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATc                                               >  2:1835949/1‑150 (MQ=255)
                                                                                                          tACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATcgc                                             <  2:2753130/150‑1 (MQ=255)
                                                                                                                cGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGAc                                       <  1:3027909/150‑1 (MQ=255)
                                                                                                                 gCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACg                                      >  1:4862630/1‑150 (MQ=255)
                                                                                                                   cAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGtt                                    <  1:1134288/150‑1 (MQ=255)
                                                                                                                    aGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTc                                   >  1:2068019/1‑150 (MQ=255)
                                                                                                                     gACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCa                                  >  1:7913/1‑150 (MQ=255)
                                                                                                                       cGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATg                                >  2:4947221/1‑150 (MQ=255)
                                                                                                                         cagcaAGTCTTTCACGTCAGGATCGTGCTGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGcc                              <  1:4933815/150‑1 (MQ=255)
                                                                                                                          agcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCg                             >  2:2826722/1‑150 (MQ=255)
                                                                                                                           gcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCgg                            <  1:3812440/150‑1 (MQ=255)
                                                                                                                           gcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCgg                            >  1:317129/1‑150 (MQ=255)
                                                                                                                             aaGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTc                          >  1:235199/1‑150 (MQ=255)
                                                                                                                              aGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcg                         <  1:3003665/150‑1 (MQ=255)
                                                                                                                                gCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgcg                       <  2:2738671/150‑1 (MQ=255)
                                                                                                                                 cTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgcgc                      >  1:266220/1‑150 (MQ=255)
                                                                                                                                    tCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCggg                   <  1:3601725/150‑1 (MQ=255)
                                                                                                                                      aCGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaa                 >  2:2253418/1‑150 (MQ=255)
                                                                                                                                      aCGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCGGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaa                 <  1:4166708/150‑1 (MQ=255)
                                                                                                                                       cGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCATAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                <  2:1074438/150‑1 (MQ=255)
                                                                                                                                        gTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAAt               >  2:1203414/1‑150 (MQ=255)
                                                                                                                                         tCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAAtt              >  2:5017820/1‑150 (MQ=255)
                                                                                                                                             gATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAACTTAAGt          >  2:2933139/1‑150 (MQ=255)
                                                                                                                                             gATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGt          <  1:4335291/150‑1 (MQ=255)
                                                                                                                                             gATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGt          >  1:756060/1‑150 (MQ=255)
                                                                                                                                              aTCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGtt         >  2:615948/1‑150 (MQ=255)
                                                                                                                                              aTCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGtt         <  1:2150955/150‑1 (MQ=255)
                                                                                                                                               tCGTGCTGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTa        <  2:3226213/150‑1 (MQ=255)
                                                                                                                                                 gTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAcc      <  1:311305/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAACGCTGATAATCGGGGATCAGAATATCGACCGCGTCGTTGAAATGCGGCGACTGGATTATGAAGTCTGCCGTTGCCACGTTGGTGGCGACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTCACCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTA  >  CP009273/1022033‑1022331

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