Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,340,901 A→T 38.1% S33S (TCT→TCA gmr ← cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,340,9010AT38.1% 214.5 / 280.6 223S33S (TCT→TCAgmrcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of Rnase II stability
Reads supporting (aligned to +/- strand):  ref base A (67/71);  new base T (34/51);  total (101/122)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

TTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACA  >  CP009273/1340752‑1341050
                                                                                                                                                     |                                                                                                                                                     
  tCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGa                                                                                                                                                     <  2:3193361/150‑1 (MQ=255)
  tCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGa                                                                                                                                                     >  2:1224180/1‑150 (MQ=255)
   cGTCCTACCAGATGCACAGACAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGAGGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCACTGTACGATCAGGTTCTGTGGTTTCGGTGGTTGaa                                                                                                                                                    >  2:2883010/1‑150 (MQ=255)
      ccTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAg                                                                                                                                                 >  1:618794/1‑150 (MQ=255)
       cTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAgc                                                                                                                                                <  2:1114683/150‑1 (MQ=255)
        tACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAgcg                                                                                                                                               <  1:40037/150‑1 (MQ=255)
             gATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCaaaa                                                                                                                                          >  2:20159/1‑150 (MQ=255)
              aTGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAAc                                                                                                                                         >  1:479832/1‑150 (MQ=255)
               tGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACa                                                                                                                                        >  2:860433/1‑150 (MQ=255)
                gaacaGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGATGAAAAGCGCAAAACAt                                                                                                                                       >  2:1060049/3‑150 (MQ=255)
                gCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACAt                                                                                                                                       >  1:690803/1‑150 (MQ=255)
                gCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACAt                                                                                                                                       <  1:1153430/150‑1 (MQ=255)
                gCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACAt                                                                                                                                       <  1:1775486/150‑1 (MQ=255)
                 cacaGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATc                                                                                                                                      <  2:1656467/150‑1 (MQ=255)
                  acaGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCg                                                                                                                                     <  1:1246631/150‑1 (MQ=255)
                  acaGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCg                                                                                                                                     >  1:3507549/1‑150 (MQ=255)
                    aGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgct                                                                                                                                   >  2:4132406/1‑150 (MQ=255)
                     gAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCgctg                                                                                                                                  <  2:2286146/150‑1 (MQ=255)
                     gAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCgctg                                                                                                                                  <  2:3264040/150‑1 (MQ=255)
                       aaGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCgctgct                                                                                                                                >  2:3875506/1‑150 (MQ=255)
                            cGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCg                                                                                                                           >  2:3589526/1‑150 (MQ=255)
                            cGCTGTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACAGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCTGTGGTAGAAAAGCGCAAAACATCGCTGCTTTcct                                                                                                                           <  2:2686420/150‑2 (MQ=255)
                             gCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGt                                                                                                                          >  2:2444148/1‑150 (MQ=255)
                              cTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTc                                                                                                                         <  1:1754564/150‑1 (MQ=255)
                               tttCATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCa                                                                                                                        <  1:4135211/150‑1 (MQ=255)
                                ttCATCGACGTTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTTATGCGAGCAGCCAGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCaa                                                                                                                       >  2:341261/1‑150 (MQ=255)
                                 tCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAAc                                                                                                                      >  2:1903803/1‑150 (MQ=255)
                                  cATCGACGGTCAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAAcc                                                                                                                     >  2:1092664/1‑150 (MQ=255)
                                   aTCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCg                                                                                                                    <  2:3010552/150‑1 (MQ=255)
                                         ggTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTa                                                                                                              >  1:881822/1‑150 (MQ=255)
                                            cAGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATgg                                                                                                           <  1:2631047/150‑1 (MQ=255)
                                             aGACTCATCATCAGGCTGGAGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGa                                                                                                          <  2:3872594/150‑1 (MQ=255)
                                             aGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGa                                                                                                          >  2:7603/1‑150 (MQ=255)
                                             aGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGa                                                                                                          <  2:3231527/150‑1 (MQ=255)
                                              gACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAt                                                                                                         <  1:4133661/150‑1 (MQ=255)
                                               aCTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAtt                                                                                                        <  1:2872151/150‑1 (MQ=255)
                                               aCTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGAtt                                                                                                        <  1:1557544/150‑1 (MQ=255)
                                                cTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATtg                                                                                                       <  2:2780283/150‑1 (MQ=255)
                                                 tcatcatcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATtgt                                                                                                      >  1:2667782/1‑150 (MQ=255)
                                                     catcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGa                                                                                                  <  1:2029804/150‑1 (MQ=255)
                                                      atcaGGCTGGAGGTGATGACCGTCATTTCCNTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGAt                                                                                                 >  1:3854352/1‑150 (MQ=255)
                                                      atcaGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGAt                                                                                                 <  1:2936901/150‑1 (MQ=255)
                                                        caGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATcc                                                                                               >  2:1982492/1‑150 (MQ=255)
                                                         aGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGACAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCg                                                                                              >  2:614941/1‑150 (MQ=255)
                                                              ggAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGaaa                                                                                         <  1:3597072/150‑1 (MQ=255)
                                                                aGGTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGt                                                                                       <  1:1015968/150‑1 (MQ=255)
                                                                 ggTTATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGtt                                                                                      >  2:3389566/1‑150 (MQ=255)
                                                                        aCCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAAc                                                                               <  2:3562067/150‑1 (MQ=255)
                                                                        aCCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAAc                                                                               <  1:920710/150‑1 (MQ=255)
                                                                           gTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGtt                                                                            >  2:1670472/1‑150 (MQ=255)
                                                                            tCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtg                                                                           <  1:139429/150‑1 (MQ=255)
                                                                            tCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGtgg                                                                           >  1:4084731/1‑148 (MQ=255)
                                                                            tCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtg                                                                           >  1:1984800/1‑150 (MQ=255)
                                                                            tCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtg                                                                           <  2:2131983/150‑1 (MQ=255)
                                                                            tCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtg                                                                           >  1:1342117/1‑150 (MQ=255)
                                                                             cATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtgt                                                                          >  2:3119420/1‑150 (MQ=255)
                                                                             cATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTtgt                                                                          <  2:2083699/150‑1 (MQ=255)
                                                                               tttCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTgg                                                                        <  2:1280892/150‑1 (MQ=255)
                                                                                    cTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTcc                                                                   >  2:3074691/1‑150 (MQ=255)
                                                                                    cTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTcc                                                                   <  1:2572306/150‑1 (MQ=255)
                                                                                     tGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTccc                                                                  <  1:2613526/150‑1 (MQ=255)
                                                                                      gATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCt                                                                 <  2:254799/150‑1 (MQ=255)
                                                                                      gATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCt                                                                 >  2:1608400/1‑150 (MQ=255)
                                                                                         gCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAAc                                                              <  2:1570672/150‑1 (MQ=255)
                                                                                         gCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAAc                                                              <  2:3667588/150‑1 (MQ=255)
                                                                                            agcagcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCTCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGt                                                           >  2:3317816/1‑150 (MQ=255)
                                                                                            agcagcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGt                                                           >  2:738768/1‑150 (MQ=255)
                                                                                            agcagcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGt                                                           <  2:3519778/150‑1 (MQ=255)
                                                                                              cagtCTGTTCGGCAGATAACTGAAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGAGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGttt                                                         >  2:2734004/1‑150 (MQ=255)
                                                                                              cagcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGttt                                                         <  1:370109/150‑1 (MQ=255)
                                                                                               agcCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGtttt                                                        <  2:3366609/150‑1 (MQ=255)
                                                                                                gcCTGTTCGGCAGATAACTGCAATGTACGGTCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTc                                                       >  2:555283/1‑150 (MQ=255)
                                                                                                  cTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCat                                                     >  2:559834/1‑150 (MQ=255)
                                                                                                  cTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCat                                                     <  2:2857327/150‑1 (MQ=255)
                                                                                                    gTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatga                                                   <  1:2290136/150‑1 (MQ=255)
                                                                                                    gTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCatga                                                   >  2:404599/1‑150 (MQ=255)
                                                                                                               aaCTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcctc                                        <  1:3248012/150‑1 (MQ=255)
                                                                                                                aCTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcctcc                                       <  1:2129095/150‑1 (MQ=255)
                                                                                                                 cTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCg                                      <  1:4162412/150‑1 (MQ=255)
                                                                                                                    cAATGTACGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAAACGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGAGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGtt                                   >  2:2693739/1‑150 (MQ=255)
                                                                                                                           cGATCAGGTTCTGTGGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                            >  1:1708505/1‑150 (MQ=255)
                                                                                                                             aTCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                          >  1:1100524/1‑150 (MQ=255)
                                                                                                                                 ggTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTcc                      <  2:4165924/150‑1 (MQ=255)
                                                                                                                                      tgtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGt                 >  1:3295084/1‑150 (MQ=255)
                                                                                                                                       gtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGtt                <  2:2224062/150‑1 (MQ=255)
                                                                                                                                        tgGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTta               <  1:1456824/150‑1 (MQ=255)
                                                                                                                                        tgGTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTta               <  1:1907780/150‑1 (MQ=255)
                                                                                                                                         ggTTTCGGTGGTTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaa              <  1:3779823/150‑1 (MQ=255)
                                                                                                                                         ggTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaa              >  2:2155943/1‑150 (MQ=255)
                                                                                                                                         ggTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaa              <  2:1437787/150‑1 (MQ=255)
                                                                                                                                         ggTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaa              <  2:1411259/150‑1 (MQ=255)
                                                                                                                                          gTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa             >  2:2577754/1‑150 (MQ=255)
                                                                                                                                              cGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa         >  2:587201/1‑150 (MQ=255)
                                                                                                                                                 tggtTGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAgg      <  2:1824557/150‑1 (MQ=255)
                                                                                                                                                    tAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaac   >  1:3411915/1‑150 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAATTAAAAGGaaca  >  1:1899974/1‑150 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAATGGaaca  >  1:1108216/1‑150 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca  >  1:2445748/1‑150 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca  >  1:745962/1‑150 (MQ=255)
                                                                                                                                                     aGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca  >  2:1392150/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACA  >  CP009273/1340752‑1341050

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: