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 1,803,201 A→G V97A (GTC→GCC)  yniB ← YniB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,803,2010AG98.0% 159.0 / ‑4.0 51V97A (GTC→GCC) yniBYniB family protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (20/30);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAG  >  NZ_CP009273/1803118‑1803283
                                                                                   |                                                                                  
gAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCaaaaa                                                                              >  1:752794/1‑90 (MQ=255)
  aTGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATc                                                                            <  2:730770/90‑1 (MQ=255)
  aTGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATc                                                                            <  2:392724/90‑1 (MQ=255)
  aTGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATc                                                                            <  2:219279/90‑1 (MQ=255)
    gAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAg                                                                          <  2:413255/90‑1 (MQ=255)
        tGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATg                                                                      <  2:832426/90‑1 (MQ=255)
        tGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATg                                                                      <  2:528617/90‑1 (MQ=255)
          aTCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGaa                                                                    >  1:713508/1‑90 (MQ=255)
              tCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATc                                                                >  1:549913/1‑90 (MQ=255)
                  cacCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCAccc                                                            <  2:403550/90‑1 (MQ=255)
                  cacCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCAccc                                                            <  2:106191/90‑1 (MQ=255)
                  cacCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCAccc                                                            <  2:119159/90‑1 (MQ=255)
                  cacCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCAccc                                                            <  2:329741/90‑1 (MQ=255)
                  cacCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCAccc                                                            <  2:904331/90‑1 (MQ=255)
                  cacCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCAccc                                                            <  1:1180394/90‑1 (MQ=255)
                   acCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCa                                                           <  1:376941/90‑1 (MQ=255)
                   acCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCa                                                           <  1:119579/90‑1 (MQ=255)
                    ccTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAg                                                          <  2:651475/90‑1 (MQ=255)
                      ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGATAATCACCCAGaa                                                        >  1:1065697/1‑90 (MQ=255)
                      ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGaa                                                        >  2:37500/1‑90 (MQ=255)
                      ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGaa                                                        >  1:869025/1‑90 (MQ=255)
                      ttCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGaa                                                        >  1:261264/1‑90 (MQ=255)
                       tCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAAc                                                       >  2:1122219/1‑90 (MQ=255)
                              agaagCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGt                                                >  1:1255001/1‑90 (MQ=255)
                              agaagCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGt                                                >  1:1268088/1‑90 (MQ=255)
                                     cgcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcc                                         <  1:815424/90‑1 (MQ=255)
                                     cgcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcc                                         <  2:350237/90‑1 (MQ=255)
                                     cgcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcc                                         <  2:1453905/90‑1 (MQ=255)
                                     cgcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcc                                         <  2:774655/90‑1 (MQ=255)
                                      gcCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcca                                        <  2:1160733/90‑1 (MQ=255)
                                            cGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTa                                  >  2:1210708/1‑90 (MQ=255)
                                                 cATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAAtt                             >  2:495215/1‑90 (MQ=255)
                                                 cATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAAtt                             >  2:593001/1‑90 (MQ=255)
                                                    ccGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAAt                                                                 >  2:431063/1‑51 (MQ=255)
                                                    ccGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAAt                                                                 <  1:431063/51‑1 (MQ=255)
                                                      gggCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcca                                        <  1:1015126/74‑1 (MQ=255)
                                                      gggCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGcca                                        >  2:1015126/1‑74 (MQ=255)
                                                        gCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTc                      >  2:396437/1‑90 (MQ=255)
                                                          aCCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCgg                    <  1:286089/90‑1 (MQ=255)
                                                          aCCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCgg                    <  1:593001/90‑1 (MQ=255)
                                                          aCCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCgg                    <  1:166863/90‑1 (MQ=255)
                                                           ccGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCa                                                           >  1:366668/1‑50 (MQ=255)
                                                           ccGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCa                                                           <  2:366668/50‑1 (MQ=255)
                                                           ccGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCggg                   <  2:713508/90‑1 (MQ=255)
                                                             ggAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGa                 >  1:831446/1‑90 (MQ=255)
                                                                aTCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGtt              <  1:246581/90‑1 (MQ=255)
                                                                aTCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGtt              <  1:756425/90‑1 (MQ=255)
                                                                aTCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGtt              <  1:996843/90‑1 (MQ=255)
                                                                        aCGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGa      >  1:484043/1‑90 (MQ=255)
                                                                          gCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAAt    >  2:1357605/1‑90 (MQ=255)
                                                                            cAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAg  <  2:836321/90‑1 (MQ=255)
                                                                                   |                                                                                  
GAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAG  >  NZ_CP009273/1803118‑1803283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ 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
GCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAGATTCAAA  >  NZ_CP009273/1803108‑1803290
                                                                                             |                                                                                         
GCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAA                                                                                     >  SRR3722078.762421/1‑100 (MQ=60)
          GAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAA                                                                           >  SRR3722078.722262/1‑100 (MQ=60)
              GAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATC                                                                       >  SRR3722078.556406/1‑100 (MQ=60)
                      CTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGATAATCACCCAGAA                                                               >  SRR3722078.1081761/1‑100 (MQ=60)
                      CTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAA                                                               >  SRR3722078.880663/1‑100 (MQ=60)
                      CTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAA                                                               >  SRR3722078.264452/1‑100 (MQ=60)
                       cggcagcgtcagatgtgtataagagacagGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAAC                                                              <  SRR3722078.436250/71‑1 (MQ=60)
                            CACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAAC                                                         <  SRR3722078.1198827/100‑1 (MQ=60)
                             ACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACG                                                        <  SRR3722078.120981/100‑1 (MQ=60)
                             ACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACG                                                        <  SRR3722078.381545/100‑1 (MQ=60)
                              CCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGT                                                       >  SRR3722078.1274872/1‑100 (MQ=60)
                              CCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGT                                                       >  SRR3722078.1288205/1‑100 (MQ=60)
                                               CGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAA                                      <  SRR3722078.826216/100‑1 (MQ=60)
                                                agacaGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAAT                                     <  SRR3722078.1030093/95‑1 (MQ=60)
                                                           CATCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGctgtctcttatacacatctgacgctgccga                          >  SRR3722078.371145/1‑70 (MQ=60)
                                                             TCCGGGCACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGA                        >  SRR3722078.842500/1‑100 (MQ=60)
                                                                    ACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATAC                 <  SRR3722078.289577/100‑1 (MQ=60)
                                                                    ACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATAC                 <  SRR3722078.599886/100‑1 (MQ=60)
                                                                    ACCGGAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATAC                 <  SRR3722078.168804/100‑1 (MQ=60)
                                                                        GAATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGA             >  SRR3722078.489793/1‑100 (MQ=60)
                                                                          ATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAAT           <  SRR3722078.249592/100‑1 (MQ=60)
                                                                          ATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAAT           <  SRR3722078.766110/100‑1 (MQ=60)
                                                                          ATCCTGCAACGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAAT           <  SRR3722078.1011358/100‑1 (MQ=60)
                                                                                   CGCCAGTCCGGCAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAGATTCAAA  >  SRR3722078.106356/1‑100 (MQ=60)
                                                                                             |                                                                                         
GCTTTTTCCAGAATGAGTTGATCTTCAACACCTTCCCGCAAGAAGCGCGCCTGGCGACTCATCCGGGCACCGGAATCCTGCAACGCCAGTCCGACAAAAATCAGGATGAAAATCACCCAGAACATAACGTTCACGCCACCATTAAAATTGGGCGTCGGGGAGTTATACCAGAATAGATTCAAA  >  NZ_CP009273/1803108‑1803290

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