Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 615,143 C→A 27.4% T171T (ACG→ACT fepC ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273615,1430CA27.4% 750.3 / 241.7 338T171T (ACG→ACTfepCiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (138/107);  new base A (48/45);  total (186/152)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

TTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTC  >  CP009273/614994‑615292
                                                                                                                                                     |                                                                                                                                                     
   ccttcccGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATTTCCAGCCAGGTCGTg                                                                                                                                                    <  2:4825552/150‑1 (MQ=255)
      tcccGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGt                                                                                                                                                 <  1:2180273/150‑1 (MQ=255)
        ccGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTtc                                                                                                                                               >  2:1059924/1‑150 (MQ=255)
         cGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTtcg                                                                                                                                              >  2:569721/1‑150 (MQ=255)
            aatgcaatCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTtcgtcg                                                                                                                                           >  1:380526/1‑150 (MQ=255)
            aatgcaatCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTtcgtcg                                                                                                                                           <  1:4624940/150‑1 (MQ=255)
             atgcaatCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGa                                                                                                                                          <  1:4661915/150‑1 (MQ=255)
              tgcaatCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGag                                                                                                                                         >  2:1524465/1‑150 (MQ=255)
                caatCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGagca                                                                                                                                       >  2:796561/1‑150 (MQ=255)
                    cAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCAt                                                                                                                                   >  1:4980516/1‑150 (MQ=255)
                    cAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCAt                                                                                                                                   >  1:4392627/1‑150 (MQ=255)
                    cAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCAt                                                                                                                                   >  1:1058436/1‑150 (MQ=255)
                    cAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCAt                                                                                                                                   >  1:4106730/1‑150 (MQ=255)
                     aaaTGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATa                                                                                                                                  >  1:852666/1‑150 (MQ=255)
                      aaTGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATaa                                                                                                                                 <  1:799715/150‑1 (MQ=255)
                        tggctggcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAAtt                                                                                                                               >  2:3624408/1‑150 (MQ=255)
                           ctggcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGcc                                                                                                                            <  2:2346631/150‑1 (MQ=255)
                           ctggcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGcc                                                                                                                            <  1:3080635/150‑1 (MQ=255)
                            tggcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCg                                                                                                                           >  1:2728046/1‑150 (MQ=255)
                             ggcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGt                                                                                                                          <  2:276433/150‑1 (MQ=255)
                             ggcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGt                                                                                                                          <  2:4393385/150‑1 (MQ=255)
                             ggcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGt                                                                                                                          >  1:3938498/1‑150 (MQ=255)
                              gcGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGtt                                                                                                                         <  2:5035242/150‑1 (MQ=255)
                               cGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGttt                                                                                                                        >  2:13436/1‑150 (MQ=255)
                               cGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGttt                                                                                                                        <  1:2252369/150‑1 (MQ=255)
                                   aCGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGGTCGCTGAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTg                                                                                                                    <  1:1491552/150‑1 (MQ=255)
                                    cGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTgg                                                                                                                   >  1:1986701/1‑150 (MQ=255)
                                     gACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTggg                                                                                                                  >  2:3427610/1‑150 (MQ=255)
                                     gACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTggg                                                                                                                  <  2:4722057/150‑1 (MQ=255)
                                            cTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCAc                                                                                                           >  2:3048161/1‑150 (MQ=255)
                                            cTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCAc                                                                                                           >  2:539260/1‑150 (MQ=255)
                                               aTTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCAt                                                                                                        >  2:4348078/1‑150 (MQ=255)
                                               aTTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCAt                                                                                                        <  1:4074351/150‑1 (MQ=255)
                                                 tAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATcg                                                                                                      >  2:1991784/1‑150 (MQ=255)
                                                    gATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCCCCATCGCGa                                                                                                   <  1:5146114/150‑1 (MQ=255)
                                                    gATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGa                                                                                                   <  1:3024011/150‑1 (MQ=255)
                                                       cGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATcc                                                                                                <  2:2115607/150‑1 (MQ=255)
                                                         tGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAc                                                                                              <  1:3397655/150‑1 (MQ=255)
                                                          gcagcaCCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAcg                                                                                             >  2:38864/1‑150 (MQ=255)
                                                           cagcaCCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAcgc                                                                                            <  2:3278852/150‑1 (MQ=255)
                                                            agcaCCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAcgct                                                                                           >  1:1227349/1‑150 (MQ=255)
                                                            agcaCCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAcgct                                                                                           >  1:2958271/1‑150 (MQ=255)
                                                              caCCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAcgctcg                                                                                         >  2:258348/1‑150 (MQ=255)
                                                              caCCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAcgctcg                                                                                         >  1:2831958/1‑150 (MQ=255)
                                                                ccgccgCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCAcgctcgct                                                                                       <  2:3556815/150‑1 (MQ=255)
                                                                   ccgccAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGc                                                                                    <  1:4797578/150‑1 (MQ=255)
                                                                   ccgccAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGc                                                                                    <  2:162137/150‑1 (MQ=255)
                                                                    cgccAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCg                                                                                   >  1:3095831/1‑150 (MQ=255)
                                                                     gccAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGt                                                                                  >  2:5122045/1‑150 (MQ=255)
                                                                      ccAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGtt                                                                                 >  2:2321929/1‑150 (MQ=255)
                                                                         gggTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTc                                                                              >  2:4803619/1‑150 (MQ=255)
                                                                           gTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCg                                                                            <  2:1178784/150‑1 (MQ=255)
                                                                             ataGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGcc                                                                          <  2:901958/150‑1 (MQ=255)
                                                                              taGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCa                                                                         >  1:2208707/1‑150 (MQ=255)
                                                                               aGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAg                                                                        <  1:1424963/150‑1 (MQ=255)
                                                                                 ccTTTGTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGaa                                                                      <  2:3894480/150‑1 (MQ=255)
                                                                                  cTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGaaa                                                                     <  1:2099279/150‑1 (MQ=255)
                                                                                  cTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGaaa                                                                     <  2:1748183/150‑1 (MQ=255)
                                                                                    ttttCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAgg                                                                   <  1:4538275/150‑1 (MQ=255)
                                                                                    ttttCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAgg                                                                   >  1:4950782/1‑150 (MQ=255)
                                                                                        cgcgGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTAt                                                               >  2:2201247/1‑150 (MQ=255)
                                                                                          cgGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATcc                                                             <  1:1076353/150‑1 (MQ=255)
                                                                                            gTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCAc                                                           >  2:4870104/1‑150 (MQ=255)
                                                                                               cAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCt                                                        >  2:719572/1‑150 (MQ=255)
                                                                                                  ttCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTg                                                     <  2:1389726/150‑1 (MQ=255)
                                                                                                  ttCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTg                                                     <  2:413713/150‑1 (MQ=255)
                                                                                                  gtCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTg                                                     <  2:3707340/149‑1 (MQ=255)
                                                                                                   tCGCTTAACAACTCCAGCAAATCACTCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCATCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGa                                                    <  1:4099212/150‑1 (MQ=255)
                                                                                                    cGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGAt                                                   <  1:186072/150‑1 (MQ=255)
                                                                                                    cGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGAt                                                   >  2:1890963/1‑150 (MQ=255)
                                                                                                    cGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGAt                                                   <  2:4356753/150‑1 (MQ=255)
                                                                                                     gCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATc                                                  >  2:2455664/1‑150 (MQ=255)
                                                                                                      cTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACTCTTAGATCt                                                 >  1:1606155/1‑150 (MQ=255)
                                                                                                      cTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCt                                                 <  2:2418970/150‑1 (MQ=255)
                                                                                                       ttAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCta                                                <  2:4106730/150‑2 (MQ=255)
                                                                                                         aacaacTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCAGAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTCTCCACGCTTTGATCTGcc                                              >  2:1466887/1‑150 (MQ=255)
                                                                                                            aacTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGa                                           <  1:3488361/150‑1 (MQ=255)
                                                                                                                 cAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGt                                      <  2:3084071/150‑1 (MQ=255)
                                                                                                                  aGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGtt                                     >  2:699303/1‑150 (MQ=255)
                                                                                                                  aGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGtt                                     >  2:3129184/1‑150 (MQ=255)
                                                                                                                   gCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTa                                    <  2:4055424/150‑1 (MQ=255)
                                                                                                                     aaaTCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTAtt                                  >  2:649692/1‑150 (MQ=255)
                                                                                                                     aaaTCAATCTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTAtt                                  <  1:397268/150‑1 (MQ=255)
                                                                                                                       atcaatcTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTcc                                >  1:4840736/1‑150 (MQ=255)
                                                                                                                       atcaatcTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTcc                                >  1:1802737/1‑150 (MQ=255)
                                                                                                                         caatcTGATGACTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCg                              >  2:4758871/1‑150 (MQ=255)
                                                                                                                              tgatgaCTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGcc                         <  2:4220164/150‑1 (MQ=255)
                                                                                                                              tgatgaCTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGcc                         <  1:992054/150‑1 (MQ=255)
                                                                                                                               gatgaCTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCt                        >  1:2761561/1‑150 (MQ=255)
                                                                                                                                  gaCTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCa                     >  2:2561752/1‑150 (MQ=255)
                                                                                                                                    cTGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCAtt                   >  1:2328759/1‑150 (MQ=255)
                                                                                                                                     tGATATCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTg                  >  2:126927/1‑150 (MQ=255)
                                                                                                                                      gATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGc                 <  1:569721/150‑1 (MQ=255)
                                                                                                                                         atCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCttt              <  1:539317/150‑1 (MQ=255)
                                                                                                                                         atCCAGCCAGGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAGAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCttt              >  2:4575232/1‑150 (MQ=255)
                                                                                                                                          tCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCTTTc             <  2:4004410/150‑1 (MQ=255)
                                                                                                                                           ccagccagGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCTTTCg            >  1:2425192/1‑150 (MQ=255)
                                                                                                                                             agccagGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCTTTCGtt          <  1:4848696/150‑1 (MQ=255)
                                                                                                                                             agccagGTAGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCTTTCGtt          <  1:384794/150‑1 (MQ=255)
                                                                                                                                                 agGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCTTTCGTTACCg      >  1:37506/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCCCTTCCCGCAATGCAATCAAATGGCTGGCGTAACGACAGGCCTGATTAAGATCGTGCAGCACCGCCGCCAGGGTATAGCCTTTTTCGCGGTTCAGTTCGCTTAACAACTCCAGCAAATCAATCTGATGACTGATATCCAGCCAGGTCGTGGGTTCGTCGAGCAGCATAATTGCCGTTTCCTGGGCCAGCACCATCGCGATCCACGCTCGCTGGCGTTGTCCGCCAGAAAGGGTATCCACGCTTTGATCTGCCAGATGAGTTATTCCCGTGGCCTGCATTGCTTTCGTTACCGCTTC  >  CP009273/614994‑615292

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: