Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,127,560 T→G 5.9% Q525P (CAG→CCG)  wzc ← colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,127,5600TG 5.9% 491.6 / 21.6 169Q525P (CAG→CCG) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (77/82);  new base G (3/7);  total (80/89)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 171 total aligned reads displayed.

AGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCTTTC  >  CP009273/2127416‑2127706
                                                                                                                                                |                                                                                                                                                  
  cTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGcc                                                                                                                                               >  1:460797/1‑148 (MQ=255)
     gTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGc                                                                                                                                            <  2:1496354/148‑1 (MQ=255)
      tGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGca                                                                                                                                           >  2:2850896/1‑148 (MQ=255)
        ggTGTAGCCTTTGCGCATATCGCAGTCGTTCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAACTGTTTTACCACTTTACGGGCTAAACCCGGTCATCATCAACCCATTGTTCTGTGCCTGcatc                                                                                                                                         <  1:2762019/148‑1 (MQ=255)
        ggTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGcatc                                                                                                                                         <  2:1363309/148‑1 (MQ=255)
               ccTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGa                                                                                                                                  >  1:780391/1‑148 (MQ=255)
                 tttGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAg                                                                                                                                >  2:1793811/1‑148 (MQ=255)
                 tttGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAg                                                                                                                                <  1:1279799/148‑1 (MQ=255)
                 tttGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAg                                                                                                                                >  2:1609929/1‑148 (MQ=255)
                        atatCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAAc                                                                                                                         <  2:1959558/148‑1 (MQ=255)
                            cGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGt                                                                                                                     >  1:1230217/1‑148 (MQ=255)
                             gCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTa                                                                                                                    <  2:67239/148‑1 (MQ=255)
                                gTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGc                                                                                                                 >  2:1453534/1‑148 (MQ=255)
                                 tCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCa                                                                                                                <  1:2193227/148‑1 (MQ=255)
                                   gATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGa                                                                                                              <  1:3008041/148‑1 (MQ=255)
                                    aTCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGAc                                                                                                             >  1:2408040/1‑148 (MQ=255)
                                     tCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACt                                                                                                            >  1:474868/1‑148 (MQ=255)
                                       aacaacACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTAc                                                                                                          <  2:1533808/148‑1 (MQ=255)
                                        acaacaCGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACg                                                                                                         >  2:192769/1‑148 (MQ=255)
                                         caacaCGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACgg                                                                                                        <  1:2867355/148‑1 (MQ=255)
                                         caacaCGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACgg                                                                                                        >  2:2607991/1‑148 (MQ=255)
                                           acacGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCGGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGAt                                                                                                      >  2:2378210/1‑148 (MQ=255)
                                            cacGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATg                                                                                                     >  1:1984347/1‑148 (MQ=255)
                                            cacGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATg                                                                                                     >  1:420407/1‑148 (MQ=255)
                                            cacGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATg                                                                                                     >  2:1649721/1‑148 (MQ=255)
                                             acGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATgg                                                                                                    <  2:2746872/148‑1 (MQ=255)
                                             acGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATgg                                                                                                    <  1:1878488/148‑1 (MQ=255)
                                               gcgTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCt                                                                                                  <  2:2518419/148‑1 (MQ=255)
                                               gcgTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCt                                                                                                  <  2:3483/148‑1 (MQ=255)
                                               gcgTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCt                                                                                                  <  2:3492/148‑1 (MQ=255)
                                                cgTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCtt                                                                                                 <  2:1652324/148‑1 (MQ=255)
                                                  tttATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCa                                                                                               <  2:2063441/148‑1 (MQ=255)
                                                   ttATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCaa                                                                                              >  2:2912695/1‑148 (MQ=255)
                                                   ttATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCaa                                                                                              <  1:2380904/148‑1 (MQ=255)
                                                   ttATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCaa                                                                                              >  2:1263349/1‑148 (MQ=255)
                                                   ttATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCGGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCaa                                                                                              <  2:1013247/148‑1 (MQ=255)
                                                       tGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGc                                                                                          >  1:2325498/1‑148 (MQ=255)
                                                       tGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGc                                                                                          >  1:1343417/1‑148 (MQ=255)
                                                               cTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCgg                                                                                  <  2:1792046/148‑1 (MQ=255)
                                                               cTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCgg                                                                                  >  2:1334627/1‑148 (MQ=255)
                                                                tGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAAAACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACTCAGACTACGGATGGCTTCAATCGCCAGATCGGt                                                                                 >  1:2302558/1‑148 (MQ=255)
                                                                 gATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGtt                                                                                <  2:2366891/148‑1 (MQ=255)
                                                                   tCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTgg                                                                              >  2:1954284/1‑148 (MQ=255)
                                                                      ccgccgCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGACGGCTTCAATCGCCAGATCGGTTGGAtt                                                                           <  1:3299863/148‑1 (MQ=255)
                                                                        gccgccAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTcc                                                                         <  1:172880/148‑1 (MQ=255)
                                                                              aGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCg                                                                   >  2:2687063/1‑148 (MQ=255)
                                                                              aGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCg                                                                   >  2:2849104/1‑148 (MQ=255)
                                                                                gTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGcc                                                                 >  2:793592/1‑148 (MQ=255)
                                                                                gTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGcc                                                                 <  1:540908/148‑1 (MQ=255)
                                                                                gTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGcc                                                                 <  1:540851/148‑1 (MQ=255)
                                                                                gTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGcc                                                                 >  2:2452551/1‑148 (MQ=255)
                                                                                 ttGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCa                                                                <  2:3085697/148‑1 (MQ=255)
                                                                                 ttGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCa                                                                >  2:774372/1‑148 (MQ=255)
                                                                                     cgcAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAg                                                            <  2:590988/148‑1 (MQ=255)
                                                                                      gcAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGc                                                           <  1:2092812/148‑1 (MQ=255)
                                                                                          aCAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTCGCTGGc                                                       >  2:1774064/1‑148 (MQ=255)
                                                                                          aCAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGc                                                       >  1:693601/1‑148 (MQ=255)
                                                                                          aCAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGc                                                       >  2:2748088/1‑148 (MQ=255)
                                                                                          aCAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGc                                                       >  1:1158576/1‑148 (MQ=255)
                                                                                           cAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGct                                                      >  1:1119260/1‑148 (MQ=255)
                                                                                           cAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGct                                                      <  2:1728040/148‑1 (MQ=255)
                                                                                            aaaGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGctc                                                     <  2:2018031/148‑1 (MQ=255)
                                                                                               ggTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGt                                                  >  2:2335272/1‑148 (MQ=255)
                                                                                                gTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGtt                                                 >  2:556101/1‑148 (MQ=255)
                                                                                                gTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGtt                                                 >  1:1639384/1‑148 (MQ=255)
                                                                                                 ttttACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGttt                                                <  2:2117044/148‑1 (MQ=255)
                                                                                                 ttttACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGttt                                                >  1:2686385/1‑148 (MQ=255)
                                                                                                     aCCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAg                                            <  2:1230217/148‑1 (MQ=255)
                                                                                                     aCCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAg                                            <  2:2094650/148‑1 (MQ=255)
                                                                                                     aCCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAg                                            >  1:2643826/1‑148 (MQ=255)
                                                                                                        aaTTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGt                                         >  1:2119879/1‑148 (MQ=255)
                                                                                                        aaTTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGt                                         >  2:695522/1‑148 (MQ=255)
                                                                                                        aaTTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGt                                         >  2:1024030/1‑148 (MQ=255)
                                                                                                           tGACGGGCTAACCCCGGTCATCATCAACACATTGTTCGGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTa                                      <  2:739201/148‑1 (MQ=255)
                                                                                                            gACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTaa                                     <  1:3246909/148‑1 (MQ=255)
                                                                                                             aCGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAAt                                    <  1:2372605/148‑1 (MQ=255)
                                                                                                              cGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATc                                   >  2:2958984/1‑148 (MQ=255)
                                                                                                               gggCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATcc                                  <  2:976119/148‑1 (MQ=255)
                                                                                                                ggCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATccc                                 <  2:2147105/148‑1 (MQ=255)
                                                                                                                ggCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATccc                                 <  2:2146149/148‑1 (MQ=255)
                                                                                                                  cTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCtt                               <  1:1610467/148‑1 (MQ=255)
                                                                                                                     aCCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGa                            <  1:623169/148‑1 (MQ=255)
                                                                                                                     aCCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGa                            <  1:2192302/148‑1 (MQ=255)
                                                                                                                      ccccGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGAt                           >  1:1899010/1‑148 (MQ=255)
                                                                                                                       cccGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATg                          <  1:1387697/148‑1 (MQ=255)
                                                                                                                             catcatcaACACATTGTTCGGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTg                    <  2:39599/148‑1 (MQ=255)
                                                                                                                              atcatcaACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGa                   <  1:899045/148‑1 (MQ=255)
                                                                                                                                catcaACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACg                 >  2:1728680/1‑148 (MQ=255)
                                                                                                                                 atcaACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGc                >  1:1091489/1‑148 (MQ=255)
                                                                                                                                 atcaACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGc                <  1:1453534/148‑1 (MQ=255)
                                                                                                                                     acacATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATc            <  2:1016949/148‑1 (MQ=255)
                                                                                                                                     acacATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATc            <  1:1609929/148‑1 (MQ=255)
                                                                                                                                     acacATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATc            <  2:1016944/148‑1 (MQ=255)
                                                                                                                                      cacaTTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCa           >  2:925248/1‑148 (MQ=255)
                                                                                                                                          ttGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCAcgcg       >  2:408912/1‑148 (MQ=255)
                                                                                                                                          ttGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCAcgcg       >  2:408267/1‑148 (MQ=255)
                                                                                                                                           tGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACgcgc      >  1:1208206/1‑148 (MQ=255)
                                                                                                                                             ttCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCtt    <  2:2548874/148‑1 (MQ=255)
                                                                                                                                               cTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCTTTc  <  1:1649721/148‑1 (MQ=255)
                                                                                                                                               cTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCTTTc  <  1:1607902/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGGTTGGCGCAGACAAAGGTTTTACCAATTGACGGGCTAACCCCGGTCATCATCAACACATTGTTCTGTGCCTGCATCATCGCGAAGTGCAAACTGGTACGCAGACTACGGATGGCTTCAATCGCCAGATCGGTTGGATTCCCCACCGCCAGTAGCTGGCTCTGTTTATAGCGTTTAATCCCTTTGATGGTTTTGACGCTATCACGCGCTTTC  >  CP009273/2127416‑2127706

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: