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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,972,691 A→G intergenic (‑237/‑543) flhD ← / → uspC flagellar transcriptional regulator FlhD/universal stress protein UspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,972,6910AG96.4% 80.0 / ‑3.5 28intergenic (‑237/‑543)flhD/uspCflagellar transcriptional regulator FlhD/universal stress protein UspC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (15/12);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTG  >  NZ_CP009273/1972605‑1972776
                                                                                      |                                                                                     
tACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgggggg                                                                                    >  1:48589/1‑90 (MQ=255)
         aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGa                                                                           <  2:430968/90‑1 (MQ=255)
                  tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCaata                                                                  >  2:90046/1‑90 (MQ=255)
                        cTTATCTAAGATTTTTCCCAAATCGACGCAACGGTACTCGTCACAACACGCACATACAACGGGGGGGGGCTGCGATTTCCAATAATGCGt                                                            <  2:257680/90‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATg                                                               <  2:384524/86‑1 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATg                                                               >  1:384524/1‑86 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                           >  1:219866/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                           >  2:456985/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                           >  2:110946/1‑90 (MQ=255)
                         ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTg                                                           >  1:380961/1‑90 (MQ=255)
                          tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGa                                                          >  1:105932/1‑90 (MQ=255)
                             cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGc                                                       <  2:400524/90‑1 (MQ=255)
                                      ttCCTAAATCGACGCAACTGTCCTCGTCACTACACGCACATACAACGGGGGGTGACCGGGGTTTTCATTAATGCGTGAGCCAGCTcacac                                              >  1:396286/1‑90 (MQ=255)
                                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACaaaa                                          <  2:380961/90‑1 (MQ=255)
                                                gACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCa                                    >  1:91402/1‑90 (MQ=255)
                                                 aCGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCaa                                   <  2:331218/90‑1 (MQ=255)
                                                  cGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAt                                  <  1:374907/90‑1 (MQ=255)
                                                   gCAACTGTACTCGTCACTACCCGCAAACACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAtt                                 <  1:456985/90‑1 (MQ=255)
                                                       cTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACtt                             <  2:457490/90‑1 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                      >  1:350041‑M1/1‑85 (MQ=255)
                                                              cGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                      >  1:230558‑M1/1‑85 (MQ=255)
                                                                 cACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtg                   >  2:385831‑M1/1‑82 (MQ=255)
                                                                   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                 >  1:348971‑M1/1‑80 (MQ=255)
                                                                      cacGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac              <  2:219866‑M1/90‑14 (MQ=255)
                                                                        cGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa            >  1:251402‑M1/1‑75 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc       <  1:296270‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc       <  1:385831‑M1/90‑21 (MQ=255)
                                                                             aTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggggcgaacaagtccc       >  2:526811‑M1/1‑70 (MQ=255)
                                                                                  aCGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata  >  2:119340‑M1/1‑65 (MQ=255)
                                                                                      |                                                                                     
TACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTG  >  NZ_CP009273/1972605‑1972776

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 29 ≤ 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
TAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTG  >  NZ_CP009273/1972594‑1972781
                                                                                                 |                                                                                          
TAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGG                                                                                          <  SRR3722109.355577/100‑1 (MQ=60)
 AATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGG                                                                                         >  SRR3722109.49164/1‑100 (MQ=60)
   cgtcagatgtgtataagagacaGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGG                                                                                       <  SRR3722109.546151/78‑1 (MQ=60)
                          GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTG                                                                >  SRR3722109.222501/1‑100 (MQ=60)
                          GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTG                                                                >  SRR3722109.390007/1‑100 (MQ=60)
                          GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTG                                                                >  SRR3722109.386375/1‑100 (MQ=60)
                           TCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGA                                                               >  SRR3722109.107259/1‑100 (MQ=60)
                                            attgacattcgtcggcagcgtcagatgtgtataagagacaGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAAC                                              <  SRR3722109.451929/60‑1 (MQ=60)
                                                 TTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCA                                         >  SRR3722109.92537/1‑100 (MQ=60)
                                                             CGCAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgga                             <  SRR3722109.380209/100‑4 (MQ=60)
                                                              GCAACTGTACTCGTCACTACCCGCAAACACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaa                            <  SRR3722109.463835/100‑5 (MQ=60)
                                                               CAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                           >  SRR3722109.354781/1‑95 (MQ=60)
                                                               CAACTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                           >  SRR3722109.233291/1‑95 (MQ=60)
                                                                    GTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                      >  SRR3722109.353678/1‑90 (MQ=60)
                                                                         CGTCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                 >  SRR3722109.254372/1‑85 (MQ=60)
                                                                                     CACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg     >  SRR3722109.501387/1‑73 (MQ=60)
                                                                                     CACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg     >  SRR3722109.505254/1‑73 (MQ=60)
                                                                                     CACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg     >  SRR3722109.546835/1‑73 (MQ=60)
                                                                                        ATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga  <  SRR3722109.300093/100‑31 (MQ=60)
                                                                                        ATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga  <  SRR3722109.391343/100‑31 (MQ=60)
                                                                                                 |                                                                                          
TAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTG  >  NZ_CP009273/1972594‑1972781

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