Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 557,784 G→T 7.9% 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,7840GT 7.9% 1195.8 / 50.1 316intergenic (+28/‑8)sfmD/sfmHputative outer membrane export usher protein/FimA homolog
Reads supporting (aligned to +/- strand):  ref base G (147/144);  new base T (12/13);  total (159/157)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 319 total aligned reads displayed.

TACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTT  >  CP009273/557635‑557933
                                                                                                                                                     |                                                                                                                                                     
tACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAg                                                                                                                                                       <  2:757904/150‑1 (MQ=255)
     gTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGattat                                                                                                                                                  >  1:2551604/1‑150 (MQ=255)
       tCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGattatta                                                                                                                                                <  1:3280576/150‑1 (MQ=255)
       tCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGattatta                                                                                                                                                <  1:1867037/150‑1 (MQ=255)
          gTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTg                                                                                                                                             <  2:1863045/150‑1 (MQ=255)
           tGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTgg                                                                                                                                            >  2:2335757/1‑150 (MQ=255)
            ggggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGc                                                                                                                                           >  2:138380/1‑150 (MQ=255)
            ggggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGc                                                                                                                                           <  1:3136715/150‑1 (MQ=255)
            ggggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGc                                                                                                                                           >  2:4656904/1‑150 (MQ=255)
             gggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCg                                                                                                                                          >  2:2654766/1‑150 (MQ=255)
              ggAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGa                                                                                                                                         <  2:3026042/150‑1 (MQ=255)
               gagaCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGAt                                                                                                                                        <  1:4947108/150‑1 (MQ=255)
               gagaCGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGAt                                                                                                                                        >  2:491790/1‑150 (MQ=255)
                    ggCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGAAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCAt                                                                                                                                   >  2:1937547/1‑150 (MQ=255)
                    ggCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGAAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCAt                                                                                                                                   >  2:3152698/1‑150 (MQ=255)
                    ggCGCGAACTCACGCTGCATCGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCAt                                                                                                                                   <  2:3460634/150‑1 (MQ=255)
                        cgAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGttt                                                                                                                               >  1:3248949/1‑150 (MQ=255)
                        cgAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGttt                                                                                                                               >  2:834198/1‑150 (MQ=255)
                        cgAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGttt                                                                                                                               <  2:4303173/150‑1 (MQ=255)
                        cgAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGttt                                                                                                                               >  2:3334300/1‑150 (MQ=255)
                         gAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTtg                                                                                                                              >  1:4442835/1‑150 (MQ=255)
                          aaCTCACGCTGCATTTCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCAGTCACTGTTATTTAGGCAGTTTGCACACATCCTGGCTCATAAAGGTAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTtgt                                                                                                                             <  2:4466914/150‑1 (MQ=255)
                          aaCTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTtgt                                                                                                                             >  1:4807289/1‑150 (MQ=255)
                           aCTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTtgtg                                                                                                                            <  1:716135/150‑1 (MQ=255)
                            cTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTgtgt                                                                                                                           <  1:924500/150‑1 (MQ=255)
                             tCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTc                                                                                                                          <  2:4383456/150‑1 (MQ=255)
                               aCGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTg                                                                                                                        <  1:3610425/150‑1 (MQ=255)
                                   tGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCaa                                                                                                                    >  1:4014457/1‑150 (MQ=255)
                                      aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                                 <  2:2734553/150‑1 (MQ=255)
                                      aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                                 >  1:3117274/1‑150 (MQ=255)
                                      aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                                 >  2:2632843/1‑150 (MQ=255)
                                      aTTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACa                                                                                                                 >  1:1384223/1‑150 (MQ=255)
                                       ttGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACat                                                                                                                >  2:582228/1‑150 (MQ=255)
                                       ttGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACat                                                                                                                <  1:1078166/150‑1 (MQ=255)
                                        tGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACata                                                                                                               >  1:1604225/1‑150 (MQ=255)
                                           gcACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATcc                                                                                                            <  2:1557722/150‑1 (MQ=255)
                                             aCTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTg                                                                                                          >  1:3741781/1‑150 (MQ=255)
                                             aCTATGTAGTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAATGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACTTATCCTg                                                                                                          <  2:4020271/150‑1 (MQ=255)
                                                aTGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGc                                                                                                       <  1:4727368/150‑1 (MQ=255)
                                                  gTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTa                                                                                                     >  2:3045080/1‑150 (MQ=255)
                                                   tATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTAc                                                                                                    <  1:4103013/150‑1 (MQ=255)
                                                     ttGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAACGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTg                                                                                                  <  1:3322448/150‑1 (MQ=255)
                                                      tGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGt                                                                                                 <  2:1741683/150‑1 (MQ=255)
                                                       gCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGtt                                                                                                >  2:2248417/1‑150 (MQ=255)
                                                          gAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTtgt                                                                                             <  1:2086540/150‑1 (MQ=255)
                                                             gCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCa                                                                                          <  1:2173275/150‑1 (MQ=255)
                                                               aaaGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCaaa                                                                                        >  1:802540/1‑150 (MQ=255)
                                                                  gCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAg                                                                                     <  1:3439916/150‑1 (MQ=255)
                                                                    ttACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATTCTGGGCTACTGTTTGTGCAAATAGTa                                                                                   >  2:1635160/1‑150 (MQ=255)
                                                                    ttACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTa                                                                                   >  1:4926125/1‑150 (MQ=255)
                                                                     tACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTAc                                                                                  <  2:472726/150‑1 (MQ=255)
                                                                        aGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATATAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTgg                                                                               <  2:1547704/150‑1 (MQ=255)
                                                                          cAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCg                                                                             >  1:2450150/1‑150 (MQ=255)
                                                                           aaGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGt                                                                            <  1:1722733/150‑1 (MQ=255)
                                                                            aGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTa                                                                           >  1:1458345/1‑150 (MQ=255)
                                                                            aGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTa                                                                           >  1:3878050/1‑150 (MQ=255)
                                                                             gCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAg                                                                          <  1:1315971/150‑1 (MQ=255)
                                                                             gCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAg                                                                          <  1:1460604/150‑1 (MQ=255)
                                                                              ccGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGc                                                                         >  1:152910/1‑150 (MQ=255)
                                                                               cGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCa                                                                        <  2:1808194/150‑1 (MQ=255)
                                                                                     tGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGAt                                                                  <  2:542392/150‑1 (MQ=255)
                                                                                        tATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGaa                                                               >  2:4662947/1‑150 (MQ=255)
                                                                                          tttCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAaca                                                             >  2:1333281/1‑150 (MQ=255)
                                                                                          tttCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAaca                                                             <  1:4615611/150‑1 (MQ=255)
                                                                                          tgtCGGCAGTTTGCACTCATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGTTTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTTTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAaca                                                             <  2:3230077/148‑1 (MQ=255)
                                                                                           ttCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAacac                                                            <  1:4029269/150‑1 (MQ=255)
                                                                                               gCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATg                                                        <  1:2272281/150‑1 (MQ=255)
                                                                                               gCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATg                                                        <  2:1689152/150‑1 (MQ=255)
                                                                                                cAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGa                                                       >  2:1445902/1‑150 (MQ=255)
                                                                                                 aGTTTGCACACATCCTGGCTCATAAAGGACATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAt                                                      >  2:3383068/1‑150 (MQ=255)
                                                                                                    ttGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctc                                                   >  1:4632548/1‑150 (MQ=255)
                                                                                                     tGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCATATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctct                                                  >  1:2732968/1‑150 (MQ=255)
                                                                                                      gCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGAtctctc                                                 >  2:1056563/1‑150 (MQ=255)
                                                                                                        acacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCaa                                               >  2:2048809/1‑150 (MQ=255)
                                                                                                        acacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCaa                                               <  1:2034531/150‑1 (MQ=255)
                                                                                                         cacaTCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCaaa                                              <  1:2654769/150‑1 (MQ=255)
                                                                                                               cTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCtt                                        <  1:2642207/150‑1 (MQ=255)
                                                                                                                tGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCttt                                       <  1:1015113/150‑1 (MQ=255)
                                                                                                                 ggCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTa                                      <  2:689585/150‑1 (MQ=255)
                                                                                                                  gCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTaa                                     >  1:4187396/1‑150 (MQ=255)
                                                                                                                    tCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATTCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATa                                   <  2:3797905/150‑1 (MQ=255)
                                                                                                                       tAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCa                                >  2:158113/1‑150 (MQ=255)
                                                                                                                        aaaGTAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAc                               <  2:4831555/150‑1 (MQ=255)
                                                                                                                        aaaGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCAc                               >  2:4174117/1‑150 (MQ=255)
                                                                                                                          aGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATTCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCa                             >  1:4677655/1‑150 (MQ=255)
                                                                                                                          aGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCa                             >  2:114966/1‑150 (MQ=255)
                                                                                                                            gAAATTATCAATAAGATAATCTGCAGATTATTTTTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCATATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAAt                           <  1:666847/150‑1 (MQ=255)
                                                                                                                             naaTTATCAATAAGATACTCTGCAGATTATTATTGGCGATGGCATGGTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATa                          >  2:2480307/2‑150 (MQ=255)
                                                                                                                                 tATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAcca                      >  1:125763/1‑150 (MQ=255)
                                                                                                                                  aTCAATAAGATACTCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccag                     >  2:3720073/1‑150 (MQ=255)
                                                                                                                                   tCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagc                    >  1:2883458/1‑150 (MQ=255)
                                                                                                                                      aTAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagccag                 >  2:2574138/1‑150 (MQ=255)
                                                                                                                                        aaGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAggg               <  1:4224371/150‑1 (MQ=255)
                                                                                                                                         aGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGc              >  2:1789338/1‑150 (MQ=255)
                                                                                                                                            tAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCTAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGa           >  1:1004894/1‑150 (MQ=255)
                                                                                                                                            tAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGa           >  2:2840054/1‑150 (MQ=255)
                                                                                                                                              aTCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAtt         >  2:4480169/1‑150 (MQ=255)
                                                                                                                                                 tGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgtt      <  1:4606420/150‑1 (MQ=255)
                                                                                                                                                    aGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTgttgtt   >  1:985180/1‑150 (MQ=255)
                                                                                                                                                     gATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGttt  >  1:569355/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACTGGTTCAGTGGGGAGACGGCGCGAACTCACGCTGCATTGCGCACTATGTATTGCCGAAGCAAAGCTTACAGCAAGCCGTCACTGTTATTTCGGCAGTTTGCACACATCCTGGCTCATAAAGGAAATTATCAATAAGATAATCTGCAGATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTT  >  CP009273/557635‑557933

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: