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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,064,259 T→G V93G (GTG→GGG)  BW25113_RS10440 → GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,064,2590TG95.5% 145.4 / 0.0 44coding (278/711 nt)BW25113_RS10440GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (22/20);  minor base . (1/1);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAG  >  NZ_CP009273/2064172‑2064341
                                                                                       |                                                                                  
tACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGc                                                                                  >  2:299047/1‑90 (MQ=255)
    gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                              >  1:18593/1‑90 (MQ=255)
    gtgGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGtt                                                                              >  1:800179/1‑90 (MQ=255)
      ggCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTtg                                                                            >  1:272701/1‑90 (MQ=255)
        cTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTtgtg                                                                          >  2:474039/1‑90 (MQ=255)
          gATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCGGGGGCTGTTTgtggt                                                                        >  1:320255/1‑90 (MQ=255)
           aTTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTgtggtg                                                                       >  2:899317/1‑90 (MQ=255)
              aaaGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACg                                                                    >  1:174816/1‑90 (MQ=255)
                    gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                              <  2:546035/90‑1 (MQ=255)
                    gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                              >  2:25022/1‑90 (MQ=255)
                    gATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGcc                                                              >  1:722003/1‑90 (MQ=255)
                                ctgtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAg                                                  <  1:792213/90‑1 (MQ=255)
                                  gtctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGcc                                                >  1:393642/1‑90 (MQ=255)
                                   tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                               <  2:393642/90‑1 (MQ=255)
                                   tctgtGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGccc                                               <  1:25022/90‑1 (MQ=255)
                                      gtgGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGc                                            >  1:689263/1‑90 (MQ=255)
                                      gtgGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGGGGTGGCGCGGGCCGACAAAACGGAGCCCTGc                                            >  2:770040/1‑90 (MQ=255)
                                        ggATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCa                                          <  1:116313/90‑1 (MQ=255)
                                         gATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCAt                                         <  1:900337/90‑1 (MQ=255)
                                          aTGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATg                                        >  2:729506/1‑90 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCTGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     >  2:474701/1‑90 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                      >  1:521331/1‑90 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                      >  2:386807/1‑90 (MQ=255)
                                            gAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGaa                                      >  1:324769/1‑90 (MQ=255)
                                             tgtcttttTGCCGACCCACCCTGCAGTGCGGACACCAGCAGGGGCTGTTTGTGGTCACGCAGGCCGCCAAAACGGAGCCCTGCCATGAAt                                     <  2:704917/86‑1 (MQ=255)
                                             aGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     <  2:429363/90‑1 (MQ=255)
                                             aGTATTTCTGGCGAACCATCAGGCAGTGAGGCCATCCGCAGGGGTTGTTTGTGGTGACCCAGGCCGACAAACCGGAGCCCTGCCATGAAt                                     <  2:113389/90‑1 (MQ=255)
                                              gTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATg                                    >  1:743353/1‑90 (MQ=255)
                                                        cGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGcc                          <  2:722003/90‑1 (MQ=255)
                                                           cacaTCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     <  2:687420/76‑1 (MQ=255)
                                                           cacaTCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAt                                     >  1:687420/1‑76 (MQ=255)
                                                                   gCAGTGCGGACATCAGCTGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:474701/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  2:800179/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:885886/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  2:521331/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:907203/90‑1 (MQ=255)
                                                                    cAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  1:729506/90‑1 (MQ=255)
                                                                       tGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGGCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTtct           >  1:541037/1‑90 (MQ=255)
                                                                           gACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTg       >  1:75877/1‑90 (MQ=255)
                                                                              aTCAGCGGGGGCTGTTTGTGGGGCCGCAGGCCGACAAAACGGGGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcac    <  1:440292/90‑1 (MQ=255)
                                                                               tCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcaca   <  2:781671/90‑1 (MQ=255)
                                                                                cagcagGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcacag  >  1:213163/1‑90 (MQ=255)
                                                                                cagcagGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTgcacag  <  1:116479/90‑1 (MQ=255)
                                                                                    agGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCt              <  2:787254/74‑1 (MQ=255)
                                                                                       |                                                                                  
TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAG  >  NZ_CP009273/2064172‑2064341

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATA  >  NZ_CP009273/2064161‑2064351
                                                                                                  |                                                                                            
ACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGG                                                                                             >  SRR3722072.748629/1‑100 (MQ=60)
     ACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTT                                                                                        >  SRR3722072.18943/1‑100 (MQ=60)
     ACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTT                                                                                        >  SRR3722072.814160/1‑100 (MQ=60)
       CTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTG                                                                                      >  SRR3722072.277363/1‑100 (MQ=60)
           TACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCGGGGGCTGTTTGTGGT                                                                                  >  SRR3722072.325852/1‑100 (MQ=60)
               GTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACG                                                                              >  SRR3722072.177838/1‑100 (MQ=60)
                     GATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCC                                                                        >  SRR3722072.734843/1‑100 (MQ=60)
                                   ACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCC                                                          >  SRR3722072.400724/1‑100 (MQ=60)
                                       TGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGC                                                      >  SRR3722072.701593/1‑100 (MQ=60)
                                           CTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATG                                                  <  SRR3722072.806077/100‑1 (MQ=60)
                                             GTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAA                                                >  SRR3722072.330454/1‑100 (MQ=60)
                                             GTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAA                                                >  SRR3722072.530711/1‑100 (MQ=60)
                                              TCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAAT                                               <  SRR3722072.25462/100‑1 (MQ=60)
                                               CTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATG                                              >  SRR3722072.756508/1‑100 (MQ=60)
                                                   GGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGAT                                          <  SRR3722072.118454/100‑1 (MQ=60)
                                                    GATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATA                                         <  SRR3722072.915831/100‑1 (MQ=60)
                                                            TTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCTGT                                 >  SRR3722072.699720/1‑100 (MQ=60)
                                                                        CATCCTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGGCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCT                     >  SRR3722072.550781/1‑100 (MQ=60)
                                                                            CTGCAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTG                 >  SRR3722072.77287/1‑100 (MQ=60)
                                                                              GCAGTGCGGACATCAGCTGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCAC              <  SRR3722072.483353/100‑1 (MQ=60)
                                                                               CAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCAC              <  SRR3722072.901127/100‑1 (MQ=60)
                                                                               CAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCAC              <  SRR3722072.922813/100‑1 (MQ=60)
                                                                               CAGTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCAC              <  SRR3722072.742450/100‑1 (MQ=60)
                                                                                 GTGCGGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAG            >  SRR3722072.216782/1‑100 (MQ=60)
                                                                                     GGACATCAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACctgtct        >  SRR3722072.801047/1‑94 (MQ=60)
                                                                                         ATCAGCGGGGGCTGTTTGTGGGGCCGCAGGCCGACAAAACGGGGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAA    <  SRR3722072.448338/100‑1 (MQ=60)
                                                                                           CAGCAGGGGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATA  <  SRR3722072.118623/100‑1 (MQ=60)
                                                                                                  |                                                                                            
ACTTGACCTGGTACTGTGGCTGATTAAAGCCGATGACCGTGCCCTGTCTGTGGATGAGTATTTCTGGCGACACATCCTGCAGTGCGGACATCAGCAGGTGCTGTTTGTGGTGACGCAGGCCGACAAAACGGAGCCCTGCCATGAATGGGATATGGCCGGTATTCAGCCTTCTCCTGCACAGGCACAGAATA  >  NZ_CP009273/2064161‑2064351

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