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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,836,838 (C)5→4 coding (312/1335 nt) adeQ ← adenine permease AdeQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,836,8340C.97.7% 175.3 / ‑2.0 44coding (316/1335 nt)adeQadenine permease AdeQ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (16/27);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

GGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAATAC  >  NZ_CP009273/3836750‑3836923
                                                                                    |                                                                                         
ggATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGccccac                                                                                      >  1:600809/1‑88 (MQ=255)
ggATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGccccac                                                                                      >  2:783866/1‑88 (MQ=255)
ggATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGccccac                                                                                      >  2:567041/1‑88 (MQ=255)
ggATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGccccac                                                                                      >  1:743039/1‑88 (MQ=255)
        atcAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTTCGCCCAGCCCGGTCTGCCAGGAGATGCCCACCGTCccacgacctcga                                                                             <  2:469993/88‑1 (MQ=255)
        aaCAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCccacgaccacga                                                                             <  2:889031/90‑1 (MQ=255)
        aaCAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCccacgaccacga                                                                             <  1:567041/90‑1 (MQ=255)
        aaCAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCccacgaccacga                                                                             <  2:343468/90‑1 (MQ=255)
        aaCAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCccacgaccacga                                                                             <  2:741034/90‑1 (MQ=255)
          cAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAg                                                                           <  2:829153/90‑1 (MQ=255)
           aGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAgg                                                                          >  1:713002/1‑90 (MQ=255)
            gCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGc                                                                         >  2:197791/1‑90 (MQ=255)
            gCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGc                                                                         >  1:241464/1‑90 (MQ=255)
            gCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGc                                                                         >  2:348566/1‑90 (MQ=255)
            gCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGc                                                                         >  2:775326/1‑90 (MQ=255)
               tGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCaaa                                                                      >  1:287917/1‑90 (MQ=255)
               tGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCaaa                                                                      >  2:835088/1‑90 (MQ=255)
                  gCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAaagaa                                                                   <  2:740637/90‑1 (MQ=255)
                    aaaaaTAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAgg                                                                 <  1:835088/90‑1 (MQ=255)
                      aaaTAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCg                                                               >  1:457617/1‑90 (MQ=255)
                          aGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCa                                                           <  1:487712/90‑1 (MQ=255)
                          aGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCa                                                           <  2:713002/90‑1 (MQ=255)
                           gTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAgg                                                                 >  2:400296/1‑83 (MQ=255)
                           gTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAgg                                                                 <  1:400296/83‑1 (MQ=255)
                            tCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCCACGACCACGAAGGCAAAGAAGGCGTTCAg                                                          <  1:264838/90‑1 (MQ=255)
                            tCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGc                                                         <  1:815470/90‑1 (MQ=255)
                                aCTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCa                                                     >  2:13613/1‑90 (MQ=255)
                                           gAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGcc                                          <  1:163093/90‑1 (MQ=255)
                                                 tGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGcc                                    <  1:503842/90‑1 (MQ=255)
                                                  gcgcCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCa                                   <  1:50333/90‑1 (MQ=255)
                                                  gcgcCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCa                                   <  1:383657/90‑1 (MQ=255)
                                                    gcCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACg                                 >  2:821556/1‑90 (MQ=255)
                                                     cccATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACgg                                >  1:679900/1‑90 (MQ=255)
                                                        aTCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTa                             <  2:477860/90‑1 (MQ=255)
                                                         tCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTaa                            >  1:190225/1‑90 (MQ=255)
                                                         tCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTaa                            <  1:783866/90‑1 (MQ=255)
                                                              gTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAg                       <  2:679900/90‑1 (MQ=255)
                                                               tCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGc                      <  1:253035/90‑1 (MQ=255)
                                                                  gCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTc                                                            <  2:185015/49‑1 (MQ=39)
                                                                  gCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTc                                                            >  1:185015/1‑49 (MQ=39)
                                                                    cAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAata                 <  2:203917/90‑1 (MQ=255)
                                                                       gagaTGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATcc              <  1:161543/90‑1 (MQ=255)
                                                                         gaTGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATcccc            <  2:640842/90‑1 (MQ=255)
                                                                         gaTGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATcccc            <  1:251883/90‑1 (MQ=255)
                                                                         gaTGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATcccc            <  1:225968/90‑1 (MQ=255)
                                                                          aTGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGACCCACCCCCACGGGAAAGTTAGCAAATATcccct           >  1:219673/1‑89 (MQ=255)
                                                                            gCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGc                      <  2:641401/77‑1 (MQ=255)
                                                                            gCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGc                      >  1:641401/1‑77 (MQ=255)
                                                                                    cccccACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAATAc  <  2:340826/90‑1 (MQ=255)
                                                                                    |                                                                                         
GGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAATAC  >  NZ_CP009273/3836750‑3836923

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ 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
CAGTACCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAATAC  >  NZ_CP009273/3836743‑3836923
                                                                                               |                                                                                     
CAGTACCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGAC                                                                                  <  SRR3722072.178144/100‑1 (MQ=60)
        GATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGG                                                                          >  SRR3722072.725708/1‑100 (MQ=60)
         ATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGC                                                                         >  SRR3722072.245511/1‑100 (MQ=60)
            CGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAA                                                                      >  SRR3722072.292886/1‑100 (MQ=60)
               AACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAA                                                                   <  SRR3722072.577189/100‑1 (MQ=60)
                   GCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCG                                                               >  SRR3722072.465996/1‑100 (MQ=60)
                           AAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCC                                                       <  SRR3722072.849556/100‑1 (MQ=60)
                           AAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCC                                                       <  SRR3722072.407515/100‑1 (MQ=60)
                                 AGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGC                                                 <  SRR3722072.496528/100‑1 (MQ=60)
                                   TCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCC                                                <  SRR3722072.269346/100‑1 (MQ=60)
                                   TCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCG                                               <  SRR3722072.829677/100‑1 (MQ=60)
                                                  GAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGG                                >  SRR3722072.692063/1‑100 (MQ=60)
                                                  GAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGG                                <  SRR3722072.165931/100‑1 (MQ=60)
                                                   ataatcggacgtctcgtgggctcggagatgtgtataagagacagCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGG                                <  SRR3722072.364856/56‑1 (MQ=60)
                                                      ATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAA                            >  SRR3722072.193526/1‑100 (MQ=60)
                                                        TGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGT                          <  SRR3722072.512936/100‑1 (MQ=60)
                                                         GCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTT                         <  SRR3722072.390512/100‑1 (MQ=60)
                                                         GCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTT                         <  SRR3722072.51256/100‑1 (MQ=60)
                                                               ATCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCtgtctcttatacacatctccgagcccacga                   >  SRR3722072.188229/1‑70 (MQ=60)
                                                                TCCCGGTCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAAT                  <  SRR3722072.797601/100‑1 (MQ=60)
                                                                      TCTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCC            <  SRR3722072.257321/100‑1 (MQ=60)
                                                                       CTGCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGACCCACCCCCACGGGAAAGTTAGCAAATATCCCCT           >  SRR3722072.223409/1‑100 (MQ=60)
                                                                         GCCAGGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCctg         >  SRR3722072.652841/1‑97 (MQ=60)
                                                                             GGAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAA     >  SRR3722072.928494/1‑100 (MQ=60)
                                                                              GAGATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAAT    <  SRR3722072.164357/100‑1 (MQ=60)
                                                                                GATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAATAC  <  SRR3722072.256145/100‑1 (MQ=60)
                                                                                GATGCCCATCGCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAATAC  <  SRR3722072.229788/100‑1 (MQ=60)
                                                                                               |                                                                                     
CAGTACCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCCCCACGACCACGAAGGCAAAGAAGGCGTTCAGCCCCATTGCCGGAGCCAGCGCCACGGGTAAGTTAGCAAATATCCCCATCGCAATAC  >  NZ_CP009273/3836743‑3836923

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