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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 945,384 (T)5→4 coding (261/591 nt) ycaK → NAD(P)H oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273945,3800T.97.6% 179.3 / ‑2.1 41coding (257/591 nt)ycaKNAD(P)H oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (14/26);  minor base G (1/0);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATG  >  NZ_CP009273/945295‑945466
                                                                                     |                                                                                      
gACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGttttcc                                                                                   >  2:790628/1‑89 (MQ=255)
   tGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCctct                                                                               <  1:1073765/91‑1 (MQ=255)
        gaaTATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCtggt                                                                          <  1:1112764/91‑1 (MQ=255)
        gaaTATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCtggt                                                                          <  2:210845/91‑1 (MQ=255)
          atatGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCtggtgg                                                                        >  2:524840/1‑91 (MQ=255)
                  aaaCGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGGGGTACCGCtt                                                                >  2:52993/1‑91 (MQ=255)
                        tATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGtttt                                                                                     <  2:277571/65‑1 (MQ=255)
                        tATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGtttt                                                                                     >  1:277571/1‑65 (MQ=255)
                        tATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGc                                                          >  1:605596/1‑91 (MQ=255)
                        tATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGc                                                          <  2:804490/91‑1 (MQ=255)
                         aTTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACCTGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCa                                                         <  1:1051516/91‑1 (MQ=255)
                          ttCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCaa                                                        <  2:1091281/91‑1 (MQ=255)
                          ttCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCaa                                                        <  1:59027/91‑1 (MQ=255)
                              ccAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCt                                                    >  1:645546/1‑91 (MQ=255)
                                      tcatcaGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGat                                            <  1:239664/91‑1 (MQ=255)
                                      tcatcaGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGat                                            <  2:755126/91‑1 (MQ=255)
                                      tcatcaGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGat                                            >  1:830394/1‑91 (MQ=255)
                                        atcaGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCtggt                                                                          <  1:953164/59‑1 (MQ=255)
                                        atcaGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCtggt                                                                          >  2:953164/1‑59 (MQ=255)
                                             cTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCaa                                                        <  1:262006/72‑1 (MQ=255)
                                             cTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCaa                                                        >  2:262006/1‑72 (MQ=255)
                                              tttATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGAc                                    <  2:114022/91‑1 (MQ=255)
                                                tATTCAGAGCTTCATGAACATGACACGTTAGTGGTAGTTTTCCTCTCTGGTGGTGCAGCTTCCCGGCAATGCTAAAAGGATAAATTGACag                                  <  2:57450/91‑1 (MQ=255)
                                                tATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACag                                  <  1:599092/91‑1 (MQ=255)
                                                tATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACag                                  <  2:662671/91‑1 (MQ=255)
                                                tATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACag                                  <  2:605596/91‑1 (MQ=255)
                                                 aTTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACaga                                 >  2:336080/1‑91 (MQ=255)
                                                  ttCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACagag                                <  2:712835/91‑1 (MQ=255)
                                                      gagCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATg                            >  1:954112/1‑91 (MQ=255)
                                                           gctTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGat                                            <  2:127981/70‑1 (MQ=255)
                                                           gctTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGat                                            >  1:127981/1‑70 (MQ=255)
                                                           gctTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGaata                       >  2:116172/1‑91 (MQ=255)
                                                           gctTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGaata                       >  1:17105/1‑91 (MQ=255)
                                                           gctTGAACATGACACGTTAGTGGGGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGaata                       >  1:149831/1‑91 (MQ=255)
                                                            ctTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGaataa                      <  2:875264/91‑1 (MQ=255)
                                                                 aCATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGggc                 <  2:217247/91‑1 (MQ=255)
                                                                 aCATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGggc                 <  1:336080/91‑1 (MQ=255)
                                                                 aCATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGggc                 <  1:219476/91‑1 (MQ=255)
                                                                 aCATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGggc                 <  2:612756/91‑1 (MQ=255)
                                                                   aTGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGggctg               <  1:546936/91‑1 (MQ=255)
                                                                     gACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACag                                  <  1:1004556/70‑1 (MQ=255)
                                                                     gACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACag                                  >  2:1004556/1‑70 (MQ=255)
                                                                     gACACGTTAGTGGTAGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGggctggc             >  1:718311/1‑91 (MQ=255)
                                                                        acGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTa          <  2:502008/91‑1 (MQ=255)
                                                                               tggtggTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGAt   <  1:442184/91‑1 (MQ=255)
                                                                                ggtggtTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATg  >  1:567898/1‑91 (MQ=255)
                                                                                ggtggtTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATg  <  2:91876/91‑1 (MQ=255)
                                                                                     |                                                                                      
GACTGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATG  >  NZ_CP009273/945295‑945466

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTC  >  NZ_CP009273/945298‑945483
                                                                                      |                                                                                                   
TGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTAC                                                                                      <  SRR3722125.1097943/101‑1 (MQ=60)
     cggagatgtgtataagagacagCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTT                                                                                 <  SRR3722125.974395/79‑1 (MQ=60)
     GAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTT                                                                                 <  SRR3722125.1137865/101‑1 (MQ=60)
           GGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTctgtctcttatacaca                                                                           >  SRR3722125.283185/1‑85 (MQ=60)
           GGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGC                                                                           >  SRR3722125.618641/1‑101 (MQ=60)
                 ACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCT                                                                     >  SRR3722125.659470/1‑101 (MQ=60)
                      ATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACCTGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAG                                                                <  SRR3722125.1075136/101‑1 (MQ=60)
                       aagagacagGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGG                                                               <  SRR3722125.267283/92‑1 (MQ=60)
                       TTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGG                                                               <  SRR3722125.60116/101‑1 (MQ=60)
                         CTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGAT                                                             >  SRR3722125.848632/1‑101 (MQ=60)
                                   TCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAG                                                   <  SRR3722125.244467/101‑1 (MQ=60)
                                         GCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATG                                             >  SRR3722125.975365/1‑101 (MQ=60)
                                             tcggagatgtgtataagagacagCACGTTAGTCGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAAT                                         <  SRR3722125.1021166/78‑1 (MQ=60)
                                             ataagagacagGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAAT                                         <  SRR3722125.1027019/90‑1 (MQ=60)
                                             TATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAAT                                         <  SRR3722125.611985/101‑1 (MQ=60)
                                              ATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGctgtctcttat                                        >  SRR3722125.130388/1‑90 (MQ=60)
                                              ATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATA                                        >  SRR3722125.17414/1‑101 (MQ=60)
                                              ATTCAGAGCTGCTTGAACATGACACGTTAGTGGGGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATA                                        >  SRR3722125.152679/1‑101 (MQ=60)
                                                        GCTTGAACATGACACGTTAGTGGTAGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGC                              >  SRR3722125.733814/1‑101 (MQ=60)
                                                              ACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGG                        <  SRR3722125.343048/101‑1 (MQ=60)
                                                              ACATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGG                        <  SRR3722125.223809/101‑1 (MQ=60)
                                                                ATGACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAG                      <  SRR3722125.558649/101‑1 (MQ=60)
                                                                   ACACGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATG                   >  SRR3722125.580042/1‑101 (MQ=60)
                                                                      CGTTAGTGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGC                >  SRR3722125.257567/1‑101 (MQ=60)
                                                                            TGGTGGTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAAT          <  SRR3722125.451517/101‑1 (MQ=60)
                                                                                   TTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCAT    <  SRR3722125.543790/101‑1 (MQ=60)
                                                                                   TTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCAT    <  SRR3722125.883256/101‑1 (MQ=60)
                                                                                     TTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTC  <  SRR3722125.597798/101‑1 (MQ=60)
                                                                                      |                                                                                                   
TGGAAGAATATGGATAAACGTTATTCTCCAGAGGTTCATCAGCTTTATTCAGAGCTGCTTGAACATGACACGTTAGTGGTGGTTTTTCCTCTCTGGTGGTACAGCTTCCCGGCAATGCTAAAAGGATATATTGACAGAGTATGGAATAATGGGCTGGCTTATGGAGATGGGCACAAATTACCATTC  >  NZ_CP009273/945298‑945483

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