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 886,189 A→G G7G (GGT→GGC grxA ← glutaredoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273886,1890AG100.0% 144.6 / NA 42G7G (GGT→GGCgrxAglutaredoxin 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/15);  total (27/15)

ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGT  >  NZ_CP009273/886103‑886270
                                                                                      |                                                                                 
aCATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                <  1:981647/90‑1 (MQ=255)
aCATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCa                                                                                <  1:263310/90‑1 (MQ=255)
   tactgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCaaaa                                                                             >  2:769358/1‑90 (MQ=255)
     ctgatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAAt                                                                           >  2:809805/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                         >  1:805811/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                         >  2:124564/1‑90 (MQ=255)
       gatactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATaa                                                                         >  2:483216/1‑90 (MQ=255)
         tactgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACg                                                                       >  1:616424/1‑90 (MQ=255)
          actgaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACgg                                                                      >  2:1083772/1‑90 (MQ=255)
             gaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                                   >  1:134772/1‑90 (MQ=255)
             gaAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGttt                                                                   <  2:781859/90‑1 (MQ=255)
               aaaTCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGc                                                                 <  1:157304/90‑1 (MQ=255)
                 atcatcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCat                                                               <  2:790868/90‑1 (MQ=255)
                   catcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCatta                                                             >  1:562227/1‑90 (MQ=255)
                     tcGCGTTCATTGCTCGATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCattatt                                                           >  1:956437/1‑90 (MQ=255)
                         gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTtctc                                                       >  2:877449/1‑90 (MQ=255)
                                gCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCAt                                                <  2:968679/90‑1 (MQ=255)
                                 cTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                                               <  1:690292/90‑1 (MQ=255)
                                 cTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATa                                               <  1:1165508/90‑1 (MQ=255)
                                     aTTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGAtt                                           >  2:313114/1‑90 (MQ=255)
                                      tttCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCAAAGAttt                                          <  1:1199620/90‑1 (MQ=255)
                                          tcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGc                                      >  2:757718/1‑90 (MQ=255)
                                          tcAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGc                                      >  1:379137/1‑90 (MQ=255)
                                           cAGCCAGATCTTCTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                     <  2:462723/90‑1 (MQ=255)
                                                         gcacgcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGc                       <  2:805811/90‑1 (MQ=255)
                                                           acgcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTa                     >  2:1044537/1‑90 (MQ=255)
                                                             gcacACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTaaa                   >  1:1032869/1‑90 (MQ=255)
                                                                      aaGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCg          >  2:305656/1‑90 (MQ=255)
                                                                        gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTAt                                        <  1:1185130/58‑1 (MQ=255)
                                                                        gggCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTAt                                        >  2:1185130/1‑58 (MQ=255)
                                                                          gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                     <  2:583881/59‑1 (MQ=255)
                                                                          gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGcc                                     >  1:583881/1‑59 (MQ=255)
                                                                          gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGa      >  1:693998/1‑90 (MQ=255)
                                                                          gCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGa      >  1:1051348/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCAt                                                <  2:420700/46‑1 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCAt                                                >  1:420700/1‑46 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa    >  2:674236/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa    >  2:315849/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa    >  1:525906/1‑90 (MQ=255)
                                                                            aaCCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGaaa    >  1:69726/1‑90 (MQ=255)
                                                                              cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt  >  1:319112/1‑90 (MQ=255)
                                                                              cccGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGt  <  1:562082/90‑1 (MQ=255)
                                                                                      |                                                                                 
ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGT  >  NZ_CP009273/886103‑886270

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ 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
TCTACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACC  >  NZ_CP009273/886100‑886280
                                                                                         |                                                                                           
TCTACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAA                                                                                   >  SRR3722089.817022/1‑100 (MQ=60)
  TACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACG                                                                                 >  SRR3722089.624137/1‑100 (MQ=60)
   ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGG                                                                                <  SRR3722089.266683/100‑1 (MQ=60)
   ACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGG                                                                                <  SRR3722089.996193/100‑1 (MQ=60)
      TACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTT                                                                             >  SRR3722089.136410/1‑100 (MQ=60)
            TACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTA                                                                       >  SRR3722089.568807/1‑100 (MQ=60)
              CTGAAAATCATCGCGTTCATTGCTCGATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATT                                                                     >  SRR3722089.970502/1‑100 (MQ=60)
                  AAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTC                                                                 <  SRR3722089.159253/100‑1 (MQ=60)
                                   GCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGC                                                >  SRR3722089.383767/1‑100 (MQ=60)
                                    CTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCC                                               <  SRR3722089.699292/100‑1 (MQ=60)
                                    CTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCC                                               <  SRR3722089.1183227/100‑1 (MQ=60)
                                         TTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCAAAGATTTATGCCTGTAA                                          <  SRR3722089.1217859/100‑1 (MQ=60)
                                           gtcagatgtgtataagagacaGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATG                                        <  SRR3722089.1203170/79‑1 (MQ=60)
                                                      TCTTTTGCACGCACACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAA                             >  SRR3722089.1048299/1‑100 (MQ=60)
                                                                   CACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCctgtctcttatacacatctga                >  SRR3722089.590889/1‑79 (MQ=60)
                                                                   CACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGA                >  SRR3722089.703042/1‑100 (MQ=60)
                                                                   CACAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGA                >  SRR3722089.1067108/1‑100 (MQ=60)
                                                                     CAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTctcttatacacatctgacgctgccgacgac              >  SRR3722089.425779/1‑70 (MQ=60)
                                                                     CAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAA              >  SRR3722089.70572/1‑100 (MQ=60)
                                                                     CAGTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAA              >  SRR3722089.532047/1‑100 (MQ=60)
                                                                       GTAAGGGCAACCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGT            >  SRR3722089.323152/1‑100 (MQ=60)
                                                                                 CCCGAACGGCCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACC  <  SRR3722089.568658/100‑1 (MQ=60)
                                                                                         |                                                                                           
TCTACATACTGATACTGAAAATCATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACC  >  NZ_CP009273/886100‑886280

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