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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,060,988 (A)6→5 intergenic (+240/‑53) yccE → / → agp YccE family protein/bifunctional glucose‑1‑phosphatase/inositol phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,060,9830A.100.0% 114.8 / NA 27intergenic (+235/‑58)yccE/agpYccE family protein/bifunctional glucose‑1‑phosphatase/inositol phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (12/15);  total (12/15)

GGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAG  >  NZ_CP009273/1060894‑1061065
                                                                                         |                                                                                  
ggCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGa                                                                                    >  2:590990/1‑90 (MQ=255)
      gTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAACCCATTTTGCGAACCTGTTCCCGGaaacagt                                                                              >  1:335274/1‑86 (MQ=255)
          gCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCata                                                                         >  1:267786/1‑90 (MQ=255)
                    ttGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCATCAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTcac                                                               >  2:392308/1‑90 (MQ=255)
                      gCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCacac                                                             <  2:631378/90‑1 (MQ=255)
                      gCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCacac                                                             <  2:503441/90‑1 (MQ=255)
                      gCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCacac                                                             <  1:806286/90‑1 (MQ=255)
                      gCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCacac                                                             <  1:1201946/90‑1 (MQ=255)
                       ccTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACAct                                                            <  2:635372/90‑1 (MQ=255)
                           tttttGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGCTCCCGGAAAAAGTCATATTTCTGTCACACTCttt                                                        <  1:1227737/90‑1 (MQ=255)
                                   gATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattg                                                <  1:751909/90‑1 (MQ=255)
                                   gATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattg                                                <  1:1213737/90‑1 (MQ=255)
                                    aTGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGtgattga                                               >  1:752291/1‑90 (MQ=255)
                                      gTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCTAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATa                                             <  1:1216654/90‑1 (MQ=255)
                                      gTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTTTGTCACACTCTTTAGTGATTGATa                                             <  1:762800/90‑1 (MQ=255)
                                      gTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATa                                             <  1:433681/90‑1 (MQ=255)
                                      gTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATa                                             <  1:979245/90‑1 (MQ=255)
                                           taGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACaaa                                        >  1:8753/1‑90 (MQ=255)
                                               aaaaCAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAagag                                    >  1:943262/1‑90 (MQ=255)
                                                aaaCAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAgg                                   >  2:748131/1‑90 (MQ=255)
                                                    aGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGcc                               <  2:8753/90‑1 (MQ=255)
                                                    aGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGcc                               <  1:706651/90‑1 (MQ=255)
                                                     ggCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCa                              >  2:303934/1‑90 (MQ=255)
                                                                     ttttGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGc              >  1:376988/1‑90 (MQ=255)
                                                                     ttttGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGc              >  1:1206158/1‑90 (MQ=255)
                                                                                 gTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGcc                               <  2:454371/61‑1 (MQ=255)
                                                                                 gTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGcc                               >  1:454371/1‑61 (MQ=255)
                                                                                 gTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAg  >  2:1148550/1‑90 (MQ=255)
                                                                                 gTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAg  >  2:637286/1‑90 (MQ=255)
                                                                                         |                                                                                  
GGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAG  >  NZ_CP009273/1060894‑1061065

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
GGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAGCTGTGGCAGGGATAGTTT  >  NZ_CP009273/1060894‑1061083
                                                                                              |                                                                                               
GGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATA                                                                                           >  SRR3722089.271214/1‑100 (MQ=60)
                      GCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGA                                                                     <  SRR3722089.1220218/100‑1 (MQ=60)
                      GCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGA                                                                     <  SRR3722089.817504/100‑1 (MQ=60)
                          GTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGA                                                                 >  SRR3722089.762293/1‑100 (MQ=60)
                           TTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGCTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGAT                                                                <  SRR3722089.1246389/100‑1 (MQ=60)
                                 AAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAA                                                          >  SRR3722089.8852/1‑100 (MQ=60)
                                   GATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAG                                                        <  SRR3722089.1232207/100‑1 (MQ=60)
                                   GATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAG                                                        <  SRR3722089.761900/100‑1 (MQ=60)
                                     TGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAG                                                      >  SRR3722089.957088/1‑100 (MQ=60)
                                      GTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCTAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGG                                                     <  SRR3722089.1235162/100‑1 (MQ=60)
                                      GTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTTTGTCACACTCTTTAGTGATTGATAACAAAAGAGG                                                     <  SRR3722089.773019/100‑1 (MQ=60)
                                      GTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGG                                                     <  SRR3722089.993730/100‑1 (MQ=60)
                                      GTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGG                                                     <  SRR3722089.438899/100‑1 (MQ=60)
                                                    AGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAAC                                       <  SRR3722089.715856/100‑1 (MQ=60)
                                                           CAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGC                                >  SRR3722089.1224499/1‑100 (MQ=60)
                                                           CAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGC                                >  SRR3722089.381596/1‑100 (MQ=60)
                                                                       TTGCGAACCTGTTCCCGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACctgtctcttatacacatct                    >  SRR3722089.459812/1‑81 (MQ=60)
                                                                                      CGGAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAGCTGTGGCAGGGATAG     <  SRR3722089.446136/100‑1 (MQ=60)
                                                                                          AAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAGCTGTGGCAGGGATAGTTT  <  SRR3722089.598152/100‑1 (MQ=60)
                                                                                              |                                                                                               
GGCATAGTTTGCGTCAAACCTTGCCTGTTTTTGAAGATGTATATAGAAAAACAGGCGTTCAACAAGCCATTTTGCGAACCTGTTCCCGGAAAAAAGTCATATTTCTGTCACACTCTTTAGTGATTGATAACAAAAGAGGTGCCAGGAATGAACAAAACGCTAATCGCCGCAGCTGTGGCAGGGATAGTTT  >  NZ_CP009273/1060894‑1061083

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