Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 607,931 G→A 100% S7F (TCC→TTC)  fepA ← iron‑enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273607,9310GA97.2% 899.1 / 10.7 249S7F (TCC→TTC) fepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base A (140/102);  total (143/106)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 271 total aligned reads displayed.

ATCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATG  >  CP009273/607783‑608076
                                                                                                                                                    |                                                                                                                                                 
  cTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGgaat                                                                                                                                                <  2:1900821/150‑1 (MQ=255)
   tGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTAATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatt                                                                                                                                               <  2:408941/150‑2 (MQ=255)
    gCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatga                                                                                                                                              <  1:924384/150‑1 (MQ=255)
    gCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatga                                                                                                                                              <  2:2504218/150‑1 (MQ=255)
    gCGGTGATGGTCGAAACGCCAGGAGCCTGTAAGTTCTGCTCGTCGGCGGTAAAGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTAGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatga                                                                                                                                              <  1:3039705/150‑1 (MQ=255)
      ggTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaat                                                                                                                                            >  1:1206744/1‑150 (MQ=255)
      ggTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaat                                                                                                                                            >  1:2927487/1‑150 (MQ=255)
      ggTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaatgaat                                                                                                                                            >  2:3523463/1‑150 (MQ=255)
        tGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtct                                                                                                                                          >  2:3991698/1‑150 (MQ=255)
           tGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtcttct                                                                                                                                       >  1:4148805/1‑150 (MQ=255)
           tGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtcttct                                                                                                                                       <  1:1007974/150‑1 (MQ=255)
           tGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtcttct                                                                                                                                       <  2:4905756/150‑1 (MQ=255)
             gTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTg                                                                                                                                     >  1:4593150/1‑150 (MQ=255)
             gTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTg                                                                                                                                     >  1:1155115/1‑150 (MQ=255)
             gTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTg                                                                                                                                     <  2:3939097/150‑1 (MQ=255)
              tCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGt                                                                                                                                    <  2:3372624/150‑1 (MQ=255)
               cGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCATGAAATGAATCTTCTTGtt                                                                                                                                   <  1:3428891/150‑1 (MQ=255)
                gAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTc                                                                                                                                  >  1:2853925/1‑150 (MQ=255)
                 aaaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCa                                                                                                                                 >  2:2410879/1‑150 (MQ=255)
                  aaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAt                                                                                                                                >  2:1487747/1‑150 (MQ=255)
                   aCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCAtt                                                                                                                               >  1:183061/1‑150 (MQ=255)
                    cGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTg                                                                                                                              <  1:537805/150‑1 (MQ=255)
                       cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                           <  2:1651011/150‑1 (MQ=255)
                       cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                           >  1:443184/1‑150 (MQ=255)
                       cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                           >  1:2493/1‑150 (MQ=255)
                         ggCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTa                                                                                                                         <  1:133941/150‑1 (MQ=255)
                           cgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTAtt                                                                                                                       >  1:1680527/1‑150 (MQ=255)
                            gcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTc                                                                                                                      >  1:837184/1‑150 (MQ=255)
                             ccTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTcc                                                                                                                     <  2:1732739/150‑1 (MQ=255)
                                ggcaGTTCGGCTCGGCGGCGGGAACGACACGAGGATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAAGCCCCAGATTGACCACCCCGGCCAGGAAATGAATCTTCTGGTTCATTGTTTTATTCCTGc                                                                                                                  <  2:3248752/147‑1 (MQ=255)
                                 tAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCa                                                                                                                 <  2:1711736/150‑1 (MQ=255)
                                  aaGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAt                                                                                                                >  2:2227087/1‑150 (MQ=255)
                                  aaGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAt                                                                                                                >  2:752751/1‑150 (MQ=255)
                                       cTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTg                                                                                                           >  2:5122211/1‑150 (MQ=255)
                                       cTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTg                                                                                                           >  2:287419/1‑150 (MQ=255)
                                        tGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGc                                                                                                          >  2:3105778/1‑150 (MQ=255)
                                          cTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCa                                                                                                        <  2:2187898/150‑1 (MQ=255)
                                            cggcggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACg                                                                                                      >  2:3462143/1‑150 (MQ=255)
                                              gcggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGaa                                                                                                    <  2:4393809/150‑1 (MQ=255)
                                              gcggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGaa                                                                                                    >  2:3786984/1‑150 (MQ=255)
                                               cggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAt                                                                                                   <  1:3263670/150‑1 (MQ=255)
                                                    gTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCaa                                                                                              <  1:4807382/150‑1 (MQ=255)
                                                      aaCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACt                                                                                            >  1:950714/1‑150 (MQ=255)
                                                      aaCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACt                                                                                            <  1:3072344/150‑1 (MQ=255)
                                                         gACTATAGTATCGTCATGTGAAACAGGAGTTTCGGTCGGCTCTTGGGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCATGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTc                                                                                         <  2:1346945/150‑1 (MQ=255)
                                                           cAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCgg                                                                                       <  2:1872986/150‑1 (MQ=255)
                                                           cAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCgg                                                                                       >  2:2355288/1‑150 (MQ=255)
                                                            aaTAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCggg                                                                                      >  2:1960484/1‑150 (MQ=255)
                                                                gTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATACCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATg                                                                                  >  1:3191257/1‑150 (MQ=255)
                                                                 tATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATgg                                                                                 <  1:3720132/150‑1 (MQ=255)
                                                                 tATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATgg                                                                                 >  1:2003578/1‑150 (MQ=255)
                                                                 tATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATgg                                                                                 <  1:2604766/150‑1 (MQ=255)
                                                                 tATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATgg                                                                                 >  2:1901970/1‑150 (MQ=255)
                                                                  aTCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgt                                                                                >  1:2775393/1‑150 (MQ=255)
                                                                  aTCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgt                                                                                >  1:1290943/1‑150 (MQ=255)
                                                                    cGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcg                                                                              >  1:3605743/1‑150 (MQ=255)
                                                                    cGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcg                                                                              >  1:3679391/1‑150 (MQ=255)
                                                                    cGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcg                                                                              <  2:1126135/150‑1 (MQ=255)
                                                                    cGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcg                                                                              >  1:4983471/1‑150 (MQ=255)
                                                                    cGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcg                                                                              <  1:4257591/150‑1 (MQ=255)
                                                                    cGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGACACGAATTGCAACTGTCGGGCATGgtcg                                                                              >  1:1496353/1‑150 (MQ=255)
                                                                          gtgAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGtcatca                                                                        >  2:202803/1‑150 (MQ=255)
                                                                            gAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAac                                                                      >  2:4702095/1‑150 (MQ=255)
                                                                             aaaCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                                     >  1:5083952/1‑150 (MQ=255)
                                                                             aaaCAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                                     <  2:4151688/150‑1 (MQ=255)
                                                                                cAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacga                                                                  >  2:3269145/1‑150 (MQ=255)
                                                                                 aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                 <  1:4304336/150‑1 (MQ=255)
                                                                                   gAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGc                                                               >  2:280462/1‑150 (MQ=255)
                                                                                    aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                              <  1:1914580/150‑1 (MQ=255)
                                                                                      tATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATc                                                            >  2:502457/1‑150 (MQ=255)
                                                                                         cggtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCg                                                         <  1:4743532/150‑1 (MQ=255)
                                                                                          ggtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGc                                                        >  1:1198418/1‑150 (MQ=255)
                                                                                            tcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTa                                                      >  1:4369178/1‑150 (MQ=255)
                                                                                              ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                    <  1:2287900/150‑1 (MQ=255)
                                                                                              ggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                                    <  2:4359757/150‑1 (MQ=255)
                                                                                                ctctTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTACATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgc                                                  >  1:1455032/1‑150 (MQ=255)
                                                                                                 tctTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgcg                                                 >  2:4845136/1‑150 (MQ=255)
                                                                                                   ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaa                                               >  2:1434476/1‑150 (MQ=255)
                                                                                                      tgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAAc                                            <  1:263245/150‑1 (MQ=255)
                                                                                                         cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                         >  2:181354/1‑150 (MQ=255)
                                                                                                         cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                         >  2:2631216/1‑150 (MQ=255)
                                                                                                         cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                         <  1:1863746/150‑1 (MQ=255)
                                                                                                            cgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGa                                      >  1:2133988/1‑150 (MQ=255)
                                                                                                             gcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGAt                                     >  1:486536/1‑150 (MQ=255)
                                                                                                                 ccccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTg                                 <  2:2094028/150‑1 (MQ=255)
                                                                                                                     aTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                             <  2:1112138/150‑1 (MQ=255)
                                                                                                                     aTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                             >  2:4508296/1‑150 (MQ=255)
                                                                                                                     aTAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                             >  1:60123/1‑150 (MQ=255)
                                                                                                                        aaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGaca                          >  2:1783491/1‑150 (MQ=255)
                                                                                                                         aTCCCCTGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCAGCTACCGCGAAAACCTTTGATCCTGAAAGacac                         <  1:3796791/150‑1 (MQ=255)
                                                                                                                         aTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGacac                         >  2:4708249/1‑150 (MQ=255)
                                                                                                                           ccccAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgc                       >  1:4740354/1‑150 (MQ=255)
                                                                                                                             ccAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcag                     <  2:2421682/150‑1 (MQ=255)
                                                                                                                               aGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtg                   >  1:4091354/1‑150 (MQ=255)
                                                                                                                                        aacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt          >  1:4064147/1‑150 (MQ=255)
                                                                                                                                        aacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtt          >  2:1196058/1‑150 (MQ=255)
                                                                                                                                         acaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtta         <  2:1412444/150‑1 (MQ=255)
                                                                                                                                            aGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaatt      >  1:2213422/1‑150 (MQ=255)
                                                                                                                                            aGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaatt      <  2:4773462/150‑1 (MQ=255)
                                                                                                                                               ccAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTtaattaat   >  2:3229205/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATG  >  CP009273/607783‑608076

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: