Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 729,149 G→A 15.0% 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,1490GA15.0% 833.3 / 124.1 288D1371N (GAT→AAT) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (88/153);  new base A (10/33);  total (98/190)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 288 total aligned reads displayed.

TTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTT  >  CP009273/729001‑729297
                                                                                                                                                    |                                                                                                                                                    
 tgtgACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGAt                                                                                                                                                    >  2:2667706/1‑150 (MQ=255)
 tgtgACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGAt                                                                                                                                                    <  2:3953634/150‑1 (MQ=255)
  gtgACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGAtt                                                                                                                                                   <  2:1710201/150‑1 (MQ=255)
  gtgACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAAtt                                                                                                                                                   >  2:4740204/1‑150 (MQ=255)
   tgACCTCGTTGAAACCGCTAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTTGTCTAATGATTc                                                                                                                                                  <  2:3333890/150‑1 (MQ=255)
    gACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCt                                                                                                                                                 <  1:2865781/150‑1 (MQ=255)
     aCCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCtt                                                                                                                                                <  1:2448193/150‑1 (MQ=255)
      ccTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCttt                                                                                                                                               <  1:697899/150‑1 (MQ=255)
      ccTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCttt                                                                                                                                               <  2:1651390/150‑1 (MQ=255)
       cTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTa                                                                                                                                              <  2:2879885/150‑1 (MQ=255)
           tnGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATCAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGaaa                                                                                                                                          >  1:1413924/1‑150 (MQ=255)
             gAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTCGAAACa                                                                                                                                        >  1:3519537/1‑150 (MQ=255)
                aCCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGt                                                                                                                                     <  1:2920650/150‑1 (MQ=255)
                aCCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGt                                                                                                                                     >  2:3747416/1‑150 (MQ=255)
                 ccGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTa                                                                                                                                    <  2:4341705/150‑1 (MQ=255)
                  cGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTATGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTaa                                                                                                                                   <  1:3630979/150‑1 (MQ=255)
                   gATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAg                                                                                                                                  <  2:4263665/150‑1 (MQ=255)
                    ataataGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGa                                                                                                                                 <  1:529615/150‑1 (MQ=255)
                    ataataGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGa                                                                                                                                 <  1:1756622/150‑1 (MQ=255)
                      aataGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATg                                                                                                                               <  1:3500516/150‑1 (MQ=255)
                       ataGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATTATTCTTTAGAAACATGTAAGATGa                                                                                                                              <  1:4671709/150‑1 (MQ=255)
                         aGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAg                                                                                                                            >  2:613010/1‑150 (MQ=255)
                             gtgtgGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCt                                                                                                                        <  1:4339886/150‑1 (MQ=255)
                              tgtgGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTc                                                                                                                       <  2:3040241/150‑1 (MQ=255)
                                tgGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCg                                                                                                                     >  2:2342817/1‑150 (MQ=255)
                                 ggAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCgg                                                                                                                    >  2:3360424/1‑150 (MQ=255)
                                  gagaACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGa                                                                                                                   <  2:3533417/150‑1 (MQ=255)
                                    gaACTGTTTAAAAGGGAAGTTTATAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGaaa                                                                                                                 <  1:448559/150‑1 (MQ=255)
                                    gaACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGaaa                                                                                                                 >  2:2360773/1‑150 (MQ=255)
                                     aaCTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAAt                                                                                                                >  1:4186051/1‑150 (MQ=255)
                                       cTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATaa                                                                                                              <  2:1163923/150‑1 (MQ=255)
                                        tGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAAc                                                                                                             >  1:2616418/1‑150 (MQ=255)
                                           ttAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACAtt                                                                                                          <  1:1494817/150‑1 (MQ=255)
                                             aaaaGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTgg                                                                                                        <  2:2445498/150‑1 (MQ=255)
                                             aaaaGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTgg                                                                                                        >  2:3639950/1‑150 (MQ=255)
                                              aaaGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGa                                                                                                       <  1:1366344/150‑1 (MQ=255)
                                                  ggAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGtt                                                                                                   >  2:2778194/1‑150 (MQ=255)
                                                  ggAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGtt                                                                                                   <  2:5055375/150‑1 (MQ=255)
                                                   gAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGttt                                                                                                  <  2:686899/150‑1 (MQ=255)
                                                    aaGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGtttt                                                                                                 <  2:3811323/150‑1 (MQ=255)
                                                    aaGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGtttt                                                                                                 >  1:470571/1‑150 (MQ=255)
                                                     aGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTtta                                                                                                <  1:4816549/150‑1 (MQ=255)
                                                             gagGTAAGGCCGCGTTATTCTGTATTGCCTGATAGAGTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGACATAACATTGGACGTTTTATTAGAttt                                                                                        <  2:3929852/150‑1 (MQ=255)
                                                               ggTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGt                                                                                      <  1:4109214/150‑1 (MQ=255)
                                                                gTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTa                                                                                     <  2:3216317/150‑1 (MQ=255)
                                                                  aaGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTAtt                                                                                   >  1:892576/1‑150 (MQ=255)
                                                                     gCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCAc                                                                                <  1:2148348/150‑1 (MQ=255)
                                                                      ccGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCAcc                                                                               <  1:638481/150‑1 (MQ=255)
                                                                      ccGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCAcc                                                                               >  1:2592479/1‑150 (MQ=255)
                                                                       cgcgTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCg                                                                              <  2:2231561/150‑1 (MQ=255)
                                                                         cgTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGa                                                                            <  2:4944759/150‑1 (MQ=255)
                                                                          gTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAg                                                                           <  1:3815041/150‑1 (MQ=255)
                                                                          gTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAg                                                                           <  2:1300084/150‑1 (MQ=255)
                                                                            tATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTg                                                                         >  2:4426668/1‑150 (MQ=255)
                                                                             aTTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGa                                                                        >  1:2077395/1‑150 (MQ=255)
                                                                             aTTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGa                                                                        >  1:196191/1‑150 (MQ=255)
                                                                                cTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGt                                                                     <  2:4186051/150‑1 (MQ=255)
                                                                                 tGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTa                                                                    >  1:4833197/1‑150 (MQ=255)
                                                                                    aTTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATg                                                                 >  2:2090462/1‑150 (MQ=255)
                                                                                      tGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGaa                                                               <  1:5097051/150‑1 (MQ=255)
                                                                                         cTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAata                                                            >  2:3371995/1‑150 (MQ=255)
                                                                                         cTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACAAGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGAGAGGTAATGAAata                                                            >  1:2742965/1‑150 (MQ=255)
                                                                                           gATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAatata                                                          >  1:1012683/1‑150 (MQ=255)
                                                                                            atatatTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAatatat                                                         <  2:4609543/150‑1 (MQ=255)
                                                                                              atatTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATAtatta                                                       >  1:3089925/1‑150 (MQ=255)
                                                                                                 ttCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttt                                                    <  1:4029866/150‑1 (MQ=255)
                                                                                                 ttCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttt                                                    <  2:983468/150‑1 (MQ=255)
                                                                                                   cacaCCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttttt                                                  <  1:1976497/150‑1 (MQ=255)
                                                                                                   cacaCCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttttt                                                  <  2:1798214/150‑1 (MQ=255)
                                                                                                     cacCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttttttt                                                <  2:175163/150‑1 (MQ=255)
                                                                                                     cacCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAttttttt                                                <  1:1706826/150‑1 (MQ=255)
                                                                                                       catataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTAt                                              >  2:4707024/3‑150 (MQ=255)
                                                                                                        cTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTAtt                                             <  2:964720/150‑1 (MQ=255)
                                                                                                          ataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCa                                           <  1:4208037/150‑1 (MQ=255)
                                                                                                          ataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCa                                           <  1:3569139/150‑1 (MQ=255)
                                                                                                           taAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCag                                          <  1:1213629/150‑1 (MQ=255)
                                                                                                           taAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCag                                          <  2:1515508/150‑1 (MQ=255)
                                                                                                             aaTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCagag                                        >  1:3541634/1‑150 (MQ=255)
                                                                                                              aTTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGc                                       <  1:3220609/150‑1 (MQ=255)
                                                                                                               tttAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCt                                      <  1:4961780/150‑1 (MQ=255)
                                                                                                                ttAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCta                                     >  1:1314539/1‑150 (MQ=255)
                                                                                                                 tAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCtat                                    >  1:5062521/1‑150 (MQ=255)
                                                                                                                       tATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtggct                              <  2:1975321/150‑1 (MQ=255)
                                                                                                                        aTTTAAAAACTGCCAAGACTGGTCTAATGATTATTTAGAAACAGGTAAGATGAAGTGCACCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGACATATATTATTTTTTTATTCAGAGCTATAtggctg                             >  2:4839976/1‑150 (MQ=255)
                                                                                                                         tttAAAAACTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtggctgg                            <  2:2561707/150‑1 (MQ=255)
                                                                                                                            aaaaaCTGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCtt                         <  2:2521829/150‑1 (MQ=255)
                                                                                                                              aaaCTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGATGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTa                       >  1:234793/1‑150 (MQ=255)
                                                                                                                               aaCTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTAt                      >  2:884607/1‑150 (MQ=255)
                                                                                                                                  tGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGc                   >  2:1307626/1‑150 (MQ=255)
                                                                                                                                  tGCCAAGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGc                   <  2:2838577/150‑1 (MQ=255)
                                                                                                                                   gCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCt                  <  1:2121456/150‑1 (MQ=255)
                                                                                                                                     cAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCttt                >  1:4415605/1‑150 (MQ=255)
                                                                                                                                      aaGACTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCtttt               <  2:816349/150‑1 (MQ=255)
                                                                                                                                       aGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTa              >  2:4514485/1‑150 (MQ=255)
                                                                                                                                          cTGGTCTAATAATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATa           <  1:1841945/150‑1 (MQ=255)
                                                                                                                                            ggTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACt         >  1:3704404/1‑150 (MQ=255)
                                                                                                                                             gTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTa        <  1:3150508/150‑1 (MQ=255)
                                                                                                                                              tCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTAt       >  2:1155435/1‑150 (MQ=255)
                                                                                                                                               cTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTAtt      >  2:1400613/1‑150 (MQ=255)
                                                                                                                                                tAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTc     <  1:3404904/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGTGACCTCGTTGAAACCGATAATAGTTGTGTGGAGAACTGTTTAAAAGGGAAGTTTAAAGAGGTAAGGCCGCGTTATTCTGTATTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTT  >  CP009273/729001‑729297

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: