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
A4 F13 I0 R1 140 52.5 2164708 98.8% 2138731 146.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 4,001,773 T→G 36.8% intergenic (+90/‑44) dmoA‑II → / → allE dimethyl‑sulfide monooxygenase/S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,001,7730TG36.8% 19.9 / 5.3 24intergenic (+90/‑44)dmoA‑II/allEdimethyl‑sulfide monooxygenase/S‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (4/11);  new base G (5/4);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGC  >  NC_002947_CJ‑RC/4001628‑4001908
                                                                                                                                                  |                                                                                                                                       
gCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGt                                                                                                                                         <  5:95598/146‑1 (MQ=255)
       cggcCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCa                                                                                                                                 <  4:116051/147‑1 (MQ=255)
            aGCTGCCCGAACGGCACGTTGGACGGCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGGGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCa                                                                                                                            <  6:218961/147‑1 (MQ=255)
              cTGCCCGAACGGCACGTTGGACTTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCTCTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATg                                                                                                                          <  6:141180/147‑1 (MQ=255)
              cTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCCCCGGCCCGGGGGGTTGAAGGCATg                                                                                                                          >  6:200059/1‑147 (MQ=255)
               tGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATgg                                                                                                                         <  3:84629/147‑1 (MQ=255)
                 cccGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCCCCCCGGCCGGTGCGGTTAAAGGCAAgg                                                                                                                         >  8:25904/1‑145 (MQ=255)
                  ccGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATggg                                                                                                                        <  1:221527/145‑1 (MQ=255)
                   cGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGgag                                                                                                                     <  1:57399/147‑1 (MQ=255)
                    gAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCCCCCCCGCCCGGGGCGGTCAAGGCCAgggggg                                                                                                                     >  6:128055/1‑144 (MQ=255)
                          cACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGaacaa                                                                                                              <  3:227764/147‑1 (MQ=255)
                            cGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCAGTTCAAGGCATGGGGAGAGAacaaca                                                                                                            <  4:129364/147‑1 (MQ=255)
                             gTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCCGCCCGGGCGGGTGAAGGCAAGGGGTGAGGAcaacaa                                                                                                           >  8:225009/1‑147 (MQ=255)
                                              cATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCCCCGGCCCCGGCGGTTAAAGGAAATGGGACTCAAAAAAAAGTCGAGATCGTcactgt                                                                                          >  6:208955/1‑142 (MQ=255)
                                                 tCAGAGATCAAATACCGGTTATTCCAACATTTATTTTGTGGATGCCACGCACTCGACCCTATTGTGTGCCCACAACCCCTGCGCCGCACCGGCGGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGa                                                                                       <  7:254340/147‑1 (MQ=255)
                                                                 gggTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGc                                                                       <  3:215035/147‑1 (MQ=255)
                                                                          aaCATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACg                                                               <  8:187780/146‑1 (MQ=255)
                                                                          aaCATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCCGCCCGGGCGGTTAAATGGAAGGGGTGAGAAAAAAAAATAGAGATAGTCAAATTTGAACCATTAATAAATAGCCCCCCACgc                                                              >  8:137023/1‑146 (MQ=255)
                                                                           aCATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCCGTGCGGGTGAAGGCAAGGGGTGAGAACAACAATAGGACAACGTACATTTTGAACCATTAAAACATAAcccccccccct                                                             >  8:70088/1‑142 (MQ=255)
                                                                            cATTTAATTTGTGGAGGCCCCGCACTCGACCCTATTCTGGGGCCACAGCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTg                                                            <  1:31716/147‑1 (MQ=255)
                                                                             aTTTAATTTGTGGAGGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTgg                                                           <  5:80209/147‑1 (MQ=255)
                                                                                                    cTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTGAAAGGCATGGGGGGAGAACAACAAGTCGAGATCGTCCAGGTCGAACCATTAATACTTCGCCCCCCCCGGTGGGCCCCCGGCGCAAACCgcgcgcc                                    >  8:149354/1‑144 (MQ=255)
                                                                                                        aCCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGctgctg                                  <  2:19898/145‑1 (MQ=255)
                                                                                                                               cTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAg         <  3:226274/147‑1 (MQ=255)
                                                                                                                                       gcACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCTTCCATGTCGCACCATTCATACTTCGCCCCCCACAGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGc  >  5:172787/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
GCCGCGGCGGCCAGCTGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGC  >  NC_002947_CJ‑RC/4001628‑4001908

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

GATK/CNVnator alignment

N/A