Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,505,141 G→T 100% R452S (CGT→AGT)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,505,1410GT98.5% 857.0 / 2.6 264R452S (CGT→AGT) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (100/159);  total (102/162)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 264 total aligned reads displayed.

ACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGC  >  CP009273/4504995‑4505286
                                                                                                                                                  |                                                                                                                                                 
aCGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCa                                                                                                                                                 >  2:3068026/1‑149 (MQ=255)
aCGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCa                                                                                                                                                 >  1:3187443/1‑149 (MQ=255)
 cGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCAt                                                                                                                                                <  2:3516686/149‑1 (MQ=255)
      aCCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCCCGGTCGTTAGGGCTTGAGCCGGACGGCAACAGCCCGCTGCTGGTTGCGGTGTAGTAACTCATTTCAt                                                                                                                                           >  1:3793570/1‑149 (MQ=255)
       cccGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATg                                                                                                                                          <  2:2843560/149‑1 (MQ=255)
          ggCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGt                                                                                                                                       >  2:1656114/1‑149 (MQ=255)
          ggCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGt                                                                                                                                       <  2:781064/149‑1 (MQ=255)
           gCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGtt                                                                                                                                      >  2:825303/1‑149 (MQ=255)
           gCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATGTCATGCGtt                                                                                                                                      >  2:825289/1‑149 (MQ=255)
             gTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGa                                                                                                                                    <  2:2212442/149‑1 (MQ=255)
             gTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGa                                                                                                                                    <  1:771288/149‑1 (MQ=255)
                  ggTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcat                                                                                                                               <  1:1269447/149‑1 (MQ=255)
                   gTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcatt                                                                                                                              >  1:249462/1‑149 (MQ=255)
                    tCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcattc                                                                                                                             >  2:3590300/1‑149 (MQ=255)
                      cAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCaa                                                                                                                           >  2:2442415/1‑149 (MQ=255)
                       aGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCaaa                                                                                                                          >  2:1719274/1‑149 (MQ=255)
                        gTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAAt                                                                                                                         <  2:573011/149‑1 (MQ=255)
                         ttGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATa                                                                                                                        <  1:3789459/149‑1 (MQ=255)
                         ttGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATa                                                                                                                        <  1:938007/149‑1 (MQ=255)
                         ttGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATa                                                                                                                        >  2:3147117/1‑149 (MQ=255)
                             cGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGa                                                                                                                    <  1:1062113/149‑1 (MQ=255)
                             cGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGa                                                                                                                    <  1:3031986/149‑1 (MQ=255)
                             cGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGa                                                                                                                    <  1:1595594/149‑1 (MQ=255)
                              gATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGAt                                                                                                                   <  2:877721/149‑1 (MQ=255)
                               aTGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATa                                                                                                                  <  1:2523971/149‑1 (MQ=255)
                                 gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgc                                                                                                                >  2:2083441/1‑149 (MQ=255)
                                 gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgc                                                                                                                <  2:3372401/149‑1 (MQ=255)
                                  tCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcc                                                                                                               <  2:1031012/149‑1 (MQ=255)
                                   cGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgccc                                                                                                              <  1:591521/149‑1 (MQ=255)
                                     aTTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccac                                                                                                            >  2:804318/1‑149 (MQ=255)
                                      ttttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacg                                                                                                           >  2:142767/1‑149 (MQ=255)
                                      ttttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacg                                                                                                           >  2:1755279/1‑149 (MQ=255)
                                       tttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgc                                                                                                          >  1:2929844/1‑149 (MQ=255)
                                        ttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgcc                                                                                                         <  1:19966/149‑1 (MQ=255)
                                        ttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgcc                                                                                                         <  1:2765166/149‑1 (MQ=255)
                                         tGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgccc                                                                                                        <  2:1567085/149‑1 (MQ=255)
                                             aTCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACtt                                                                                                    >  1:844763/1‑149 (MQ=255)
                                                cAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGATCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGt                                                                                                 >  1:3605082/1‑149 (MQ=255)
                                                  gATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTgc                                                                                               <  1:3966282/149‑1 (MQ=255)
                                                  gATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTgc                                                                                               >  1:1943388/1‑149 (MQ=255)
                                                     aCCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGcgcg                                                                                            >  2:2482835/1‑149 (MQ=255)
                                                         ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAg                                                                                        >  1:4098109/1‑149 (MQ=255)
                                                          gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAgg                                                                                       <  2:3064068/149‑1 (MQ=255)
                                                            gTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGc                                                                                     >  1:3866407/1‑149 (MQ=255)
                                                            gTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGc                                                                                     >  1:3762563/1‑149 (MQ=255)
                                                             tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGcc                                                                                    >  2:516444/1‑149 (MQ=255)
                                                               cgcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGa                                                                                  <  2:1961624/149‑1 (MQ=255)
                                                                gcCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGAt                                                                                 >  2:3223704/1‑149 (MQ=255)
                                                                   tCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCAt                                                                              <  1:1826663/149‑1 (MQ=255)
                                                                     ggTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaa                                                                            >  2:2404358/1‑149 (MQ=255)
                                                                     ggTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaa                                                                            <  2:2862496/149‑1 (MQ=255)
                                                                      gTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaaa                                                                           <  2:3713944/149‑1 (MQ=255)
                                                                        gCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGa                                                                         >  1:3129216/1‑149 (MQ=255)
                                                                           ggAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGAACATAAAGATCt                                                                      >  1:572970/1‑149 (MQ=255)
                                                                            gAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTg                                                                     <  2:1206788/149‑1 (MQ=255)
                                                                            gAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGATGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTg                                                                     <  1:1112059/149‑1 (MQ=255)
                                                                              aCGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGc                                                                   <  1:2066060/149‑1 (MQ=255)
                                                                              aCGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGc                                                                   <  1:976633/149‑1 (MQ=255)
                                                                                 cgTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGt                                                                <  2:404307/149‑1 (MQ=255)
                                                                                    aTCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGc                                                             <  2:3187443/149‑1 (MQ=255)
                                                                                     tCGCGGTCGTAAGGGCTTGATCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg                                                            <  2:3605082/149‑1 (MQ=255)
                                                                                       gcgGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTg                                                          <  2:481255/149‑1 (MQ=255)
                                                                                         ggTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGc                                                        <  1:2810017/149‑1 (MQ=255)
                                                                                          gTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCt                                                       >  2:1608771/1‑149 (MQ=255)
                                                                                              tAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAAt                                                   <  2:124988/149‑1 (MQ=255)
                                                                                              tAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAAt                                                   <  2:794862/149‑1 (MQ=255)
                                                                                               aaGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATa                                                  <  2:414831/149‑1 (MQ=255)
                                                                                               aaGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATa                                                  >  2:1274959/1‑149 (MQ=255)
                                                                                                aGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                 <  2:3818910/149‑1 (MQ=255)
                                                                                                aGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                 <  1:3952083/149‑1 (MQ=255)
                                                                                                aGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                 <  1:804318/149‑1 (MQ=255)
                                                                                                aGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                 <  1:1737439/149‑1 (MQ=255)
                                                                                                          ccGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCa                                       <  2:1378835/149‑1 (MQ=255)
                                                                                                            ggACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagt                                     <  2:1601899/149‑1 (MQ=255)
                                                                                                             gACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagta                                    >  1:2404129/1‑149 (MQ=255)
                                                                                                              aCGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagtag                                   >  2:2510131/1‑149 (MQ=255)
                                                                                                                 gCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTa                                >  2:161463/1‑149 (MQ=255)
                                                                                                                    aCTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAcgc                             <  2:398371/149‑1 (MQ=255)
                                                                                                                            ctgctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATATCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGagag                     <  1:2660802/149‑1 (MQ=255)
                                                                                                                            ctgctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGagag                     <  2:2718814/149‑1 (MQ=255)
                                                                                                                             tgctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAgg                    <  1:3623540/149‑1 (MQ=255)
                                                                                                                             tgctgGTGGCGGTGTAGTAACTCATGTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAgg                    <  2:2631806/149‑1 (MQ=255)
                                                                                                                              gctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAggg                   <  1:1640217/149‑1 (MQ=255)
                                                                                                                              gctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAggg                   <  1:3590300/149‑1 (MQ=255)
                                                                                                                               ctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGt                  >  2:286220/1‑149 (MQ=255)
                                                                                                                                 ggtggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGa                <  2:4003591/149‑1 (MQ=255)
                                                                                                                                 ggtggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGa                <  2:2656377/149‑1 (MQ=255)
                                                                                                                                 ggtggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGa                <  1:2250174/149‑1 (MQ=255)
                                                                                                                                  gtggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCTGCGAGAGGGTGAt               <  1:532746/149‑1 (MQ=255)
                                                                                                                                   tggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAAGAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATg              <  2:3960428/149‑1 (MQ=255)
                                                                                                                                       ggTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGtt          >  1:330860/1‑149 (MQ=255)
                                                                                                                                       ggTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGtt          <  1:3130272/149‑1 (MQ=255)
                                                                                                                                       ggTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGtt          <  1:2965992/149‑1 (MQ=255)
                                                                                                                                          gtagtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGc       <  1:1434194/149‑1 (MQ=255)
                                                                                                                                          gtagtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGc       >  1:2149885/1‑149 (MQ=255)
                                                                                                                                          gtagtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGc       >  1:2878116/1‑149 (MQ=255)
                                                                                                                                          gtagtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGc       <  1:417585/149‑1 (MQ=255)
                                                                                                                                            agtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATGCCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCt     <  2:3704324/149‑1 (MQ=255)
                                                                                                                                              taACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTg   <  1:1719274/149‑1 (MQ=255)
                                                                                                                                               aaCTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGc  >  1:2601990/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGC  >  CP009273/4504995‑4505286

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: