Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,784 G→T 50.6% intergenic (+28/‑8) sfmD → / → sfmH putative outer membrane export usher protein/FimA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,7840GT50.6% 89.4 / inf 247intergenic (+28/‑8)sfmD/sfmHputative outer membrane export usher protein/FimA homolog
Reads supporting (aligned to +/- strand):  ref base G (69/53);  new base T (78/47);  total (147/100)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

Only 100 of 250 total aligned reads displayed.

CTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTT  >  CP009273/557637‑557929
                                                                                                                                                   |                                                                                                                                                 
cTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGat                                                                                                                                                 <  1:3268586/150‑1 (MQ=255)
      cAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAtattattat                                                                                                                                           <  2:1111235/150‑1 (MQ=255)
       aGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGattattatt                                                                                                                                          <  1:4178642/150‑1 (MQ=255)
       aGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGattattatt                                                                                                                                          <  2:3995006/150‑1 (MQ=255)
        gTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTg                                                                                                                                         >  1:2023135/1‑150 (MQ=255)
         tGGGGAGACGGCGCGAACTCAAGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCTATAAGATAATCTGCATATTATTATTgg                                                                                                                                        <  2:2914652/150‑1 (MQ=255)
          ggggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGc                                                                                                                                       <  1:2310418/150‑1 (MQ=255)
              agaCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATg                                                                                                                                   >  2:21237/1‑150 (MQ=255)
               gaCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATgg                                                                                                                                  >  1:1351685/1‑150 (MQ=255)
                aCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGc                                                                                                                                 <  1:2770687/150‑1 (MQ=255)
                  ggCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCAt                                                                                                                               >  1:4101073/1‑150 (MQ=255)
                  ggCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCAt                                                                                                                               >  1:1881089/1‑150 (MQ=255)
                    agcgAACTCACGCTTCATGGCGAACTATGTATTGCCGAAGCAAAGCTTACAGCTAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATATAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGt                                                                                                                             <  2:76631/149‑1 (MQ=255)
                     gcgAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGtt                                                                                                                            <  2:2476924/150‑1 (MQ=255)
                     gcgAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGtt                                                                                                                            <  2:175940/150‑1 (MQ=255)
                         aCTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTAGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTtgtg                                                                                                                        >  2:1363503/1‑150 (MQ=255)
                          cTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTgtgt                                                                                                                       >  1:2498988/1‑150 (MQ=255)
                           tCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTc                                                                                                                      >  2:3898286/1‑150 (MQ=255)
                           tCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTc                                                                                                                      <  2:3797707/150‑1 (MQ=255)
                            cACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCt                                                                                                                     >  2:764893/1‑150 (MQ=255)
                                 tGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCaa                                                                                                                >  2:64441/1‑150 (MQ=255)
                                 tGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCaa                                                                                                                >  2:2368623/1‑150 (MQ=255)
                                   cATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAAc                                                                                                              >  1:3780762/1‑150 (MQ=255)
                                    aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                             <  2:3079294/150‑1 (MQ=255)
                                    aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                             <  2:797498/150‑1 (MQ=255)
                                      tGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACata                                                                                                           >  1:1074471/1‑150 (MQ=255)
                                         gcACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATcc                                                                                                        <  1:123350/150‑1 (MQ=255)
                                              aTGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGc                                                                                                   >  1:968576/1‑150 (MQ=255)
                                              aTGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGc                                                                                                   >  1:1250335/1‑150 (MQ=255)
                                              aTGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGc                                                                                                   >  1:1250243/1‑150 (MQ=255)
                                                           gCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCa                                                                                      <  2:4101073/150‑1 (MQ=255)
                                                             aaaGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCaaa                                                                                    >  1:4381844/1‑150 (MQ=255)
                                                                gCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAg                                                                                 <  2:1007510/150‑1 (MQ=255)
                                                                 cTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGt                                                                                >  1:2618257/1‑150 (MQ=255)
                                                                 cTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGt                                                                                <  2:2666867/150‑1 (MQ=255)
                                                                  ttACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTa                                                                               >  1:1879381/1‑150 (MQ=255)
                                                                   tACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTAc                                                                              <  1:4225273/150‑1 (MQ=255)
                                                                      aGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTgg                                                                           <  2:3615562/150‑1 (MQ=255)
                                                                         aaGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGt                                                                        >  1:3656944/1‑150 (MQ=255)
                                                                         aaGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGt                                                                        >  1:4017633/1‑150 (MQ=255)
                                                                          aGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTa                                                                       >  1:1056235/1‑150 (MQ=255)
                                                                              gTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCag                                                                   >  1:4347878/1‑150 (MQ=255)
                                                                              gTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCag                                                                   >  2:3870430/1‑150 (MQ=255)
                                                                               tCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCaga                                                                  <  2:1404126/150‑1 (MQ=255)
                                                                                  cTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCagaaga                                                               >  1:1628957/1‑150 (MQ=255)
                                                                                   tGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGAt                                                              >  1:3285618/1‑150 (MQ=255)
                                                                                      tATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGaa                                                           >  2:329217/1‑150 (MQ=255)
                                                                                         ttCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAacac                                                        <  1:934883/150‑1 (MQ=255)
                                                                                          tCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACt                                                       >  2:478654/1‑150 (MQ=255)
                                                                                           cGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTa                                                      >  2:3695931/1‑150 (MQ=255)
                                                                                           cGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTa                                                      <  2:3148402/150‑1 (MQ=255)
                                                                                           cGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTa                                                      <  2:1691627/150‑1 (MQ=255)
                                                                                             gCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATg                                                    <  1:338829/150‑1 (MQ=255)
                                                                                               aGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAt                                                  >  1:2127965/1‑150 (MQ=255)
                                                                                               aGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAt                                                  >  1:580179/1‑150 (MQ=255)
                                                                                                 tttGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtct                                                <  1:3898286/150‑1 (MQ=255)
                                                                                                   tGCACACATCCTGGCTCATATTGTAAATTATCAATAAGATAATATGCAGATTATTATTGGCGATGGCATGTTTGTGTATGTCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGTACACTATGAtctct                                              <  2:2053122/150‑1 (MQ=255)
                                                                                                   tGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctct                                              >  2:36551/1‑150 (MQ=255)
                                                                                                   tGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctct                                              <  2:2610731/150‑1 (MQ=255)
                                                                                                    gCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCTTGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctctc                                             <  1:4436633/150‑1 (MQ=255)
                                                                                                    gCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctctc                                             <  1:859478/150‑1 (MQ=255)
                                                                                                       cacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCaaa                                          <  1:764893/150‑1 (MQ=255)
                                                                                                        acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                         >  2:859533/1‑150 (MQ=255)
                                                                                                        acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                         <  1:4069722/150‑1 (MQ=255)
                                                                                                        acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                         <  1:3831078/150‑1 (MQ=255)
                                                                                                         caTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAata                                        >  2:513094/1‑150 (MQ=255)
                                                                                                            ccTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCt                                     >  2:3518687/1‑150 (MQ=255)
                                                                                                             cTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCtt                                    <  1:68518/150‑1 (MQ=255)
                                                                                                               ggCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTa                                  >  1:3403879/1‑150 (MQ=255)
                                                                                                               ggCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGACTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTa                                  >  2:4372516/1‑150 (MQ=255)
                                                                                                               ggCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTa                                  >  2:425492/1‑150 (MQ=255)
                                                                                                                gCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                 <  2:2023107/150‑1 (MQ=255)
                                                                                                                gCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                 >  1:1820456/1‑150 (MQ=255)
                                                                                                                gCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                 >  1:2742444/1‑150 (MQ=255)
                                                                                                                 cTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAAt                                <  1:4456036/150‑1 (MQ=255)
                                                                                                                  tCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATa                               <  1:3255911/150‑1 (MQ=255)
                                                                                                                  tCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATa                               <  1:2001864/150‑1 (MQ=255)
                                                                                                                   cATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAg                              <  1:3810138/150‑1 (MQ=255)
                                                                                                                   cATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAg                              <  1:3810217/150‑1 (MQ=255)
                                                                                                                    aTAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGc                             <  2:1014867/150‑1 (MQ=255)
                                                                                                                     tAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCa                            >  1:2892994/1‑150 (MQ=255)
                                                                                                                       aaGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAcc                          >  2:3630877/1‑150 (MQ=255)
                                                                                                                           aaaTTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATa                      <  1:526700/150‑1 (MQ=255)
                                                                                                                           aaaTTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATa                      >  2:2274165/1‑150 (MQ=255)
                                                                                                                            aaTTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATaa                     >  1:2845644/1‑150 (MQ=255)
                                                                                                                              ttATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcc                   >  2:3029786/1‑150 (MQ=255)
                                                                                                                              ttATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcc                   >  2:3396137/1‑150 (MQ=255)
                                                                                                                               tATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcca                  >  2:2938156/1‑150 (MQ=255)
                                                                                                                                aTCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccag                 <  2:2845644/150‑1 (MQ=255)
                                                                                                                                 tCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagc                >  2:4527028/1‑150 (MQ=255)
                                                                                                                                 tCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagc                >  2:1844692/1‑150 (MQ=255)
                                                                                                                                    aTAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAAccagccag             >  2:1192173/1‑150 (MQ=255)
                                                                                                                                       aGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGc          >  1:4002596/1‑150 (MQ=255)
                                                                                                                                         aTAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAg        >  2:1678416/1‑150 (MQ=255)
                                                                                                                                          tAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGa       >  1:2559505/1‑150 (MQ=255)
                                                                                                                                              cTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgt   >  1:3432969/1‑150 (MQ=255)
                                                                                                                                              cTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgt   >  2:190788/1‑150 (MQ=255)
                                                                                                                                              cTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgt   <  1:625937/150‑1 (MQ=255)
                                                                                                                                               tGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgtt  >  1:2601598/1‑150 (MQ=255)
                                                                                                                                               tGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgtt  <  2:1614866/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTT  >  CP009273/557637‑557929

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: