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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,965,782 T→G E389D (GAA→GAC tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,965,7820TG97.6% 129.8 / ‑3.6 41E389D (GAA→GACtarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (17/23);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAA  >  NZ_CP009273/1965696‑1965864
                                                                                      |                                                                                  
aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                                 >  1:900824/1‑90 (MQ=255)
aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                                 >  1:719863/1‑90 (MQ=255)
 gggCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCAc                                                                                >  2:28771/1‑90 (MQ=255)
           ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                      >  1:1272584/1‑90 (MQ=255)
           ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                      <  2:1029551/90‑1 (MQ=255)
             ctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCgg                                                                    >  1:216379/1‑90 (MQ=255)
                     gccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAcc                                                            >  2:207375/1‑90 (MQ=255)
                      ccTGGGCGCTGCGCCTGGCAAGTTTACGCACTCCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACcg                                                           <  2:714324/90‑1 (MQ=21)
                       cTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACcgc                                                          >  1:707704/1‑90 (MQ=255)
                       cTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACcgc                                                          >  2:275440/1‑90 (MQ=255)
                              cTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTc                                                   >  1:214981/1‑90 (MQ=255)
                                gCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAg                                                 >  2:902208/1‑90 (MQ=255)
                                    cTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                             <  2:719863/90‑1 (MQ=255)
                                    cTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                             <  1:645522/90‑1 (MQ=255)
                                    cTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                             <  2:900824/90‑1 (MQ=255)
                                     tGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGa                                            <  1:28771/90‑1 (MQ=255)
                                           gATTACGCCCTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGatat                                      <  1:50132/90‑1 (MQ=37)
                                           gATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGatat                                      <  2:214981/90‑1 (MQ=255)
                                              tACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAg                                   >  1:1058903/1‑90 (MQ=255)
                                               aCGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGt                                  >  1:380630/1‑90 (MQ=255)
                                                  cACTTCCCCCGCCACCACGGCAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                              <  2:719911/90‑1 (MQ=255)
                                                   ccTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                              <  2:654883/89‑1 (MQ=255)
                                                   ccTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                              <  2:163162/89‑1 (MQ=255)
                                                   aCTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                              <  2:715758/90‑1 (MQ=255)
                                                   aCTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                              <  2:216379/90‑1 (MQ=255)
                                                   aCTCCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                              <  1:694395/90‑1 (MQ=255)
                                                    cTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGa                             >  2:821362/1‑90 (MQ=255)
                                                     ttCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGaa                            <  1:1547235/90‑1 (MQ=255)
                                                          ccGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCaa                       <  2:421613/90‑1 (MQ=255)
                                                           cGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAAt                      <  2:82843/90‑1 (MQ=255)
                                                            gCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATa                     >  2:367219/1‑90 (MQ=255)
                                                              caccacGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                      <  2:298506/39‑1 (MQ=25)
                                                              caccacGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                      >  1:298506/1‑39 (MQ=37)
                                                                     gCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa            <  2:806183/90‑1 (MQ=255)
                                                                     gCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa            <  1:703886/90‑1 (MQ=255)
                                                                     gCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGAAAGGCAATACCGTCGATa            <  1:1106605/90‑1 (MQ=255)
                                                                      cAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATaa           <  1:625282/90‑1 (MQ=255)
                                                                       aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc          <  1:321810/90‑1 (MQ=255)
                                                                             aCGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGAt    >  1:1140054/1‑90 (MQ=255)
                                                                               ggCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATaa  >  2:339943/1‑90 (MQ=255)
                                                                               ggCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATaa  >  2:1502211/1‑90 (MQ=255)
                                                                                      |                                                                                  
AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAA  >  NZ_CP009273/1965696‑1965864

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ 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
GTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATT  >  NZ_CP009273/1965686‑1965875
                                                                                                |                                                                                             
GTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCA                                                                                            >  SRR3722078.912918/1‑100 (MQ=60)
GTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCA                                                                                            >  SRR3722078.728768/1‑100 (MQ=60)
  CTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGTCCCTTGTCACC                                                                                          >  SRR3722078.483661/1‑100 (MQ=60)
           GGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGC                                                                                 >  SRR3722078.1292780/1‑100 (MQ=60)
             GCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGG                                                                               >  SRR3722078.219005/1‑100 (MQ=60)
                       CTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGC                                                                     >  SRR3722078.716289/1‑100 (MQ=60)
                              CGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTC                                                              >  SRR3722078.217582/1‑100 (MQ=60)
                                              CTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAG                                              <  SRR3722078.652943/100‑1 (MQ=60)
                                              CTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAG                                              >  SRR3722078.1074828/1‑100 (MQ=60)
                                               TGGCAAGATTCCGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGT                                             <  SRR3722078.29128/100‑1 (MQ=60)
                                               TGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGT                                             >  SRR3722078.385279/1‑100 (MQ=60)
                                                     GATTACGCCCTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAA                                       <  SRR3722078.50713/100‑1 (MQ=60)
                                                             cctcCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATAC                               <  SRR3722078.702669/96‑1 (MQ=60)
                                                              CTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCctgtctcttatacacatctgacgctgccgacgatggcggtc                              >  SRR3722078.302157/1‑59 (MQ=60)
                                                               TTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCG                             <  SRR3722078.1572442/100‑1 (MQ=60)
                                                                             CGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGAT               >  SRR3722078.1157569/1‑100 (MQ=60)
                                                                               GCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAA             <  SRR3722078.712387/100‑1 (MQ=60)
                                                                               GCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGAAAGGCAATACCGTCGATAACGCTGATAA             <  SRR3722078.1123491/100‑1 (MQ=60)
                                                                                CAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAAT            <  SRR3722078.632478/100‑1 (MQ=60)
                                                                                 AAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATG           <  SRR3722078.325744/100‑1 (MQ=60)
                                                                                          GCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATT  >  SRR3722078.844871/1‑100 (MQ=60)
                                                                                                |                                                                                             
GTCTTCAATGAGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGATAATGTCGGCAATT  >  NZ_CP009273/1965686‑1965875

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