Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,277,142 C→A 47.0% P214Q (CCG→CAG)  kbaY → tagatose 6‑phosphate aldolase 1, kbaY subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,277,1420CA47.0% 151.3 / 314.4 247P214Q (CCG→CAG) kbaYtagatose 6‑phosphate aldolase 1, kbaY subunit
Reads supporting (aligned to +/- strand):  ref base C (67/64);  new base A (49/67);  total (116/131)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 247 total aligned reads displayed.

AACTGACTGGCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATTATATGCGCGTC  >  CP009273/3276995‑3277287
                                                                                                                                                   |                                                                                                                                                 
aaCTGACTGGCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCg                                                                                                                                                  >  1:1998728/1‑149 (MQ=255)
  ctgactgGCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGa                                                                                                                                                <  1:3871622/149‑1 (MQ=255)
    gactgGCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATg                                                                                                                                              >  1:2688526/1‑149 (MQ=255)
       tgGCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAAt                                                                                                                                           <  2:1571320/149‑1 (MQ=255)
        ggCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAAtt                                                                                                                                          >  1:494450/1‑149 (MQ=255)
         gCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAAttt                                                                                                                                         <  2:1816122/149‑1 (MQ=255)
          cGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTg                                                                                                                                        <  1:3001571/149‑1 (MQ=255)
           gTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGt                                                                                                                                       >  1:1082031/1‑149 (MQ=255)
            tCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTc                                                                                                                                      <  2:2107485/149‑1 (MQ=255)
             cGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTcc                                                                                                                                     <  1:1855090/149‑1 (MQ=255)
              gACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCg                                                                                                                                    >  1:3095461/1‑149 (MQ=255)
                  gCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGc                                                                                                                                >  1:2460687/1‑149 (MQ=255)
                       gCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTAt                                                                                                                           <  1:2092261/149‑1 (MQ=255)
                        cGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTAtt                                                                                                                          <  2:1469301/149‑1 (MQ=255)
                          gTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGa                                                                                                                        >  1:1682991/1‑149 (MQ=255)
                          gTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGa                                                                                                                        >  1:2838574/1‑149 (MQ=255)
                          gTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGa                                                                                                                        >  1:1594236/1‑149 (MQ=255)
                             gCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGTTGCTGAATGGTGCCAGCGATCTACCGGATGAATTTGTCCGTCGCACTATTGAACt                                                                                                                     >  1:943194/1‑149 (MQ=255)
                             gCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACt                                                                                                                     <  1:3667377/149‑1 (MQ=255)
                             gCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACt                                                                                                                     <  1:83346/149‑1 (MQ=255)
                               gATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTg                                                                                                                   >  1:1491675/1‑149 (MQ=255)
                               gATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTg                                                                                                                   <  1:95495/149‑1 (MQ=255)
                                aTTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGATAGTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTgg                                                                                                                  <  2:4046641/149‑1 (MQ=255)
                                aTTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTgg                                                                                                                  >  2:2766123/1‑149 (MQ=255)
                                aTTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTgg                                                                                                                  >  1:3734155/1‑149 (MQ=255)
                                aTTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTgg                                                                                                                  >  2:3485999/1‑149 (MQ=255)
                                     tACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTca                                                                                                             >  2:1965166/1‑149 (MQ=255)
                                      aCGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTcac                                                                                                            >  2:3024450/1‑149 (MQ=255)
                                         gcgcACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACaaa                                                                                                         >  1:1464208/1‑149 (MQ=255)
                                         gcgcACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACaaa                                                                                                         >  1:3752845/1‑149 (MQ=255)
                                             aCGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTg                                                                                                     <  1:1213767/149‑1 (MQ=255)
                                               ggCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGaa                                                                                                   >  2:2287138/1‑149 (MQ=255)
                                                gCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAAc                                                                                                  <  2:259582/149‑1 (MQ=255)
                                                  ttATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGt                                                                                                >  1:906344/1‑149 (MQ=255)
                                                  ttATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGt                                                                                                >  2:259357/1‑149 (MQ=255)
                                                  ttATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGt                                                                                                <  2:2424833/149‑1 (MQ=255)
                                                   tataCAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGtt                                                                                               <  2:1228063/149‑1 (MQ=255)
                                                   tataCAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGtt                                                                                               >  1:1374788/1‑149 (MQ=255)
                                                        aGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCcac                                                                                          <  1:3286709/149‑1 (MQ=255)
                                                         gCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCcaca                                                                                         <  2:1746104/149‑1 (MQ=255)
                                                           aaaaCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGa                                                                                       <  2:1310375/149‑1 (MQ=255)
                                                           aaaaCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGa                                                                                       >  1:2186725/1‑149 (MQ=255)
                                                            aaaCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGaa                                                                                      <  2:1044548/149‑1 (MQ=255)
                                                             aaCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAAt                                                                                     >  1:3956886/1‑149 (MQ=255)
                                                             aaCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAAt                                                                                     >  2:562566/1‑149 (MQ=255)
                                                              aCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAAtt                                                                                    >  2:3839595/1‑149 (MQ=255)
                                                              aCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAAtt                                                                                    >  2:3073110/1‑149 (MQ=255)
                                                              aCGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAAtt                                                                                    >  2:3248774/1‑149 (MQ=255)
                                                                gCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGTTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCTCACTATTGAACTTGGCGTCACAAAAGTGAACTTTGCCACAGAATTaa                                                                                  >  1:3512143/1‑149 (MQ=255)
                                                                  cGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTaaaa                                                                                <  2:3149626/149‑1 (MQ=255)
                                                                   gAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGAGGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTaaaaa                                                                               <  2:419151/149‑1 (MQ=255)
                                                                    aaGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAAt                                                                              <  2:3253878/149‑1 (MQ=255)
                                                                      gattgattTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAg                                                                            <  1:2607655/149‑1 (MQ=255)
                                                                       attgattTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGc                                                                           >  2:687458/1‑149 (MQ=255)
                                                                          gattTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCtt                                                                        >  2:477032/1‑149 (MQ=255)
                                                                          gattTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCtt                                                                        >  1:3144933/1‑149 (MQ=255)
                                                                           attTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTc                                                                       <  1:2292990/149‑1 (MQ=255)
                                                                           agttCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTc                                                                       <  2:2158070/147‑1 (MQ=255)
                                                                            tttCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCg                                                                      >  1:3167891/1‑149 (MQ=255)
                                                                            tttCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGATGCCCCAGAATTAAAAATAGCCTTCg                                                                      >  1:3168162/1‑149 (MQ=255)
                                                                                cAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTgg                                                                  <  1:2980181/149‑1 (MQ=255)
                                                                                 aGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGgc                                                                 >  1:974909/1‑149 (MQ=255)
                                                                                     gCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCgg                                                             <  2:2460687/149‑1 (MQ=255)
                                                                                      cTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGTATTAAAAATAGCCTTCGCTGGCGCGGt                                                            <  2:905253/149‑1 (MQ=255)
                                                                                          cGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTaaa                                                        <  1:120305/149‑1 (MQ=255)
                                                                                           ggAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAg                                                       >  1:3079310/1‑149 (MQ=255)
                                                                                            gAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGc                                                      <  2:3242113/149‑1 (MQ=255)
                                                                                              aaTTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCt                                                    <  2:355728/149‑1 (MQ=255)
                                                                                              aaTTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCt                                                    <  1:3524027/149‑1 (MQ=255)
                                                                                               aTTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTg                                                   >  2:2099007/1‑149 (MQ=255)
                                                                                               aTTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTg                                                   >  2:20640/1‑149 (MQ=255)
                                                                                                ttCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTgg                                                  >  2:1955337/1‑149 (MQ=255)
                                                                                                     gAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGc                                             >  1:3510266/1‑149 (MQ=255)
                                                                                                     gAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGc                                             >  1:1878169/1‑149 (MQ=255)
                                                                                                       aGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCgg                                           >  1:3847653/1‑149 (MQ=255)
                                                                                                           gtggATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGaaaa                                       <  1:3959191/149‑1 (MQ=255)
                                                                                                           gtggATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGaaaa                                       <  1:757988/149‑1 (MQ=255)
                                                                                                           gtggATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGaaaa                                       >  1:4028344/1‑149 (MQ=255)
                                                                                                              gATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATcc                                    >  1:2186207/1‑149 (MQ=255)
                                                                                                                tGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGc                                  >  1:2748731/1‑149 (MQ=255)
                                                                                                                 gTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCa                                 >  2:817025/1‑149 (MQ=255)
                                                                                                                   tCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAgg                               <  1:669208/149‑1 (MQ=255)
                                                                                                                     ctctGGTGCTGCATGGTGACAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTTAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGt                             <  1:56446/149‑1 (MQ=255)
                                                                                                                       ctGGTTCTGCATGGTGCAAGCGATGTTGAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGATCGTTGCAACAGAATTAAAAATAGCCTTCGATGGCGCGGTTATAGCCTGGTTTGCGGAAAATCCGCATGGTaa                           <  1:4056062/149‑1 (MQ=255)
                                                                                                                       ctGGTGCTGCATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTaa                           <  1:3578684/149‑1 (MQ=255)
                                                                                                                          gTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGa                        <  2:3034638/149‑1 (MQ=255)
                                                                                                                            gctgcATGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATc                      <  2:1561032/149‑1 (MQ=255)
                                                                                                                               gcATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTc                   >  1:3036159/1‑149 (MQ=255)
                                                                                                                                  tGGTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGtt                >  2:3784698/1‑149 (MQ=255)
                                                                                                                                   ggTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGtta               <  2:3752829/149‑1 (MQ=255)
                                                                                                                                   ggTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGtta               <  2:2838574/149‑1 (MQ=255)
                                                                                                                                    gTGCCAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGttat              >  1:1059888/1‑149 (MQ=255)
                                                                                                                                      gCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGttatta            >  1:3373679/1‑149 (MQ=255)
                                                                                                                                       ccAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTtattat           >  2:4000142/1‑149 (MQ=255)
                                                                                                                                        cAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATtata          >  2:109882/1‑149 (MQ=255)
                                                                                                                                        cAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATtata          >  2:1663183/1‑149 (MQ=255)
                                                                                                                                        cAGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATtata          >  2:2186715/1‑149 (MQ=255)
                                                                                                                                         aGCGATGTTCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATTatat         <  2:386051/149‑1 (MQ=255)
                                                                                                                                          gCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATTATATg        <  2:2162859/149‑1 (MQ=255)
                                                                                                                                                ttCAGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATTATATGCGCGTc  <  1:2613732/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AACTGACTGGCGTCGACAGCCTGGCGGTAGCGATTGGTACGGCGCACGGCTTATACAGCAAAACGCCGAAGATTGATTTCCAGCGGCTGGCGGAAATTCGTGAAGTGGTGGATGTTCCTCTGGTGCTGCATGGTGCCAGCGATGTTCCGGATGAATTTGTCCGTCGCACTATTGAACTTGGCGTCACAAAAGTGAACGTTGCCACAGAATTAAAAATAGCCTTCGCTGGCGCGGTTAAAGCCTGGTTTGCGGAAAATCCGCAGGGTAATGATCCTCGTTATTATATGCGCGTC  >  CP009273/3276995‑3277287

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: