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 F14 I0 R1 248 61.6 2601574 98.7% 2567753 146.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 5,754,561 A→C 25.2% intergenic (‑184/‑287) thiI ← / → glnA tRNA sulfurtransferase/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC5,754,5610AC25.2% 46.2 / 6.7 32intergenic (‑184/‑287)thiI/glnAtRNA sulfurtransferase/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (15/6);  new base C (2/5);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 7.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCACGCAAACGCCCCCTT  >  NC_002947_CJ‑RC/5754413‑5754704
                                                                                                                                                      |                                                                                                                                                
ggggCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGa                                                                                                                                                  <  8:44166/149‑1 (MQ=255)
  ggCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCTTTGCACCGCCTGTAGgcgc                                                                                                                                                <  6:289705/149‑1 (MQ=37)
  ggCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACGGGCTTTGCACCGCCTGTAGgcgc                                                                                                                                                <  4:149587/149‑1 (MQ=38)
    cgcgAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCgg                                                                                                                                             <  4:192178/149‑1 (MQ=255)
       gAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTTAAACCGCGCCAACAGGCGGTGCAAAGCCCGTAGGGGCGGGtt                                                                                                                                          >  6:66323/1‑149 (MQ=255)
       gAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGtt                                                                                                                                          <  2:272242/149‑1 (MQ=255)
          ttATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAc                                                                                                                                       <  4:24744/149‑1 (MQ=255)
                cGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGGTAAAACCCCCCCGAAAGGCGGGCTAAACCCCCAAGAGCGGGttttcacccg                                                                                                                                  >  1:36584/1‑141 (MQ=255)
                  taaaTTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGAAAAGCCCGTAGGAGCGGGTTTACCCGCGAag                                                                                                                               <  7:170261/148‑1 (MQ=255)
                          tCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTTAAACCCCCCCAAAAGGCGGGGGAAAACCCCGAGGGGCGGGGTTTACCCCGCAGAGGccc                                                                                                                         >  4:273849/1‑147 (MQ=255)
                            aGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAAACCGCGCCGACAGGCGGTGCAAAGCCCGTAGGAGCGGGGTTTCCCGCGAAGAGGCCCTAg                                                                                                                      >  8:235539/1‑148 (MQ=255)
                                   cacaaaCTTTTTATCATACGCTTTTTTAAAAATGCAGCGCCCTGTTTTGGCCGCGCGTTTCCGGGCCGTTTCGGGGTTAACCCCGCGCAAACCGGGGGTCAAAAGCCGTTCGGCGGGGGTTTACCCGCGAAGAGGCCCTAGAACCTgg                                                                                                               <  8:210936/144‑1 (MQ=255)
                                    aCCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGGTAAACCCCCCCGAAAGGCGGGGGCAAGCCCCGTGGGGCGGGGTTTCACGCGCAGAGGGCCGAGAACCAGGcc                                                                                                             >  1:86037/1‑148 (MQ=255)
                                       aaCTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTTAAAACCCCCCGAAAAGCGGTGGAAAGCCCCTAGGAGCGGGGTTTACCCCGAAGAGGCCCAATAACCAGGccccc                                                                                                          >  1:220661/1‑145 (MQ=255)
                                                             taTAAATCCAGCGGCCAGGTTTGGCAGCGCCTGTCCCGGCCATTTCGGGGGTAACCCCGCGCATACCGGCGGTGAAAAGCCCTTAGGAGGGTGTTAACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTAt                                                                                     <  8:184226/147‑1 (MQ=255)
                                                                       gcggCCCGGTTTTGCCGCGCCTGTCCCGGCCCTTTGGCGGTTAACCCCGCGCCAACCGGGGGGCAAAAGCCGTTAGTAGCGGGTTAACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGacac                                                                          <  2:55209/147‑1 (MQ=255)
                                                                               gTTTTGCCCCGCCTGTCCCGCCCCTTGCGGGGTAAACCCCGCGCAAACCGGGGGGGAAATGCCCGTCGGCGCGGGTTAACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCaa                                                                    <  7:331538/147‑1 (MQ=255)
                                                                                    gCAGCGCGTGTTCCGGCCATTTCGGGGGTAACCCCGGGCTTACGGGGGGTGAAATGCCGGTAGTCGGGGTTTTACCCGCGAAGAGGCCCTAGAACCTTTATCACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACg                                                              <  2:36584/148‑1 (MQ=255)
                                                                                     cAGCGCGTGTTCCGCCCCTTTCGGGGTAAACCCCGCGCCTACTGGCGGTTAAAAGCCCGTTGGTGCGGTTTAACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGc                                                             <  8:58507/148‑1 (MQ=255)
                                                                                     cAGCGCCTGTTCCGGCCCTTTCGTGGGTAAAACCCCCCCAACAAGCGGTGGAAAGCCCCGAGGAGCGGGGTTAACCGCGCAGATGGCCTAGAAACTGGCCCACCACAAACCCCACCCCCCCAATTGGTGCGACAAGATCAAAAAACGc                                                             >  6:109591/1‑148 (MQ=255)
                                                                                              tcccGCCCCTTGCGCGGTAAACCCCGCGCAAACCGGGGGTCAAACGCCGTTTGTCGGGGGTTAACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAgg                                                    <  7:52384/146‑1 (MQ=255)
                                                                                                              gggTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGc                                    >  6:230905/1‑149 (MQ=255)
                                                                                                               ggTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGc                                    >  1:207330/1‑148 (MQ=255)
                                                                                                                gTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCt                                   >  8:328470/1‑148 (MQ=255)
                                                                                                                gTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCt                                   >  6:65282/1‑148 (MQ=255)
                                                                                                                   aaccccGCGCATACCGGCGGTGAAAAGCCGGACGGCGCGGGTTAACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTaa                                 <  7:235539/145‑1 (MQ=255)
                                                                                                                    aCCCGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAAcc                               >  8:208050/1‑148 (MQ=37)
                                                                                                                       ccGCTCCTACGGGCTTTGCACCGCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCaa                            >  2:270641/1‑149 (MQ=37)
                                                                                                                              tACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTAACCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCaaa                     >  4:166536/1‑149 (MQ=255)
                                                                                                                                  ggCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCa                 >  3:236006/1‑149 (MQ=255)
                                                                                                                                   gCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCAc                >  2:83788/1‑149 (MQ=255)
                                                                                                                                      gTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCACGc              >  5:136530/1‑148 (MQ=255)
                                                                                                                                       tGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCACGCa             >  1:312417/1‑148 (MQ=255)
                                                                                                                                          aaaGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCACGCAAACg         >  3:134444/1‑149 (MQ=255)
                                                                                                                                                 gTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACACGCTAACCCAACGCCAAATGCACGCAAACGCCCCCtt  >  2:305374/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GGGGCGCGAATTATATCGGAAATTGCTCAGGATTTGACCAACTTTTGATCAGAAGCTTTGAGATAAATGCAGAGGCCAGGTTTTGCAGCGCCTGTTCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGCGGTGCAAAGCCCGTAGGAGCGGGTTTACCCGCGAAGAGGCCCTAGAACCTGGCACACCACTCACGCCACGCACCCTATTGGTGCGACACGATCAAAAAACGCATCGTAAGGGTGCACTTTCACCCGCTAACCCAACGCCAAATGCACGCAAACGCCCCCTT  >  NC_002947_CJ‑RC/5754413‑5754704

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