Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 997,965 G→T 94.5% D235Y (GAT→TAT)  elfG → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273997,9650GT94.5% 783.8 / 35.8 256D235Y (GAT→TAT) elfGputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (9/5);  new base T (125/117);  total (134/122)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 257 total aligned reads displayed.

CGCGATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCAG  >  CP009273/997816‑998110
                                                                                                                                                     |                                                                                                                                                 
    aTTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGt                                                                                                                                               >  1:2265730/1‑150 (MQ=255)
         atGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCg                                                                                                                                          <  1:2979796/150‑1 (MQ=255)
            aaTGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTg                                                                                                                                       >  2:1100296/1‑150 (MQ=255)
               gCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCa                                                                                                                                    <  2:3338843/150‑1 (MQ=255)
                 tACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAt                                                                                                                                  >  2:2791661/1‑150 (MQ=255)
                 tACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAAt                                                                                                                                  >  1:1623765/1‑150 (MQ=255)
                     acaTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTcc                                                                                                                              >  2:601345/1‑150 (MQ=255)
                       aTACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgt                                                                                                                            >  1:1732223/1‑150 (MQ=255)
                        tACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgtg                                                                                                                           <  1:1602348/150‑1 (MQ=255)
                        tACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgtg                                                                                                                           >  1:1652488/1‑150 (MQ=255)
                        tACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgtg                                                                                                                           >  1:2160712/1‑150 (MQ=255)
                         aCCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCgtgt                                                                                                                          <  2:1336896/150‑1 (MQ=255)
                           ccGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGtgttg                                                                                                                        >  1:347170/1‑150 (MQ=255)
                               aGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAgg                                                                                                                    >  1:2922902/1‑150 (MQ=255)
                               aGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAgg                                                                                                                    >  1:2398590/1‑150 (MQ=255)
                                ggTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAggg                                                                                                                   <  1:869520/150‑1 (MQ=255)
                                 gTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGc                                                                                                                  >  2:62485/1‑150 (MQ=255)
                                  tGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCa                                                                                                                 <  1:3449705/150‑1 (MQ=255)
                                   ggatCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCAt                                                                                                                >  1:4261850/4‑150 (MQ=255)
                                   gATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCAt                                                                                                                <  1:2027755/150‑1 (MQ=255)
                                    ataccTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATa                                                                                                               >  1:1961209/4‑150 (MQ=255)
                                    aTTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATa                                                                                                               <  1:2337659/150‑1 (MQ=255)
                                       ccTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAAc                                                                                                            <  2:2133266/150‑1 (MQ=255)
                                       ccTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAAc                                                                                                            <  2:2068667/150‑1 (MQ=255)
                                        cTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAAcc                                                                                                           >  1:1102671/1‑150 (MQ=255)
                                         tGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACct                                                                                                          <  2:819494/150‑1 (MQ=255)
                                              gACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTgg                                                                                                     <  2:200464/150‑1 (MQ=255)
                                               aCTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTggg                                                                                                    >  1:2685636/1‑150 (MQ=255)
                                               aCTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTggg                                                                                                    >  1:3534536/1‑150 (MQ=255)
                                                 tACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTAAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGc                                                                                                  <  1:3361430/150‑1 (MQ=255)
                                                   cAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCa                                                                                                <  1:3830222/150‑1 (MQ=255)
                                                      gTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGaa                                                                                             >  2:2959476/1‑150 (MQ=255)
                                                       tACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAAc                                                                                            <  1:2782249/150‑1 (MQ=255)
                                                              gTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGaa                                                                                     <  1:2384361/150‑1 (MQ=255)
                                                                 ttAGTTTATTACCTGTCTGGTACGATCAAGTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACtt                                                                                  <  2:3072300/150‑1 (MQ=255)
                                                                 ttAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACtt                                                                                  <  2:809632/150‑1 (MQ=255)
                                                                  tAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTg                                                                                 <  1:1094187/150‑1 (MQ=255)
                                                                  tAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTg                                                                                 >  1:3754197/1‑150 (MQ=255)
                                                                      gtATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAAt                                                                             <  1:3880874/150‑1 (MQ=255)
                                                                          tACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCa                                                                         >  2:2546449/1‑150 (MQ=255)
                                                                          tACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCa                                                                         <  1:1037256/150‑1 (MQ=255)
                                                                           aCCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAg                                                                        >  2:1338885/1‑150 (MQ=255)
                                                                            ccTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGt                                                                       >  1:54933/1‑150 (MQ=255)
                                                                               gTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAg                                                                    <  2:4191159/150‑1 (MQ=255)
                                                                                tCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAgg                                                                   >  2:4436436/1‑150 (MQ=255)
                                                                                  tGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgtg                                                                 >  2:923103/1‑150 (MQ=255)
                                                                                   ggTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGgtgt                                                                >  2:1372395/1‑150 (MQ=255)
                                                                                        gATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAacac                                                           >  1:4380925/1‑150 (MQ=255)
                                                                                          tCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACgg                                                         >  2:4085299/1‑150 (MQ=255)
                                                                                           cAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACggg                                                        >  2:2940449/1‑150 (MQ=255)
                                                                                           cAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACggg                                                        <  1:2579549/150‑1 (MQ=255)
                                                                                            aaTTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGa                                                       <  2:1978862/150‑1 (MQ=255)
                                                                                            aaTTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGAGTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGa                                                       >  1:3974487/1‑150 (MQ=255)
                                                                                                cACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAAcg                                                   <  2:2330660/150‑1 (MQ=255)
                                                                                                     ggCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGa                                              <  1:972815/150‑1 (MQ=255)
                                                                                                       cTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGAtt                                            <  1:3732956/150‑1 (MQ=255)
                                                                                                        tCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTa                                           >  1:1213907/1‑150 (MQ=255)
                                                                                                        tCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTa                                           <  1:4226617/150‑1 (MQ=255)
                                                                                                        tCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTa                                           <  1:1894949/150‑1 (MQ=255)
                                                                                                          aTGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                         >  2:4411803/1‑150 (MQ=255)
                                                                                                          aTGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGt                                         >  2:3223590/1‑150 (MQ=255)
                                                                                                             ttCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTaa                                      <  2:4248610/150‑1 (MQ=255)
                                                                                                               ccGTCAATGCCTGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCAGTGTTGTAGGGCATAAACCTCTTGGGTCCAGAACGGCAGAACTTGCAATTCCATTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATg                                    <  1:658199/150‑1 (MQ=255)
                                                                                                               ccGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATg                                    <  1:2145036/150‑1 (MQ=255)
                                                                                                               ccGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATg                                    <  1:3572579/150‑1 (MQ=255)
                                                                                                                  tCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCa                                 <  1:3372068/150‑1 (MQ=255)
                                                                                                                  tCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCa                                 <  1:593977/150‑1 (MQ=255)
                                                                                                                   cAATGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCaa                                >  2:2708592/1‑150 (MQ=255)
                                                                                                                    aaTGCCGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAAc                               >  2:3580796/1‑150 (MQ=255)
                                                                                                                        ccGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGa                           <  2:2398547/150‑1 (MQ=255)
                                                                                                                         cGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                          >  2:3089430/1‑150 (MQ=255)
                                                                                                                         cGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                          >  2:3367492/1‑150 (MQ=255)
                                                                                                                         cGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                          >  2:3707683/1‑150 (MQ=255)
                                                                                                                         cGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                          >  1:4632752/1‑150 (MQ=255)
                                                                                                                         cGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                          >  2:474059/1‑150 (MQ=255)
                                                                                                                         cGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                          >  2:1595450/1‑150 (MQ=255)
                                                                                                                         cGGAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAg                          >  2:1045494/1‑150 (MQ=255)
                                                                                                                          ggAGAGACAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGt                         <  1:1420129/150‑1 (MQ=255)
                                                                                                                          ggAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGt                         >  2:179751/1‑150 (MQ=255)
                                                                                                                           gagagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTc                        >  1:731840/1‑150 (MQ=255)
                                                                                                                           gagagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTc                        >  1:3096149/1‑150 (MQ=255)
                                                                                                                           gagagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTc                        >  2:3889371/1‑150 (MQ=255)
                                                                                                                            agagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCt                       >  2:2584854/1‑150 (MQ=255)
                                                                                                                            agagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCt                       <  1:2215968/150‑1 (MQ=255)
                                                                                                                             gagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTg                      <  1:3589097/150‑1 (MQ=255)
                                                                                                                             gagaCAATTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTg                      >  2:358442/1‑150 (MQ=255)
                                                                                                                                aCAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACg                   >  1:4485727/1‑150 (MQ=255)
                                                                                                                                  aaTTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGc                 >  1:4374030/1‑150 (MQ=255)
                                                                                                                                  aaTTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGc                 >  1:3280064/1‑150 (MQ=255)
                                                                                                                                   aTTGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCa                >  1:946843/1‑150 (MQ=255)
                                                                                                                                    ttGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCaa               >  2:4630410/1‑150 (MQ=255)
                                                                                                                                     tGAAGTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAAc              >  2:114446/1‑150 (MQ=255)
                                                                                                                                         gTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACcaca          >  1:4166156/1‑150 (MQ=255)
                                                                                                                                         gTTGATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACcaca          <  1:427396/150‑1 (MQ=255)
                                                                                                                                            gATTTAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAc       >  1:1299392/1‑150 (MQ=255)
                                                                                                                                               ttAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAcccc    <  2:4230091/150‑1 (MQ=255)
                                                                                                                                               ttAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGAcccc    <  2:1623765/150‑1 (MQ=255)
                                                                                                                                                tAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCa   <  2:1885981/150‑1 (MQ=255)
                                                                                                                                                tAGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCa   >  1:4141775/1‑150 (MQ=255)
                                                                                                                                                 aGGTTATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCAg  >  2:1447351/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGCGATTATATGAATGCTACAACATACCCGCAGGTGATTCCTGCACGACTACAGGTACACCGGTTTTAGTGTATTACCTGTCTGGTACGATCAATTCACTTGGCTCATGTTCCGTCAATGCCGGAGAGACAATTGAAGTTGATTTAGGTGATGTCTTCGCTGCCAATTTCCGTGTTGTAGGGCATAAACCTCTTGGGGCCAGAACGGCAGAACTTGCAATTCCAGTCAGGTGTAACACGGGAAACGCGGGATTAGTTAATGTCAACCTGAGTCTGACGGCAACCACAGACCCCAG  >  CP009273/997816‑998110

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: