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
A2 F13 I0 R1 203 85.6 3649704 98.8% 3605907 146.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 4,442,726 A→C 27.6% intergenic (‑89/+49) nicP‑II ← / ← nicT porin‑like protein/metabolite transport protein NicT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,442,7260AC27.6% 26.9 / 2.5 32intergenic (‑89/+49)nicP‑II/nicTporin‑like protein/metabolite transport protein NicT
Reads supporting (aligned to +/- strand):  ref base A (8/12);  new base C (6/3);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGCCGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGG  >  NC_002947_CJ‑RC/4442583‑4442872
                                                                                                                                                |                                                                                                                                                  
aaCCGCCGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCACGCgg                                                                                                                                               >  3:147595/1‑149 (MQ=255)
 aCCGCCGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCggg                                                                                                                                              <  7:306858/149‑1 (MQ=255)
  ccgccgGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCCGGCgggg                                                                                                                                             >  5:260986/1‑148 (MQ=255)
          cAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGccc                                                                                                                                     >  8:2865/1‑149 (MQ=255)
                    gACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCCGCCCCcctggcctg                                                                                                                           >  7:145475/1‑149 (MQ=255)
                       gaagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGGCCGGGGCGCccccgcccccccgccccccg                                                                                                                         >  5:49139/1‑145 (MQ=255)
                                cgcgGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGaa                                                                                                               <  2:94184/149‑1 (MQ=255)
                                      gcAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCCGGCGGGCCCCCGCCCCCCTGCCCAGCCGTTGCGTGAAgcgcag                                                                                                         >  6:401683/1‑147 (MQ=255)
                                          tGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCCGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGcc                                                                                                     >  8:132411/1‑149 (MQ=255)
                                              aCCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCATTGCATGGCCACCGGCGTGCGGGTTGAATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTa                                                                                                 <  8:64038/149‑1 (MQ=255)
                                               ccTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGGCCGCCGGGCCCCGCCCCCCCGCCCCGCCGGTTGGGGGAGCGCGGGGCCGcgat                                                                                                >  3:136239/1‑146 (MQ=255)
                                                 ttGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcc                                                                                                                                  >  5:103702/1‑113 (MQ=255)
                                                 ttGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcc                                                                                                                                  <  6:103702/113‑1 (MQ=255)
                                                    cATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACGGGGGTGGGGGTTGCATCTAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGtt                                                                                           <  8:460209/149‑1 (MQ=255)
                                                     ttcgggCTACTCCTGTATTTGGGCAACCGCAACCGGGCGCGGCGGCACTGCTTGGCCACGGGGGGGGGGGTTGAATCGAGAGTGGGCGGCCGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTg                                                                                          <  4:110464/146‑1 (MQ=255)
                                                     aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCCGGCGGGGCCCCGCCCCCCAGCCCAGCCGGTGGGGGTAGCGCGGGGGCCCGAGCGGGTg                                                                                          >  3:205325/1‑149 (MQ=255)
                                                           cTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGc                                                                                    <  6:87964/149‑1 (MQ=255)
                                                                      tttGGGCAAAAGAAACCCGGCCCGGCGGCCATGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCCCCCCCCCGGCCCGGCCCTTGGGGGAAACGCGGGGGCCCGTGGGGTTGACCACCATGGGCggg                                                                           >  4:56960/1‑146 (MQ=255)
                                                                               caccAACCCGGCGCGGCGGCACTGCAGGGCCACGGGCGGGGGGGTTGATTCGACCGAGGGCGGCAGGGGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGAc                                                                <  1:22367/146‑1 (MQ=255)
                                                                               cAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGAc                                                                <  2:148027/149‑1 (MQ=255)
                                                                                      ccGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCACGCCGGCCCCCGCCCCCCAGCCCAGCCCTTGGGTGAAGCCCTGGGCCGCTACCCGTTGACAAGCTTGGCCCGCAAGGCCATGAACACGAag                                                         >  8:16358/1‑149 (MQ=255)
                                                                                             gggggcACTGCAGGGCCACGGGCGTGGGGTTTGCATCGAGTGTGGTCGGCTGGGGGGCCCGGGCCGGCCGGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                  <  4:147595/146‑1 (MQ=255)
                                                                                             ggcggcACTGCTTGGCCCCGGGGGTGGGGGTTGAATCGAGTGAGGGGGGCCGGGGGGCCCCGGCCGGCGTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                  <  4:328768/149‑1 (MQ=255)
                                                                                                   aCTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCa                                            <  1:4930/149‑1 (MQ=255)
                                                                                                   aCTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCa                                            <  8:391369/149‑1 (MQ=255)
                                                                                                                  gCGTGCGGGTTCCATCGATAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                  <  7:256299/128‑1 (MQ=255)
                                                                                                                               aTCGAGAGAGGTCGGCAGGCGGGCCCCCGCCCGCCGGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCGAGCGGTTGACCAGCGTGGCCGGCAAGGCGATGACCACGAAGCAGCGGAGGGACAGGCAGCCGGCCAACACCCAcagcgcagc                >  8:207461/1‑149 (MQ=255)
                                                                                                                                  gtgtgagGTCGGCTGGCGGGCCCGGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGt              <  2:237365/144‑1 (MQ=255)
                                                                                                                                             ggCTGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCg   <  6:260986/149‑1 (MQ=255)
                                                                                                                                             ggCCGGGGGGCCCGGCCCCGCCGGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCg   <  5:145010/149‑1 (MQ=255)
                                                                                                                                              gCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCgg  <  6:312413/149‑1 (MQ=255)
                                                                                                                                              gCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCgg  >  2:372335/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AACCGCCGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTTGCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGG  >  NC_002947_CJ‑RC/4442583‑4442872

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