Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 389,725 G→T 11.8% pseudogene (1300/1449 nt) yaiT → pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273389,7250GT11.8% 658.6 / 52.2 200pseudogene (1300/1449 nt)yaiTpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base G (91/85);  new base T (11/13);  total (102/98)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

ACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGT  >  CP009273/389576‑389872
                                                                                                                                                     |                                                                                                                                                   
aCGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACg                                                                                                                                                     >  2:2433491/1‑150 (MQ=255)
  gACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGt                                                                                                                                                   >  2:4888337/1‑150 (MQ=255)
    cAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTAc                                                                                                                                                 >  1:770300/1‑150 (MQ=255)
     aGCATGCGTTATGAACTGGGTGTAGATGCAGGTTAGACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACt                                                                                                                                                <  2:2969494/150‑1 (MQ=255)
         tGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAg                                                                                                                                            <  2:4020029/150‑1 (MQ=255)
          gCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAgg                                                                                                                                           <  1:3144632/150‑1 (MQ=255)
             ttATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTc                                                                                                                                        >  1:459537/1‑150 (MQ=255)
               aTGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTg                                                                                                                                      >  2:3946716/1‑150 (MQ=255)
                 gAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                                                                                                                                    >  1:4912384/1‑150 (MQ=255)
                 gAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                                                                                                                                    >  1:2330825/1‑150 (MQ=255)
                 gAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                                                                                                                                    >  2:493644/1‑150 (MQ=255)
                    cTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTa                                                                                                                                 >  2:1535792/1‑150 (MQ=255)
                    cTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTa                                                                                                                                 <  2:1911944/150‑1 (MQ=255)
                    cTGGGTGTAGATCCAGGTTATACCATCACCAACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAGCAACGATAACGAGGTGAACGGCGATTCCATCGATAACTGTCCTGGAGGGTCTGCGGtg                                                                                                                                 >  2:1535660/1‑149 (MQ=255)
                            aGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTgg                                                                                                                         >  1:2849967/1‑150 (MQ=255)
                             gATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGt                                                                                                                        >  1:4863368/1‑150 (MQ=255)
                               tGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCt                                                                                                                      <  2:233167/150‑1 (MQ=255)
                                  aGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTggg                                                                                                                   <  1:4936438/150‑1 (MQ=255)
                                   ggTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGt                                                                                                                  >  1:3221247/1‑150 (MQ=255)
                                   ggTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGt                                                                                                                  >  2:1706471/1‑150 (MQ=255)
                                    gTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTa                                                                                                                 >  2:100321/1‑150 (MQ=255)
                                      tataCCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACt                                                                                                               >  1:3882970/1‑150 (MQ=255)
                                        taCCTTCACCTACCGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCa                                                                                                             <  2:1335600/150‑1 (MQ=255)
                                        taCCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCa                                                                                                             >  2:1259575/1‑150 (MQ=255)
                                         aCCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAg                                                                                                            <  2:736357/150‑1 (MQ=255)
                                          ccTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGt                                                                                                           <  2:3224374/150‑1 (MQ=255)
                                          ccTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGt                                                                                                           <  2:3224518/150‑1 (MQ=255)
                                           cTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGtt                                                                                                          >  2:2576359/1‑150 (MQ=255)
                                              cACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAg                                                                                                       <  1:5056658/150‑1 (MQ=255)
                                               aCCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGc                                                                                                      >  2:1961600/1‑150 (MQ=255)
                                                  tACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTc                                                                                                   <  1:2735941/150‑1 (MQ=255)
                                                    cAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCAc                                                                                                 <  2:838037/150‑1 (MQ=255)
                                                    cAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACGGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCAc                                                                                                 <  1:1488262/150‑1 (MQ=255)
                                                     aGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCAcc                                                                                                <  2:4028973/150‑1 (MQ=255)
                                                       cGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCaa                                                                                              >  2:157507/1‑150 (MQ=255)
                                                        gAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAg                                                                                             >  1:1454768/1‑150 (MQ=255)
                                                        gAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAg                                                                                             >  1:3107374/1‑150 (MQ=255)
                                                         aaGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGa                                                                                            >  2:3187456/1‑150 (MQ=255)
                                                          aGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGaa                                                                                           >  2:3603687/1‑150 (MQ=255)
                                                           gATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCTCCAAGAAc                                                                                          >  2:4680949/1‑150 (MQ=255)
                                                           gATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAAc                                                                                          >  1:631223/1‑150 (MQ=255)
                                                           gATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAAc                                                                                          >  2:2298467/1‑150 (MQ=255)
                                                               aGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCa                                                                                      >  1:2141978/1‑150 (MQ=255)
                                                                ggCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAg                                                                                     <  1:3666447/150‑1 (MQ=255)
                                                                ggCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAg                                                                                     <  1:2048608/150‑1 (MQ=255)
                                                                     tGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCt                                                                                <  2:1107996/150‑1 (MQ=255)
                                                                     tGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCt                                                                                <  1:1089552/150‑1 (MQ=255)
                                                                      gACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCta                                                                               >  2:3513102/1‑150 (MQ=255)
                                                                           cGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCg                                                                          <  1:5105491/150‑1 (MQ=255)
                                                                              aCTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATg                                                                       >  1:1966100/1‑150 (MQ=255)
                                                                              aCTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATg                                                                       >  1:3483545/1‑150 (MQ=255)
                                                                               cTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGc                                                                      >  1:2574625/1‑150 (MQ=255)
                                                                                  cAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTaa                                                                   <  2:2573399/150‑1 (MQ=255)
                                                                                  cAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTaa                                                                   >  2:1267364/1‑150 (MQ=255)
                                                                                  cAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTaa                                                                   <  1:378579/150‑1 (MQ=255)
                                                                                   aaaCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAAc                                                                  <  1:3506426/150‑1 (MQ=255)
                                                                                   aaaCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAAc                                                                  <  2:3475856/150‑1 (MQ=255)
                                                                                    aaCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACt                                                                 <  1:5184669/150‑1 (MQ=255)
                                                                                     aCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTa                                                                <  1:329372/150‑1 (MQ=255)
                                                                                         gCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTc                                                            <  2:2709976/150‑1 (MQ=255)
                                                                                         gCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTc                                                            <  2:2918614/150‑1 (MQ=255)
                                                                                           ttACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCgg                                                          >  1:4269650/1‑150 (MQ=255)
                                                                                           ttACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCgg                                                          <  1:4303940/150‑1 (MQ=255)
                                                                                            tACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCggt                                                         >  1:2220165/1‑150 (MQ=255)
                                                                                                tCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCggtggtg                                                     <  2:770608/150‑1 (MQ=255)
                                                                                                tCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCggtggtg                                                     <  1:2866863/150‑1 (MQ=255)
                                                                                                   acgacgACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGgtggtggtg                                                  <  1:3367912/150‑1 (MQ=255)
                                                                                                          ctctAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGa                                           >  1:243182/1‑150 (MQ=255)
                                                                                                          ctctAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGCTGGTGGTGACGTAGa                                           >  2:5172545/1‑150 (MQ=255)
                                                                                                            ctAACAACGATAACGCTGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATc                                         <  1:4275987/150‑1 (MQ=255)
                                                                                                            ctAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATc                                         >  1:2387562/1‑150 (MQ=255)
                                                                                                              aacaacGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCaa                                       >  1:439732/1‑150 (MQ=255)
                                                                                                              aacaacGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCaa                                       >  1:120761/1‑150 (MQ=255)
                                                                                                                caacGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGa                                     >  2:1319020/1‑150 (MQ=255)
                                                                                                                 aacgataacgatGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGAc                                    <  1:4323779/150‑1 (MQ=255)
                                                                                                                 aacgataacgatGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGAc                                    >  2:4762246/1‑150 (MQ=255)
                                                                                                                 aacgataacgatGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGAc                                    <  1:3581438/150‑1 (MQ=255)
                                                                                                                 aacgataacgatGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGAc                                    <  1:4554988/150‑1 (MQ=255)
                                                                                                                    gataacgatGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTgg                                 <  1:1655107/150‑1 (MQ=255)
                                                                                                                       aacgatGTGAACGGCGATTCCATCGATAACTGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTcc                              <  1:1078189/150‑1 (MQ=255)
                                                                                                                             gtgAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAAc                        <  2:459292/150‑1 (MQ=255)
                                                                                                                              tgAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACg                       <  1:4732088/150‑1 (MQ=255)
                                                                                                                              tgAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACg                       <  1:1582910/150‑1 (MQ=255)
                                                                                                                                    gCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGgtg                 <  2:3899580/150‑1 (MQ=255)
                                                                                                                                    gCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGgtg                 <  1:4675082/150‑1 (MQ=255)
                                                                                                                                     cGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGgtgt                <  1:699106/150‑1 (MQ=255)
                                                                                                                                     cGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGgtgt                <  1:698972/150‑1 (MQ=255)
                                                                                                                                      gATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGtt               >  2:4665114/1‑150 (MQ=255)
                                                                                                                                      gATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGtt               <  1:3674194/150‑1 (MQ=255)
                                                                                                                                       aTTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTa              <  1:3905311/150‑1 (MQ=255)
                                                                                                                                          ccATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAat           >  1:960038/1‑150 (MQ=255)
                                                                                                                                            aTCGATAACGGAACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTAATCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGAAGATCAAGACTGGACCGCGAACGAGGGTGTTAAatat         >  2:4146305/1‑150 (MQ=255)
                                                                                                                                             tCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAAtata        >  1:2970778/1‑150 (MQ=255)
                                                                                                                                               gATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATAcc      <  1:2433491/150‑1 (MQ=255)
                                                                                                                                               gATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTCCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATAcc      >  2:4182717/1‑150 (MQ=255)
                                                                                                                                                aTAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCt     <  1:1871553/150‑1 (MQ=255)
                                                                                                                                                aTAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCt     >  2:951890/1‑150 (MQ=255)
                                                                                                                                                  aaCGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTgg   >  2:2870330/1‑150 (MQ=255)
                                                                                                                                                  aaCGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTgg   >  2:1293406/1‑150 (MQ=255)
                                                                                                                                                   aCGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGt  <  1:1508577/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGGGTACTCAGTTTAGCTTCACCAAGAACTTCAGCGCCTATACCGATGCTAACTACCTCGGTGGTGGTGACGTAGATCAAGACTGGTCCGCGAACGTGGGTGTTAAATATACCTGGT  >  CP009273/389576‑389872

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: