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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,215,724 T→G W657G (TGG→GGG)  metH → methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,215,7240TG93.0% 132.3 / 1.8 43W657G (TGG→GGG) metHmethionine synthase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (13/27);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAG  >  NZ_CP009273/4215639‑4215810
                                                                                     |                                                                                      
cGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGa                                                                                    <  1:141016/90‑1 (MQ=255)
  aGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  <  1:329643/90‑1 (MQ=255)
  aGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  <  1:516689/90‑1 (MQ=255)
  aGCGGTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  <  2:281340/90‑1 (MQ=255)
     gTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGa                                                                               >  1:176131/1‑90 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATaa                                                                           <  2:615182/90‑1 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  2:79114/90‑1 (MQ=255)
         aCTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATaa                                                                           <  2:102939/90‑1 (MQ=255)
            ggAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  <  1:287455/80‑1 (MQ=255)
            ggAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAg                                                                                  >  2:287455/1‑80 (MQ=255)
              aGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTc                                                                      <  1:243461/90‑1 (MQ=255)
              aGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTc                                                                      <  1:325371/90‑1 (MQ=255)
               gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGGGGGGCGCTCGGGGGAAGGGAATAAACGTCt                                                                     >  1:667063/1‑90 (MQ=255)
               gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCt                                                                     >  1:695801/1‑90 (MQ=255)
               gCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGCGGCGCTCGGGGGAAGGGAATAAACGTCt                                                                     >  1:576165/1‑90 (MQ=255)
                  tGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGa                                                                  >  1:614664/1‑90 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  2:181189/90‑1 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  1:499597/90‑1 (MQ=255)
                          aaTATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGc                                                          <  1:709906/90‑1 (MQ=255)
                           atatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCt                                                         <  1:432789/90‑1 (MQ=255)
                           atatCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCt                                                         <  1:270270/90‑1 (MQ=255)
                                gcgGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCa                                                    <  2:614053/90‑1 (MQ=255)
                                  ggCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa                                                  <  2:526715/90‑1 (MQ=255)
                                    cagcaAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAgg                                                <  1:655456/90‑1 (MQ=255)
                                           aCCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTAcc                                         <  2:89995/90‑1 (MQ=255)
                                                     ccGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCg                               >  2:329200/1‑90 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  1:130085/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  1:452884/90‑1 (MQ=255)
                                                        ccAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGc                            <  2:122343/90‑1 (MQ=255)
                                                          aaCGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAg                          <  1:495189/90‑1 (MQ=255)
                                                           aCGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAgg                         <  1:329200/90‑1 (MQ=255)
                                                             gCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGAt                       <  1:421945/90‑1 (MQ=255)
                                                               ccAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATAc                     >  1:426615/1‑90 (MQ=255)
                                                                     ggCGGAGTGGCACTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCgaaga               >  1:316165/1‑90 (MQ=255)
                                                                       cGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAg                          >  2:278035/1‑77 (MQ=255)
                                                                       cGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAg                          <  1:278035/77‑1 (MQ=255)
                                                                        ggAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGc            >  1:14700/1‑90 (MQ=255)
                                                                           gTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCg         >  2:709250/1‑90 (MQ=255)
                                                                               cgcTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcag     >  1:593419/1‑90 (MQ=255)
                                                                               cgcTCGTGGGAAGTGAATAAACGTCTGAATACTCGCTGGGCAAAGGCATTACCGGGGTTATCGAGCAGGATACCGAAGAAGCCCGCcagc    >  1:593420/1‑90 (MQ=255)
                                                                               cgcTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCg                               >  2:173323/1‑64 (MQ=255)
                                                                               cgcTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCg                               <  1:173323/64‑1 (MQ=255)
                                                                                  tCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCcagcag  <  1:709250/90‑1 (MQ=255)
                                                                                     |                                                                                      
CGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAG  >  NZ_CP009273/4215639‑4215810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
TGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCC  >  NZ_CP009273/4215633‑4215820
                                                                                           |                                                                                                
TGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGAGGAGTGGCGCTCGGGGGAAGTG                                                                                          <  SRR3722112.285554/100‑1 (MQ=60)
 GGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGA                                                                                         >  SRR3722112.178293/1‑100 (MQ=60)
      CGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAA                                                                                    <  SRR3722112.142781/100‑1 (MQ=60)
        AGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACG                                                                                  <  SRR3722112.290790/100‑1 (MQ=60)
        AGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACG                                                                                  <  SRR3722112.524264/100‑1 (MQ=60)
        AGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACG                                                                                  <  SRR3722112.333597/100‑1 (MQ=60)
           GTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGGGGGGCGCTCGGGGGAAGGGAATAAACGTCT                                                                               >  SRR3722112.677471/1‑100 (MQ=60)
           GTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCT                                                                               >  SRR3722112.706667/1‑100 (MQ=60)
           GTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGCGGCGCTCGGGGGAAGGGAATAAACGTCT                                                                               >  SRR3722112.584959/1‑100 (MQ=60)
              TACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGA                                                                            >  SRR3722112.624117/1‑100 (MQ=60)
                    AGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTC                                                                      <  SRR3722112.246347/100‑1 (MQ=60)
                    AGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTC                                                                      <  SRR3722112.329233/100‑1 (MQ=60)
                                AATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGG                                                          <  SRR3722112.721009/100‑1 (MQ=60)
                                AATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGG                                                          <  SRR3722112.506836/100‑1 (MQ=60)
                                 ATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGC                                                         <  SRR3722112.438689/100‑1 (MQ=60)
                                 ATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGC                                                         <  SRR3722112.273432/100‑1 (MQ=60)
                                          CAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAG                                                <  SRR3722112.665681/100‑1 (MQ=60)
                                                           tgtgtataagagacagGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATAC                               <  SRR3722112.175448/84‑1 (MQ=60)
                                                           CCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATAC                               >  SRR3722112.432340/1‑100 (MQ=60)
                                                              CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                            <  SRR3722112.459246/100‑1 (MQ=60)
                                                              CCAACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGA                            <  SRR3722112.131699/100‑1 (MQ=60)
                                                                cagGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAG                          <  SRR3722112.281274/97‑1 (MQ=60)
                                                                AACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAG                          <  SRR3722112.502349/100‑1 (MQ=60)
                                                                 ACGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGA                         <  SRR3722112.333143/100‑1 (MQ=60)
                                                                 ACGCCCAGCAGGCGGAGTGGCACTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGA                         >  SRR3722112.319822/1‑100 (MQ=60)
                                                                   GCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAG                       <  SRR3722112.427558/100‑1 (MQ=60)
                                                                    CCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGC                      >  SRR3722112.14889/1‑100 (MQ=60)
                                                                           GGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAG               >  SRR3722112.602488/1‑100 (MQ=60)
                                                                           GGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGAATACTCGCTGGGCAAAGGCATTACCGGGGTTATCGAGCAGGATACCGAAGAAGCCCGCCAGC              >  SRR3722112.602489/1‑100 (MQ=60)
                                                                                        TCGGGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTCCCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCC  <  SRR3722112.720338/100‑1 (MQ=60)
                                                                                           |                                                                                                
TGGCACCGAGCGTTTACTGGAGCTTGCCGAGAAATATCGCGGCAGCAAAACCGACGACACCGCCAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAGGCATTACCGAGTTTATCGAGCAGGATACCGAAGAAGCCCGCCAGCAGGCTACGCGCC  >  NZ_CP009273/4215633‑4215820

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