Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 729,149 G→A 6.6% D1371N (GAT→AAT)  rhsC → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273729,1490GA 6.6% 1045.7 / 42.1 304D1371N (GAT→AAT) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (154/129);  new base A (10/10);  total (164/140)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

GTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTT  >  CP009273/729000‑729298
                                                                                                                                                     |                                                                                                                                                     
gTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATg                                                                                                                                                       >  1:2508191/1‑150 (MQ=255)
   gtgACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGAtt                                                                                                                                                    <  1:135684/150‑1 (MQ=255)
       ccTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCttt                                                                                                                                                <  1:4205991/150‑1 (MQ=255)
        cTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTACGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTa                                                                                                                                               <  2:4287073/150‑1 (MQ=255)
        cTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTa                                                                                                                                               >  2:4264384/1‑150 (MQ=255)
         tCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAg                                                                                                                                              <  2:2025830/150‑1 (MQ=255)
          cGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGa                                                                                                                                             <  1:1290155/150‑1 (MQ=255)
            ttGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGaaa                                                                                                                                           >  1:4414234/1‑150 (MQ=255)
             tGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGTTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGACAAGACTGGTCTAATGATTCTTTAAAAAc                                                                                                                                          <  2:3403744/150‑1 (MQ=255)
                  ccGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTa                                                                                                                                     >  1:1199399/1‑150 (MQ=255)
                   cGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTaa                                                                                                                                    <  1:1767732/150‑1 (MQ=255)
                    gATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAg                                                                                                                                   <  1:3864025/150‑1 (MQ=255)
                      taataGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGAt                                                                                                                                 >  2:3042611/1‑150 (MQ=255)
                       aataGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATg                                                                                                                                >  1:3360744/1‑150 (MQ=255)
                       aataGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATg                                                                                                                                >  1:3465551/1‑150 (MQ=255)
                         taGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGaa                                                                                                                              >  1:200606/1‑150 (MQ=255)
                           gTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGt                                                                                                                            >  2:3331466/1‑150 (MQ=255)
                           gTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGt                                                                                                                            >  1:4145460/1‑150 (MQ=255)
                                gtgGAGAACTGTTTAAAAGTGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTcc                                                                                                                       <  2:2642758/150‑1 (MQ=255)
                                gtgGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTcc                                                                                                                       <  1:2813011/150‑1 (MQ=255)
                                 tgGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCg                                                                                                                      <  1:3833865/150‑1 (MQ=255)
                                      aaCTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAAt                                                                                                                 <  1:1692959/150‑1 (MQ=255)
                                       aCTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATa                                                                                                                <  2:479907/150‑1 (MQ=255)
                                          gTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACa                                                                                                             >  1:2234993/1‑150 (MQ=255)
                                          gTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACa                                                                                                             <  2:1014089/150‑1 (MQ=255)
                                           tttAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACAt                                                                                                            >  2:4385481/1‑150 (MQ=255)
                                           tttAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACAt                                                                                                            <  1:1644913/150‑1 (MQ=255)
                                                aaGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGAc                                                                                                       >  2:2379495/1‑150 (MQ=255)
                                                aaGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGAc                                                                                                       >  1:4110156/1‑150 (MQ=255)
                                                 aGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACg                                                                                                      >  1:191125/1‑150 (MQ=255)
                                                  gggAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGt                                                                                                     >  2:2600410/1‑150 (MQ=255)
                                                        tttAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTttatt                                                                                               <  1:3916819/150‑1 (MQ=255)
                                                          tAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAg                                                                                             >  1:3586664/1‑150 (MQ=255)
                                                          tAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAg                                                                                             >  2:1725910/1‑150 (MQ=255)
                                                          tAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAg                                                                                             >  2:57763/1‑150 (MQ=255)
                                                          tAAAGAGGTAAGGCCGAGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCTGAAATAACATTGGACGTTTTATTAg                                                                                             >  2:3912409/1‑150 (MQ=255)
                                                             agagGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGAtt                                                                                          >  2:2249156/1‑150 (MQ=255)
                                                              gagGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGAttt                                                                                         <  2:1839120/150‑1 (MQ=255)
                                                              gagGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGAttt                                                                                         <  1:3078844/150‑1 (MQ=255)
                                                                ggTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGt                                                                                       <  2:1901387/150‑1 (MQ=255)
                                                                   aaGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTAtt                                                                                    <  2:4155424/150‑1 (MQ=255)
                                                                    aGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTc                                                                                   >  2:1464823/1‑150 (MQ=255)
                                                                      gCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCAc                                                                                 >  2:2305598/1‑150 (MQ=255)
                                                                      gCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCAc                                                                                 >  2:4540468/1‑150 (MQ=255)
                                                                          cgTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGa                                                                             >  1:3643277/1‑150 (MQ=255)
                                                                          cgTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGa                                                                             >  1:851367/1‑150 (MQ=255)
                                                                            ttATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGt                                                                           >  2:2623989/1‑150 (MQ=255)
                                                                            ttATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGt                                                                           >  1:1754525/1‑150 (MQ=255)
                                                                             tATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTg                                                                          >  2:181829/1‑150 (MQ=255)
                                                                               ttCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGAt                                                                        <  2:1224756/150‑1 (MQ=255)
                                                                                 cTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGt                                                                      >  2:838577/1‑150 (MQ=255)
                                                                                 cTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGt                                                                      >  1:1777936/1‑150 (MQ=255)
                                                                                   gTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTaa                                                                    >  1:1051708/1‑150 (MQ=255)
                                                                                    tATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAAt                                                                   <  2:4029285/150‑1 (MQ=255)
                                                                                      ttGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGa                                                                 >  1:288001/1‑150 (MQ=255)
                                                                                         ccTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAat                                                              <  1:98356/150‑1 (MQ=255)
                                                                                           tGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAatat                                                            <  2:2730092/150‑1 (MQ=255)
                                                                                           tGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAatat                                                            <  1:4521329/150‑1 (MQ=255)
                                                                                             atatatTCACACCTATAAATTTAGGGATATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAatatat                                                          <  2:2814345/150‑1 (MQ=255)
                                                                                                 atTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATattatt                                                      <  2:2043456/150‑1 (MQ=255)
                                                                                                 atTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATattatt                                                      >  2:3646781/1‑150 (MQ=255)
                                                                                                 atTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATattatt                                                      <  1:3011768/150‑1 (MQ=255)
                                                                                                  ttCACACCTCTAAATTTAGGGATATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttt                                                     <  1:4068206/150‑1 (MQ=255)
                                                                                                   tCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAtttt                                                    <  2:3465280/150‑1 (MQ=255)
                                                                                                     acacCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAtttttt                                                  <  2:3807319/150‑1 (MQ=255)
                                                                                                     acacCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAtttttt                                                  >  2:3562268/1‑150 (MQ=255)
                                                                                                      cacCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttttttt                                                 <  1:401569/150‑1 (MQ=255)
                                                                                                       acCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTa                                                >  1:966057/1‑150 (MQ=255)
                                                                                                       acCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTa                                                <  1:1597881/150‑1 (MQ=255)
                                                                                                         cTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTAtt                                              >  1:2139460/1‑150 (MQ=255)
                                                                                                          tataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTc                                             >  2:1707604/1‑150 (MQ=255)
                                                                                                          tataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTc                                             <  1:553919/150‑1 (MQ=255)
                                                                                                            taAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCag                                           <  1:1502194/150‑1 (MQ=255)
                                                                                                            taAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCag                                           >  1:1495332/1‑150 (MQ=255)
                                                                                                            taAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCag                                           <  1:1208243/150‑1 (MQ=255)
                                                                                                             aaaTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCaga                                          <  2:2465715/150‑1 (MQ=255)
                                                                                                              aaTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCagag                                         >  2:1179778/1‑150 (MQ=255)
                                                                                                              aaTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCagag                                         <  1:689262/150‑1 (MQ=255)
                                                                                                              aaTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCagag                                         <  1:2977373/150‑1 (MQ=255)
                                                                                                               aTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGc                                        >  1:943812/1‑150 (MQ=255)
                                                                                                               aTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGc                                        >  1:4074135/1‑150 (MQ=255)
                                                                                                                tttAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCt                                       <  1:2944882/150‑1 (MQ=255)
                                                                                                                tttAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCt                                       <  1:751290/150‑1 (MQ=255)
                                                                                                                 ttAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCta                                      >  2:907129/1‑150 (MQ=255)
                                                                                                                 ttAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCta                                      >  2:1080332/1‑150 (MQ=255)
                                                                                                                  tAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCtat                                     >  1:2449117/1‑150 (MQ=255)
                                                                                                                    gggCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTatat                                   <  1:3219332/150‑1 (MQ=255)
                                                                                                                      gCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtgg                                 >  1:3665408/1‑150 (MQ=255)
                                                                                                                       cTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtggc                                <  2:2925602/150‑1 (MQ=255)
                                                                                                                         aTTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtggctg                              >  2:741149/1‑150 (MQ=255)
                                                                                                                              aaaaCTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCttt                         >  1:4470023/1‑150 (MQ=255)
                                                                                                                                aaCTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTAt                       <  2:3511125/150‑1 (MQ=255)
                                                                                                                                      cAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCttt                 >  2:655744/1‑150 (MQ=255)
                                                                                                                                        aGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTa               >  2:1665428/1‑150 (MQ=255)
                                                                                                                                              gTCTAATTATTCTTTAGAAACATTTAAGATGAAGTGCTCCTGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTa         <  2:4005287/150‑1 (MQ=255)
                                                                                                                                                cTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTAtt       >  1:2251948/1‑150 (MQ=255)
                                                                                                                                                 tAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTc      >  1:2601700/1‑150 (MQ=255)
                                                                                                                                                  aaTGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCt     <  2:4145460/150‑1 (MQ=255)
                                                                                                                                                   aTGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCtt    >  1:236427/1‑150 (MQ=255)
                                                                                                                                                     gATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCtttt  <  2:3694975/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTT  >  CP009273/729000‑729298

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: