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 F4 I96 R1 227 23.3 882352 98.0% 864704 126.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,568,721 A→G I145T (ATC→ACC)  mcrB ← 5‑methylcytosine‑specific restriction endonuclease subunit McrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,568,7210AG100.0% 89.5 / NA 25I145T (ATC→ACC) mcrB5‑methylcytosine‑specific restriction endonuclease subunit McrB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/11);  total (14/11)

TTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTT  >  NZ_CP009273/4568588‑4568856
                                                                                                                                     |                                                                                                                                       
tttGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTAtg                                                                                                                                    <  2:308003/139‑1 (MQ=255)
tttGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTAtg                                                                                                                                    <  2:221628/139‑1 (MQ=255)
                         tCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCggg                                                                                                           >  2:396533/1‑139 (MQ=255)
                                                        ccAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTggg                                                                             >  1:313905/1‑138 (MQ=255)
                                                         cAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGaa                                                                           <  2:201760/139‑1 (MQ=255)
                                                                tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAAc                                                                          >  2:89205/1‑133 (MQ=255)
                                                                tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAAc                                                                          <  1:89205/133‑1 (MQ=255)
                                                                tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg                                                                     >  2:147346/1‑138 (MQ=255)
                                                                tATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCg                                                                     >  2:324711/1‑138 (MQ=255)
                                                                          tttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                                                                       >  2:203395/1‑126 (MQ=255)
                                                                          tttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAgg                                                                       <  1:203395/126‑1 (MQ=255)
                                                                          tttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTg                                                          >  2:100127/1‑139 (MQ=255)
                                                                          tttttGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTg                                                          >  1:286754/1‑139 (MQ=255)
                                                                                     ggTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTAtttt                                                <  2:150457/138‑1 (MQ=255)
                                                                                       tGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTa                                             <  1:324711/139‑1 (MQ=255)
                                                                                              acacTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGatat                                       <  1:322370/138‑1 (MQ=255)
                                                                                                                               tCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAAc                                                                          <  2:205110/70‑1 (MQ=255)
                                                                                                                               tCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAAc                                                                          >  1:205110/1‑70 (MQ=255)
                                                                                                                                 gTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGtt   >  2:24246/1‑139 (MQ=255)
                                                                                                                                 gTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGtt   >  2:390664/1‑139 (MQ=255)
                                                                                                                                 gTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACATTTTGGGAACAgg                                                                       <  1:400295/71‑1 (MQ=255)
                                                                                                                                 gTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACATTTTGGGAACAgg                                                                       >  2:400295/1‑71 (MQ=255)
                                                                                                                                  ttGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGa        <  2:124399/133‑1 (MQ=255)
                                                                                                                                  ttGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGa        >  1:124399/1‑133 (MQ=255)
                                                                                                                                  ttGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGttt  >  1:428889/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                       
TTTGAGTATCGTCTCTATTGTGGTTTCAGGGATAAACAAATCATTTAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTT  >  NZ_CP009273/4568588‑4568856

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATG  >  NZ_CP009273/4568633‑4568864
                                                                                        |                                                                                                                                               
TAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGA                                                                                    >  SRR3722119.317099/1‑150 (MQ=60)
            CAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAG                                                                       <  SRR3722119.90082/151‑1 (MQ=60)
               cagGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGAC                                                                    <  SRR3722119.205572/148‑1 (MQ=60)
               agacgtcatacgtgattcgccttagtctcgtgggctcggagatgtgtataagagacagAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACATTTTGGGAACAGGCGTAATAGGTC                                                                    <  SRR3722119.404292/93‑1 (MQ=60)
                  GTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGT                                                                 >  SRR3722119.289711/1‑151 (MQ=60)
                                         GTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACC                                          <  SRR3722119.328025/151‑1 (MQ=60)
                                                AACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGT                                    <  SRR3722119.325662/150‑1 (MQ=60)
                                                                       TTAATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGctgtctcttatacacatctccgagcccacgagactaaggcgaatgtcgtatgccgtctt            >  SRR3722119.207307/1‑92 (MQ=60)
                                                                          ATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTT         >  SRR3722119.125638/1‑151 (MQ=60)
                                                                          ATTTATAGTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTT         >  SRR3722119.433196/1‑151 (MQ=60)
                                                                                 GTCGTTGGTTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATG  >  SRR3722119.233676/1‑151 (MQ=60)
                                                                                        |                                                                                                                                               
TAACGCATCTTCCAGACAGTATGATTCAGTTTTTGACATAGGTGGAATAACACTCTTGCCAGAATTAAATATTAATTTATAGTCGTTGATTATGTTGTCCAGCATAGAGGCAAATCGGGTGTAATCAATACCCTGTGAGACTTTTTGGGAACAGGCGTAATAGGACTGTCCGTATTTTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATG  >  NZ_CP009273/4568633‑4568864

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