Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,887 A→T 9.4% S32S (TCA→TCT)  sfmH → FimA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273557,8870AT 9.4% 819.4 / 59.4 255S32S (TCA→TCT) sfmHFimA homolog
Reads supporting (aligned to +/- strand):  ref base A (126/105);  new base T (15/9);  total (141/114)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 257 total aligned reads displayed.

CACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGA  >  CP009273/557740‑558036
                                                                                                                                                   |                                                                                                                                                     
 acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                                                                                                                                    >  2:859533/1‑150 (MQ=255)
 acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                                                                                                                                    <  1:3830961/150‑1 (MQ=255)
 acaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAat                                                                                                                                                    <  1:3832176/150‑1 (MQ=255)
     ccTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCt                                                                                                                                                >  2:3518687/1‑150 (MQ=255)
        ggCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTa                                                                                                                                             >  2:425492/1‑150 (MQ=255)
         gCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                                                                                                                            >  1:1821711/1‑150 (MQ=255)
         gCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                                                                                                                            >  1:2742444/1‑150 (MQ=255)
         gCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                                                                                                                            <  2:2023135/150‑1 (MQ=255)
           tCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATa                                                                                                                                          <  1:3255931/150‑1 (MQ=255)
           tCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATa                                                                                                                                          <  1:2001864/150‑1 (MQ=255)
            cATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAg                                                                                                                                         <  1:3810197/150‑1 (MQ=255)
            cATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAg                                                                                                                                         >  1:3411017/1‑150 (MQ=255)
             aTAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGc                                                                                                                                        <  2:1014867/150‑1 (MQ=255)
               aaaGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAc                                                                                                                                      <  2:4500219/150‑1 (MQ=255)
                aaGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCTGAAGATGAACACTATGATCTCTCAATTATCTTTAATAGCAcc                                                                                                                                     <  1:4432973/150‑1 (MQ=255)
                aaGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAcc                                                                                                                                     >  2:3630877/1‑150 (MQ=255)
                  ggAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCaa                                                                                                                                   >  2:453348/1‑150 (MQ=255)
                   gAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAAt                                                                                                                                  <  2:4205316/150‑1 (MQ=255)
                    aaaTTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATa                                                                                                                                 <  1:526700/150‑1 (MQ=255)
                     aaTTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATaa                                                                                                                                >  2:3325569/1‑150 (MQ=255)
                     aaTTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATaa                                                                                                                                >  2:1963600/1‑150 (MQ=255)
                      aTTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAc                                                                                                                               >  2:3048421/1‑150 (MQ=255)
                      aTTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAAc                                                                                                                               >  2:1513923/1‑150 (MQ=255)
                       ttATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcc                                                                                                                              >  2:3836869/1‑150 (MQ=255)
                        tATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcca                                                                                                                             >  2:725212/1‑150 (MQ=255)
                         aTCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccag                                                                                                                            >  1:337313/1‑150 (MQ=255)
                         aTCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccag                                                                                                                            <  2:2845644/150‑1 (MQ=255)
                          tCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagc                                                                                                                           >  2:1844692/1‑150 (MQ=255)
                             aTAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAAccagccag                                                                                                                        >  2:1192173/1‑150 (MQ=255)
                               aaGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAggg                                                                                                                      >  1:2790890/1‑150 (MQ=255)
                                aGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGc                                                                                                                     >  1:4002596/1‑150 (MQ=255)
                                 gATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCa                                                                                                                    <  1:1935771/150‑1 (MQ=255)
                                  aTAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAg                                                                                                                   >  2:1678416/1‑150 (MQ=255)
                                  aTAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAg                                                                                                                   >  2:1678564/1‑150 (MQ=255)
                                   tAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGa                                                                                                                  >  1:2559479/1‑150 (MQ=255)
                                       cTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgt                                                                                                              >  2:190788/1‑150 (MQ=255)
                                        tGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgtt                                                                                                             >  1:2601587/1‑150 (MQ=255)
                                        tGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgtt                                                                                                             <  2:1614844/150‑1 (MQ=255)
                                           atattattattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTgttgtt                                                                                                          >  2:1025385/3‑150 (MQ=255)
                                           atattattattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTgttgtt                                                                                                          >  1:3230134/3‑150 (MQ=255)
                                             attattattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGtttt                                                                                                        <  2:580165/150‑1 (MQ=255)
                                                    ttGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGaa                                                                                                 <  1:3510303/150‑1 (MQ=255)
                                                    ttGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGaa                                                                                                 >  1:3115381/1‑150 (MQ=255)
                                                    ttGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTACCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGaa                                                                                                 <  2:4397955/150‑1 (MQ=255)
                                                       gCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGaaaaa                                                                                              <  1:1269695/150‑1 (MQ=255)
                                                        cGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAAt                                                                                             <  1:64197/150‑1 (MQ=255)
                                                           tGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCg                                                                                          >  1:416245/1‑150 (MQ=255)
                                                              cATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCt                                                                                       <  2:2314486/150‑1 (MQ=255)
                                                              cATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCt                                                                                       <  1:230402/150‑1 (MQ=255)
                                                                tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTgg                                                                                     >  2:2890003/1‑150 (MQ=255)
                                                                tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTgg                                                                                     <  1:1102728/150‑1 (MQ=255)
                                                                tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATATTTAATAGCACCAATAAACAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTgg                                                                                     >  2:4334994/1‑150 (MQ=255)
                                                                        tCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGgtgt                                                                             >  1:709550/1‑150 (MQ=255)
                                                                         cTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGtc                                                                            <  2:1280653/150‑1 (MQ=255)
                                                                            gCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCa                                                                         >  1:165944/1‑150 (MQ=255)
                                                                             cAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAg                                                                        <  1:102920/150‑1 (MQ=255)
                                                                                  atatCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAAtt                                                                   <  2:2956229/150‑1 (MQ=255)
                                                                                      ccTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTc                                                               <  2:3772736/150‑1 (MQ=255)
                                                                                        tGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCa                                                             >  1:1872760/1‑150 (MQ=255)
                                                                                           gCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCAccc                                                          >  2:1128006/1‑150 (MQ=255)
                                                                                              aCTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGt                                                       <  1:3488052/150‑1 (MQ=255)
                                                                                              aCTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGt                                                       >  2:3191045/1‑150 (MQ=255)
                                                                                               cTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTa                                                      <  2:1903737/150‑1 (MQ=255)
                                                                                                tGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTAc                                                     <  2:4017633/150‑1 (MQ=255)
                                                                                                tGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTAc                                                     <  1:1310139/150‑1 (MQ=255)
                                                                                                  tttGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGc                                                   >  2:2189740/1‑150 (MQ=255)
                                                                                                   ttGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCt                                                  <  1:3048421/150‑1 (MQ=255)
                                                                                                    tgtgCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTg                                                 >  1:653717/1‑150 (MQ=255)
                                                                                                       gCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTg                                              <  1:2081038/150‑1 (MQ=255)
                                                                                                         aaaTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGaa                                            <  1:2987691/150‑1 (MQ=255)
                                                                                                         aaaTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGaa                                            >  2:2141508/1‑150 (MQ=255)
                                                                                                           aTAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAAtt                                          <  2:3230134/150‑1 (MQ=255)
                                                                                                             aGTACTGGCGTAGCAGAAGATGAACACAATGATCTCTCAAATATCTTTCATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATtat                                        <  2:3674343/150‑1 (MQ=255)
                                                                                                              gTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATtata                                       <  1:329217/150‑1 (MQ=255)
                                                                                                              gTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATtata                                       <  1:2625508/150‑1 (MQ=255)
                                                                                                                  tGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTatacata                                   >  1:1022161/1‑150 (MQ=255)
                                                                                                                   ggCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTatacatac                                  <  2:1821711/150‑1 (MQ=255)
                                                                                                                     cGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCg                                <  1:1700290/150‑1 (MQ=255)
                                                                                                                          cAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGtt                           <  2:85329/150‑1 (MQ=255)
                                                                                                                          cAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGtt                           <  2:4372316/150‑1 (MQ=255)
                                                                                                                           agaagaTGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTa                          >  2:950121/1‑150 (MQ=255)
                                                                                                                             aagaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATg                        >  2:1336664/1‑150 (MQ=255)
                                                                                                                              agaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGt                       >  2:3421838/1‑150 (MQ=255)
                                                                                                                              agaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGt                       >  2:3872671/1‑150 (MQ=255)
                                                                                                                              agaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGt                       >  2:2571768/1‑150 (MQ=255)
                                                                                                                               gaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTc                      <  2:3391311/150‑1 (MQ=255)
                                                                                                                               gaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTc                      <  2:1056235/150‑1 (MQ=255)
                                                                                                                               gaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTc                      >  2:3168696/1‑150 (MQ=255)
                                                                                                                                aTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCa                     >  1:2853419/1‑150 (MQ=255)
                                                                                                                                aTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCa                     >  1:2852677/1‑150 (MQ=255)
                                                                                                                                  gAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCAcc                   >  2:3262549/1‑150 (MQ=255)
                                                                                                                                  gAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCAcc                   >  2:614029/1‑150 (MQ=255)
                                                                                                                                      acTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACt               <  1:148182/150‑1 (MQ=255)
                                                                                                                                         aTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTa            <  2:2559479/150‑1 (MQ=255)
                                                                                                                                          tGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTAt           >  2:1818055/1‑150 (MQ=255)
                                                                                                                                           gATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTAtt          <  1:4289841/150‑1 (MQ=255)
                                                                                                                                             tctctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGt        >  1:3831289/1‑150 (MQ=255)
                                                                                                                                             tctctcAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGt        <  2:904870/150‑1 (MQ=255)
                                                                                                                                                  cAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAgg   >  1:1380834/1‑150 (MQ=255)
                                                                                                                                                   aaaTATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAAGTTTATTGTTCAGGa  <  2:2892994/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAGCAATTTGTCCACCCGGTACGCTGGTGAATTATACATACCGTAGTTATGTCACCAACTTTATTGTTCAGGA  >  CP009273/557740‑558036

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