Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 607,931 G→A 91.5% 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,9310GA91.5% 699.3 / 38.1 226S7F (TCC→TTC) fepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (11/8);  new base A (106/100);  total (117/109)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 233 total aligned reads displayed.

CTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTCC  >  CP009273/607785‑608079
                                                                                                                                                  |                                                                                                                                                    
cTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGaaat                                                                                                                                                   >  1:3177598/1‑150 (MQ=255)
     gTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtc                                                                                                                                              <  2:4663690/150‑1 (MQ=255)
      tGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAAtct                                                                                                                                             <  1:1714795/150‑1 (MQ=255)
      tGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtct                                                                                                                                             <  2:2989836/150‑1 (MQ=255)
       gATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtctt                                                                                                                                            <  2:5126888/150‑1 (MQ=255)
       gATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAtctt                                                                                                                                            <  2:2549316/150‑1 (MQ=255)
           gTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTg                                                                                                                                        >  2:3489584/1‑150 (MQ=255)
               aaaCGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCa                                                                                                                                    >  2:4004265/1‑150 (MQ=255)
                  cGCCAGGTGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTg                                                                                                                                 <  2:525946/150‑1 (MQ=255)
                  cGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTg                                                                                                                                 <  2:526113/150‑1 (MQ=255)
                     cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                              >  1:3387958/1‑150 (MQ=255)
                     cAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAAGAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAACTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGttt                                                                                                                              <  2:2251402/150‑1 (MQ=255)
                        gcgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAACTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                           <  2:4421241/150‑1 (MQ=255)
                        gcgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                           <  1:3964772/150‑1 (MQ=255)
                        gcgcCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTAt                                                                                                                           <  2:5002237/150‑1 (MQ=255)
                           ccTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTcc                                                                                                                        <  1:1909503/150‑1 (MQ=255)
                           ccTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTcc                                                                                                                        <  2:2641580/150‑1 (MQ=255)
                             tGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTACCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTg                                                                                                                      >  2:1816555/1‑150 (MQ=255)
                              gTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGc                                                                                                                     >  2:3058672/1‑150 (MQ=255)
                              gTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGc                                                                                                                     >  1:4453048/1‑150 (MQ=255)
                                aaGGTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAt                                                                                                                   <  1:1243810/150‑1 (MQ=255)
                                  gTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAttt                                                                                                                 <  2:3725407/150‑1 (MQ=255)
                                  gTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCGACCCCATAAATCCCCAGATTGACCAACACGGCCAGGAAAGGAATCTTCTTGTTCATTGTTTTATTCCTGCAttt                                                                                                                 <  2:2281244/150‑1 (MQ=255)
                                   ttCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCAtttt                                                                                                                <  2:1629880/150‑1 (MQ=255)
                                     cTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTg                                                                                                              <  2:5048379/150‑1 (MQ=255)
                                        cTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCa                                                                                                           >  2:3574192/1‑150 (MQ=255)
                                         tCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCAc                                                                                                          >  1:5040645/1‑150 (MQ=255)
                                            gcggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGaa                                                                                                       >  2:2172380/1‑150 (MQ=255)
                                              ggcggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAAtt                                                                                                     >  1:1933660/1‑150 (MQ=255)
                                                 ggTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCa                                                                                                  >  2:1095202/1‑150 (MQ=255)
                                                  gTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCaa                                                                                                 >  1:2692181/1‑150 (MQ=255)
                                                   tAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAAc                                                                                                <  2:1978920/150‑1 (MQ=255)
                                                    aaCGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACt                                                                                               <  1:2296740/150‑1 (MQ=255)
                                                       gACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTc                                                                                            <  2:4587891/150‑1 (MQ=255)
                                                        aCAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCg                                                                                           <  2:2701601/150‑1 (MQ=255)
                                                        aCAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCg                                                                                           <  2:1896157/150‑1 (MQ=255)
                                                         cAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCgg                                                                                          >  1:2640087/1‑150 (MQ=255)
                                                              gTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATg                                                                                     <  2:2557384/150‑1 (MQ=255)
                                                                aTCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgt                                                                                   <  1:4630399/150‑1 (MQ=255)
                                                                 tCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtc                                                                                  >  1:1791427/1‑150 (MQ=255)
                                                                    tCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGgtcgtc                                                                               >  1:1609349/1‑150 (MQ=255)
                                                                      aTGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGtcat                                                                             >  1:1667699/1‑150 (MQ=255)
                                                                          gAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAac                                                                         <  1:435559/150‑1 (MQ=255)
                                                                              cAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacga                                                                     <  1:69581/150‑1 (MQ=255)
                                                                               aGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                                    <  2:1447113/150‑1 (MQ=255)
                                                                                ggAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgacg                                                                   >  1:809258/1‑150 (MQ=255)
                                                                                  aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                                 >  1:3484095/1‑150 (MQ=255)
                                                                                  aGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                                 >  2:240638/1‑150 (MQ=255)
                                                                                   gTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCAt                                                                >  1:3272885/1‑150 (MQ=255)
                                                                                   gTATCGGTCGGCTCTTGTGCCAGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCAt                                                                >  1:4979066/1‑150 (MQ=255)
                                                                                     aTCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCTACACGACGCATcc                                                              >  2:3285530/1‑150 (MQ=255)
                                                                                       cggtcggCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCg                                                            >  1:3890830/1‑150 (MQ=255)
                                                                                         gtcggCTCTTGTGCCTGCGCTACCCCATACATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCGTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCt                                                          <  2:4348190/150‑1 (MQ=255)
                                                                                             gCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTCCcg                                                      >  2:2034229/1‑150 (MQ=255)
                                                                                               tctTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgcg                                                    <  1:2114461/150‑1 (MQ=255)
                                                                                               tctTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcgcg                                                    <  2:1001015/150‑1 (MQ=255)
                                                                                                ctTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGa                                                   <  2:1791427/150‑1 (MQ=255)
                                                                                                ctTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGa                                                   <  2:820734/150‑1 (MQ=255)
                                                                                                 ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaa                                                  >  1:3666113/1‑150 (MQ=255)
                                                                                                 ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaa                                                  >  1:3878955/1‑150 (MQ=255)
                                                                                                  tgtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaa                                                 <  2:4892888/150‑1 (MQ=255)
                                                                                                  tgtgCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaa                                                 >  1:4940084/1‑150 (MQ=255)
                                                                                                     gCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAAcc                                              <  2:2008994/150‑1 (MQ=255)
                                                                                                      ccTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCt                                             >  1:3997332/1‑150 (MQ=255)
                                                                                                       cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                            <  2:2338291/150‑1 (MQ=255)
                                                                                                       cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                                            <  1:605186/150‑1 (MQ=255)
                                                                                                        tGCGCTACCCCATAAATCCCCCGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                           <  2:709199/150‑1 (MQ=255)
                                                                                                        tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                           >  1:2699309/1‑150 (MQ=255)
                                                                                                        tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCttt                                           >  2:3667756/1‑150 (MQ=255)
                                                                                                         gctcTACCCCATAAATCCCCAGATTGACTAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATACCTGCATTTTTGCCACGAATAGCAACTGTCGGGCATGGTCTTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTg                                          <  1:455800/147‑1 (MQ=255)
                                                                                                          cgcTACCCCATACATCCCCAGATTGACCAACACGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGa                                         <  1:4524106/150‑1 (MQ=255)
                                                                                                          cgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAAGCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGa                                         <  2:4169347/150‑1 (MQ=255)
                                                                                                            cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATc                                       >  1:603665/1‑150 (MQ=255)
                                                                                                               ccccATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTg                                    <  1:2411439/150‑1 (MQ=255)
                                                                                                               ccccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTg                                    >  1:4430168/1‑150 (MQ=255)
                                                                                                               ccccATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTg                                    >  2:3882345/1‑150 (MQ=255)
                                                                                                                 ccATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaa                                  >  2:1122778/1‑150 (MQ=255)
                                                                                                                  cATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGGATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaaa                                 >  2:3490110/1‑150 (MQ=255)
                                                                                                                  cATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaaa                                 >  2:3490129/1‑150 (MQ=255)
                                                                                                                  cATAAATCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaaa                                 >  1:3172275/1‑150 (MQ=255)
                                                                                                                     aaaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGac                              >  2:167043/1‑150 (MQ=255)
                                                                                                                      aaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGaca                             <  1:1377348/150‑1 (MQ=255)
                                                                                                                      aaTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGaca                             <  1:2264467/150‑1 (MQ=255)
                                                                                                                       aTCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGacac                            >  1:3361995/1‑150 (MQ=255)
                                                                                                                        tCCCCAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg                           >  2:855389/1‑150 (MQ=255)
                                                                                                                            cAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt                       <  2:4454686/150‑1 (MQ=255)
                                                                                                                            cAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt                       >  2:4577134/1‑150 (MQ=255)
                                                                                                                            cAGATTGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt                       >  1:251837/1‑150 (MQ=255)
                                                                                                                                 tGACCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtgcagt                  >  1:2556310/1‑150 (MQ=255)
                                                                                                                                   aCCAACAAGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTg                <  1:1034116/150‑1 (MQ=255)
                                                                                                                                     caacaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCACCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGt              <  1:4089541/150‑1 (MQ=255)
                                                                                                                                       gcaaGGCCAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGtta            >  1:961788/2‑150 (MQ=255)
                                                                                                                                         aaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGttaat          >  2:2824760/1‑150 (MQ=255)
                                                                                                                                             ccAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTtaattaat      >  1:3468378/1‑150 (MQ=255)
                                                                                                                                              cAGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATg     <  2:2378434/150‑1 (MQ=255)
                                                                                                                                               aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt    >  1:473474/1‑150 (MQ=255)
                                                                                                                                               aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt    >  1:3223514/1‑150 (MQ=255)
                                                                                                                                               aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt    >  1:1381309/1‑150 (MQ=255)
                                                                                                                                               aGGAAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGt    <  2:2902642/150‑1 (MQ=255)
                                                                                                                                                 ggAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTcc  >  2:1521846/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGCGGTGATGGTCGAAACGCCAGGCGCCTGTAAGTTCTGCTCGGCGGCGGTAACGACAATAGTATCGTCATGTGAAACAGGAGTATCGGTCGGCTCTTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTGGTTAATTAATGTCC  >  CP009273/607785‑608079

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: