Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F2 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,505,986 T→C intergenic (‑283/+7) BW25113_RS17530 ← / ← trpS GntR family transcriptional regulator YhfZ/tryptophan‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,505,9860TC93.6% 132.0 / 2.2 47intergenic (‑283/+7)BW25113_RS17530/trpSGntR family transcriptional regulator YhfZ/tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (32/12);  total (34/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

TTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCC  >  NZ_CP009273/3505900‑3506068
                                                                                      |                                                                                  
tttCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGt                                                                                 <  2:1434844/90‑1 (MQ=255)
  tCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGa                                                                               >  1:176891/1‑90 (MQ=255)
  tCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGa                                                                               >  1:671289/1‑90 (MQ=255)
     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAAtt                                                                            >  2:1232325/1‑90 (MQ=255)
     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAAtt                                                                            >  2:424874/1‑90 (MQ=255)
     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAAtt                                                                            >  1:765741/1‑90 (MQ=255)
       gCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTAc                                                                          >  2:260893/1‑90 (MQ=255)
          atatGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGc                                                                       >  1:255854/1‑90 (MQ=255)
            atGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCtt                                                                     >  2:399930/1‑90 (MQ=255)
            atGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCtt                                                                     >  2:1580676/1‑90 (MQ=255)
                 cATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCca                                                                >  2:713656/1‑90 (MQ=255)
                    tACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACaa                                                             <  1:164257/90‑1 (MQ=255)
                      cTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGt                                                                                 <  1:1054788/68‑1 (MQ=255)
                      cTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGt                                                                                 >  2:1054788/1‑68 (MQ=255)
                      cTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACaaaa                                                           >  2:1389478/1‑90 (MQ=255)
                        ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAAcc                                                         >  2:905844/1‑90 (MQ=255)
                        ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAAcc                                                         >  2:598861/1‑90 (MQ=255)
                          gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCaa                                                       >  1:1357771/1‑90 (MQ=255)
                           tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAAt                                                      <  1:563888/90‑1 (MQ=255)
                           tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAAt                                                      <  2:362285/90‑1 (MQ=255)
                           tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAAt                                                      <  1:239555/90‑1 (MQ=255)
                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCg                                                    >  2:986009/1‑90 (MQ=255)
                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCg                                                    <  2:840929/90‑1 (MQ=255)
                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCt                                                  >  2:778111/1‑90 (MQ=255)
                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCt                                                  >  2:116402/1‑90 (MQ=255)
                                      aGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTaca                                           >  2:1514158/1‑90 (MQ=255)
                                       ggCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTacac                                          >  1:1077450/1‑90 (MQ=255)
                                        gCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACAcc                                         <  2:1007158/90‑1 (MQ=255)
                                         cGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCg                                        >  2:1420876/1‑90 (MQ=255)
                                           ttCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCt                                      >  2:1385454/1‑90 (MQ=255)
                                           ttCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCt                                      >  1:779560/1‑90 (MQ=255)
                                            tCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCtt                                     >  1:1187208/1‑90 (MQ=255)
                                                   gcATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGt                              <  1:230699/90‑1 (MQ=255)
                                                         ggCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGa                        >  1:513729/1‑90 (MQ=255)
                                                          gCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGaa                       <  1:778111/90‑1 (MQ=255)
                                                             tGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCg                    >  2:810329/1‑90 (MQ=255)
                                                                   aaGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcg              >  2:1307669/1‑90 (MQ=255)
                                                                   aaGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcg              >  1:682256/1‑90 (MQ=255)
                                                                   aaGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcg              >  1:808590/1‑90 (MQ=255)
                                                                    aGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             <  1:1499360/90‑1 (MQ=255)
                                                                    aGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             <  2:699210/90‑1 (MQ=255)
                                                                        cAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACAGCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGc         >  2:1195462/1‑90 (MQ=255)
                                                                        cAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGc         >  1:318379/1‑90 (MQ=255)
                                                                          aTTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             >  2:1421029/1‑84 (MQ=255)
                                                                          aTTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             <  1:1421029/84‑1 (MQ=255)
                                                                             tGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCttttt    >  2:730502/1‑90 (MQ=255)
                                                                               ccAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTcc  >  2:992969/1‑90 (MQ=255)
                                                                                      |                                                                                  
TTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCC  >  NZ_CP009273/3505900‑3506068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTC  >  NZ_CP009273/3505892‑3506067
                                                                                              |                                                                                 
GGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGA                                                                              >  SRR3722078.678954/1‑100 (MQ=60)
GGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGA                                                                              >  SRR3722078.178964/1‑100 (MQ=60)
   CTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATT                                                                           >  SRR3722078.775608/1‑100 (MQ=60)
        tataagagacagATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGC                                                                      <  SRR3722078.1070630/88‑1 (MQ=60)
        TTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGC                                                                      >  SRR3722078.258984/1‑100 (MQ=60)
                        GCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAA                                                      >  SRR3722078.1379570/1‑100 (MQ=60)
                            TACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGC                                                  <  SRR3722078.166168/100‑1 (MQ=60)
                                   TAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTAC                                           <  SRR3722078.570502/100‑1 (MQ=60)
                                   TAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTAC                                           <  SRR3722078.242482/100‑1 (MQ=60)
                                     GGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACAC                                         >  SRR3722078.1093742/1‑100 (MQ=60)
                                         GGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCT                                     >  SRR3722078.789679/1‑100 (MQ=60)
                                          GATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTT                                    >  SRR3722078.1205760/1‑100 (MQ=60)
                                                       CGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGA                       >  SRR3722078.519809/1‑100 (MQ=60)
                                                           GCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCG                   <  SRR3722078.233496/100‑1 (MQ=60)
                                                                 GGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCG             >  SRR3722078.690212/1‑100 (MQ=60)
                                                                 GGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCG             >  SRR3722078.819266/1‑100 (MQ=60)
                                                                  GCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGC            <  SRR3722078.788201/100‑1 (MQ=60)
                                                                      GAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGC        >  SRR3722078.322279/1‑100 (MQ=60)
                                                                            AGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTC  <  SRR3722078.1523811/100‑1 (MQ=60)
                                                                            AGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTC  <  SRR3722078.1444109/100‑1 (MQ=60)
                                                                                              |                                                                                 
GGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTC  >  NZ_CP009273/3505892‑3506067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: