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 2,057,664 T→C G95G (GGA→GGG cobT ← nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,057,6640TC100.0% 137.9 / NA 39G95G (GGA→GGGcobTnicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/18);  total (21/18)

GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGC  >  NZ_CP009273/2057578‑2057747
                                                                                        |                                                                                   
gggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCa                                                                                  >  2:569363/1‑90 (MQ=255)
     aGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGt                                                                             >  1:566651/1‑90 (MQ=255)
      ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTc                                                                            >  1:417686/1‑90 (MQ=255)
      ggCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTc                                                                            >  2:344790/1‑90 (MQ=255)
        cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCat                                                                          <  2:832848/90‑1 (MQ=255)
        cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCat                                                                          <  1:54828/90‑1 (MQ=255)
          cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCatat                                                                        >  1:129499/1‑90 (MQ=255)
          cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCatat                                                                        >  1:168981/1‑90 (MQ=255)
            gcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATAttt                                                                      >  1:755708/1‑90 (MQ=255)
                tATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGGTGTCCCACGGGTCATATTTTCAg                                                                  >  2:214553/1‑90 (MQ=255)
                      tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                            >  2:162628/1‑90 (MQ=255)
                      tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                            >  2:540400/1‑90 (MQ=255)
                      tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                            >  1:562532/1‑90 (MQ=255)
                         cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                         >  1:437908/1‑90 (MQ=255)
                         cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                         >  2:593146/1‑90 (MQ=255)
                         cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGGCATATTTTCAGCCTGTATgg                                                         >  2:886171/1‑90 (MQ=255)
                           aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                         >  1:285062/1‑90 (MQ=255)
                           aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                       >  2:810317/1‑90 (MQ=255)
                           aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                       >  2:45623/1‑90 (MQ=255)
                               cAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCa                                                   >  1:22549/1‑90 (MQ=255)
                                 aTTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                 <  1:683449/90‑1 (MQ=255)
                                       tGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                           <  1:886171/90‑1 (MQ=255)
                                        ggACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                 >  1:11666/1‑83 (MQ=255)
                                        ggACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAct                                                 <  2:11666/83‑1 (MQ=255)
                                             ttAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGaa                                     <  1:569363/90‑1 (MQ=255)
                                                 cgcCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAg                                 <  2:891974/90‑1 (MQ=255)
                                                    ccGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGa                              <  1:569723/90‑1 (MQ=255)
                                                      gCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                            >  2:753430/1‑90 (MQ=255)
                                                        ttGTTCTGCCAGCACACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                            <  2:285062/90‑1 (MQ=255)
                                                         tGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCt                         <  1:665596/90‑1 (MQ=255)
                                                          gTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                        <  1:110562/90‑1 (MQ=255)
                                                                   g‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGa               <  2:553082/90‑1 (MQ=255)
                                                                        cacacacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGc            <  1:450128/90‑1 (MQ=255)
                                                                               gCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc     <  1:881547/90‑1 (MQ=255)
                                                                               gCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATc     <  1:358970/90‑1 (MQ=255)
                                                                                  ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                        <  1:682157/68‑1 (MQ=255)
                                                                                  ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                        >  2:682157/1‑68 (MQ=255)
                                                                                  ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc  <  2:174463/90‑1 (MQ=255)
                                                                                  ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc  <  1:810317/90‑1 (MQ=255)
                                                                                        |                                                                                   
GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGC  >  NZ_CP009273/2057578‑2057747

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTAC  >  NZ_CP009273/2057568‑2057762
                                                                                                  |                                                                                                  
GATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCA                                                                                                 <  SRR3722072.920741/100‑1 (MQ=60)
    AGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGG                                                                                             <  SRR3722072.307583/100‑1 (MQ=60)
     GCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGT                                                                                            >  SRR3722072.576795/1‑100 (MQ=60)
      CCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTC                                                                                           >  SRR3722072.425254/1‑100 (MQ=60)
          GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATAT                                                                                       >  SRR3722072.131835/1‑100 (MQ=60)
          GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATAT                                                                                       >  SRR3722072.171902/1‑100 (MQ=60)
            GATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTT                                                                                     >  SRR3722072.769047/1‑100 (MQ=60)
                  CTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCC                                                                               <  SRR3722072.55820/100‑1 (MQ=60)
                      GCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTA                                                                           >  SRR3722072.572621/1‑100 (MQ=60)
                         GTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGG                                                                        >  SRR3722072.445894/1‑100 (MQ=60)
                           ATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGG                                                                        >  SRR3722072.289949/1‑100 (MQ=60)
                               ATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCA                                                                  >  SRR3722072.22946/1‑100 (MQ=60)
                                        TCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTG                                                         >  SRR3722072.11887/1‑100 (MQ=60)
                                           ATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAG                                                      <  SRR3722072.695680/100‑1 (MQ=60)
                                                 TGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAG                                                <  SRR3722072.901418/100‑1 (MQ=60)
                                                       TTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCC                                          <  SRR3722072.579543/100‑1 (MQ=60)
                                                              CCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTC                                   <  SRR3722072.579907/100‑1 (MQ=60)
                                                                   TGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATctgtctcttatacacatctccgagcccacgagacgtccgattatctcgtatgccgtc                              >  SRR3722072.336272/1‑43 (MQ=60)
                                                                   TGTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGA                              <  SRR3722072.677497/100‑1 (MQ=60)
                                                                    GTTCTGCCAG‑‑CACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGAC                             <  SRR3722072.112607/100‑1 (MQ=60)
                                                                      tataagagacagCACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGAC                             <  SRR3722072.694362/88‑1 (MQ=60)
                                                                                  CACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGC                 <  SRR3722072.458372/100‑1 (MQ=60)
                                                                                         GCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATA          <  SRR3722072.896718/100‑1 (MQ=60)
                                                                                         GCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATA          <  SRR3722072.365323/100‑1 (MQ=60)
                                                                                            GGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACC       <  SRR3722072.824442/100‑1 (MQ=60)
                                                                                                 TCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTAC  <  SRR3722072.766733/100‑1 (MQ=60)
                                                                                                  |                                                                                                  
GATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAG‑‑CACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTAC  >  NZ_CP009273/2057568‑2057762

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