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,122,270 (A)6→5 coding (9/924 nt) yceM → Gfo/Idh/MocA family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,122,2650A.96.5% 111.6 / ‑1.2 29coding (4/924 nt)yceMGfo/Idh/MocA family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (10/18);  minor base G (1/0);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGAC  >  NZ_CP009273/1122179‑1122350
                                                                                      |                                                                                     
ggCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaa                                                                                    >  1:1219347/1‑90 (MQ=255)
  cccGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaat                                                                                  >  2:760530/1‑89 (MQ=255)
   ccGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaatt                                                                                 <  2:557774/90‑3 (MQ=255)
     gcgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaattac                                                                               >  2:1028886/1‑86 (MQ=255)
     gcgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaattac                                                                               >  2:914872/1‑86 (MQ=255)
     gcgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGaaaaattac                                                                               >  2:1410695/1‑86 (MQ=255)
       gTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGt                                                                            <  1:1230938/90‑1 (MQ=255)
             gCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGc                                                                      <  2:1219347/90‑1 (MQ=255)
                cTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTa                                                                   <  1:152580/90‑1 (MQ=255)
                       tCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGAt                                                            <  1:168235/90‑1 (MQ=255)
                       tCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGAt                                                            <  2:724320/90‑1 (MQ=255)
                       tCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGAt                                                            <  2:933975/90‑1 (MQ=255)
                                       gAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCa                                            >  2:1254613/1‑90 (MQ=255)
                                         aaCAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCaaa                                          >  1:145868/1‑90 (MQ=255)
                                             gCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGc                                      >  2:1415479/1‑90 (MQ=255)
                                             gCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGc                                      >  1:892268/1‑90 (MQ=255)
                                              cGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCg                                     <  1:104683/90‑1 (MQ=255)
                                                  ttGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGt                                 >  2:952460/1‑90 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:182449/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:1028886/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:1410695/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  2:1471630/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  2:145868/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTAcc                             <  1:861434/90‑1 (MQ=255)
                                                      ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTCAcc                             <  1:341321/90‑1 (MQ=255)
                                                           tGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGc                        >  1:937825/1‑90 (MQ=255)
                                                                    aTCTGGGAGATTAAAGGGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGGGGTTACCGGTGCTGGCGGCAg               >  1:1393920/1‑90 (MQ=255)
                                                                            gaTTAAAGTGCAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGAc       <  1:1518565/90‑1 (MQ=255)
                                                                              ttAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTg     >  2:672924/1‑90 (MQ=255)
                                                                              ttAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTg     >  2:1234223/1‑90 (MQ=255)
                                                                              ttAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTg     >  2:886172/1‑90 (MQ=255)
                                                                                aaaGTGAAAAATTACGTATTGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGa   <  1:914872/90‑1 (MQ=255)
                                                                                aaaGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGa   <  1:952460/90‑1 (MQ=255)
                                                                                aaaGTGAAAAATTACGTACCGGCGTAGTGGGATTAGGTGGCATTGCGCTAAAAGCGGGGTTATTGGTGCTGGCGGCAGCGTCTGACGGGa   <  1:1345517/90‑1 (MQ=255)
                                                                                 aaGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGAc  >  2:648208/1‑90 (MQ=255)
                                                                                      |                                                                                     
GGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGAC  >  NZ_CP009273/1122179‑1122350

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ 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
GATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGACGTTACAAGGAG  >  NZ_CP009273/1122171‑1122361
                                                                                                   |                                                                                           
GATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAT                                                                                             <  SRR3722078.19061/100‑1 (MQ=60)
GATTTACAGGCCCGCGCGGAAGCCCTGGAAAACGAAGTGGAAGAACTGAAAAAGCGTCTTTATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAT                                                                                             <  SRR3722078.1200044/100‑1 (MQ=60)
               GTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAG                                                                             <  SRR3722078.1250338/100‑1 (MQ=60)
                        CTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAG                                                                    <  SRR3722078.154348/100‑1 (MQ=60)
                               TCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCAT                                                             <  SRR3722078.170185/100‑1 (MQ=60)
                                       GCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAA                                                     >  SRR3722078.147562/1‑100 (MQ=60)
                                           AACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGC                                                 >  SRR3722078.904263/1‑100 (MQ=60)
                                                      CGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGG                                      <  SRR3722078.105926/100‑1 (MQ=60)
                                                         CTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGC                                   >  SRR3722078.950831/1‑100 (MQ=60)
                                                              TTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCG                              <  SRR3722078.1433586/100‑1 (MQ=60)
                                                              TTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCG                              <  SRR3722078.184612/100‑1 (MQ=60)
                                                              TTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCG                              <  SRR3722078.872958/100‑1 (MQ=60)
                                                              TTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCG                              <  SRR3722078.1044158/100‑1 (MQ=60)
                                                              TTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTCACCGGTGCTGGCG                              <  SRR3722078.345498/100‑1 (MQ=60)
                                                                  TTGCTGGCCCATCTGGGAGATTAAAGGGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGGGGTTACCGGTGCTGGCGGCAG                          >  SRR3722078.1416446/1‑100 (MQ=60)
                                                                                    GATTAAAGTGCAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGACGTTAC        <  SRR3722078.1543331/100‑1 (MQ=60)
                                                                                        AAAGTGAAAAATTACGTATTGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGACGTTACAAGG    <  SRR3722078.927289/100‑1 (MQ=60)
                                                                                        AAAGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGACGTTACAAGG    <  SRR3722078.965854/100‑1 (MQ=60)
                                                                                        AAAGTGAAAAATTACGTACCGGCGTAGTGGGATTAGGTGGCATTGCGCTAAAAGCGGGGTTATTGGTGCTGGCGGCAGCGTCTGACGGGACGTTCCAAGG    <  SRR3722078.1367113/100‑1 (MQ=60)
                                                                                          AGTGAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGACGTTACAAGGAG  <  SRR3722078.1343518/100‑1 (MQ=60)
                                                                                                   |                                                                                           
GATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTGCGCAAAAAGCGTGGTTACCGGTGCTGGCGGCAGCGTCTGACTGGACGTTACAAGGAG  >  NZ_CP009273/1122171‑1122361

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