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
A0 F0 I1 R1 52 69.5 2964196 98.4% 2916768 145.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 278,744 +A intergenic (‑239/+75) PP_0224 ← / ← sctC DszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC278,7421.A83.3% 60.3 / 8.4 24intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base A (2/18);  total (5/19)
Fisher's exact test for biased strand distribution p-value = 1.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAA  >  NC_002947_CJ‑RC/278595‑278857
                                                                                                                                                      |                                                                                                                       
aTGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGt                                                                                                                          >  4:258906/1‑148 (MQ=255)
  gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCAGGTaaa                                                                                                                      >  8:162982/1‑148 (MQ=255)
         cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTc                                                                                                              <  2:21027/149‑1 (MQ=255)
           gCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCt                                                                                                            <  6:97353/149‑1 (MQ=255)
               aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTac                                                                                                          <  4:344363/147‑1 (MQ=255)
                     aCAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGc                                                                                                  <  1:333127/149‑1 (MQ=255)
                        gACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTc                                                                                               <  7:206506/149‑1 (MQ=255)
                          cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCaa                                                                                             <  5:149512/149‑1 (MQ=255)
                           gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAg                                                                                            <  3:9701/149‑1 (MQ=255)
                              tgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGTAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGAt                                                                                          <  7:317885/148‑1 (MQ=255)
                                cTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCt                                                                                        <  7:290488/148‑1 (MQ=255)
                                 tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGGTAA‑‑ACCCCTCCGAGACCCAAGACGCAAGATAtgtg                                                                                     >  8:357891/1‑149 (MQ=255)
                                   gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                      <  6:198825/147‑1 (MQ=255)
                                        gtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGt                                                                               <  4:284030/149‑1 (MQ=255)
                                              cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAg                                                                          <  1:106815/148‑1 (MQ=255)
                                              cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGc                                                                         <  2:119778/149‑1 (MQ=255)
                                                       taCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAc                                                               <  6:230037/148‑1 (MQ=255)
                                                         cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAccc                                                             <  5:233863/148‑1 (MQ=255)
                                                                gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTTAA‑‑CCCGCGCCAAAACAAATTCCGCAAGATGTGTGCACGGTAGGAGCGGGTTTTACCGCGAAAgg                                                     >  8:131166/1‑149 (MQ=255)
                                                                     ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                    >  4:151308/1‑130 (MQ=255)
                                                                     ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                    <  3:151308/130‑1 (MQ=255)
                                                                                   aactaactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTTAA‑‑ACCCCACCAACACAAAATCCGCAAGATAAGTGGACAGTATGAGCGGGGTTTACCCCGAAAAGGGCGGGACAGGgtttatag                                  >  8:186761/1‑144 (MQ=255)
                                                                                                            aaGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCAGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgc         <  7:162982/149‑1 (MQ=255)
                                                                                                                    tttGTTACCGGGGGGCCGGGCTCTCGGGGGTAAACCCCGCCCAAACACCCCCCCCTACAGTTGTTTGCATTGTCGGCGGGGTTTACCCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGaa  <  2:275363/148‑1 (MQ=255)
                                                                                                                                              gcgGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATc                                                                                         >  2:351277/1‑37 (MQ=17)
                                                                                                                                              gcgGGTAAA‑‑CCCGCTCCTACACACATGCCTCAAGATc                                                                                         <  1:351277/37‑1 (MQ=17)
                                                                                                                                                      |                                                                                                                       
ATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAA‑‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAA  >  NC_002947_CJ‑RC/278595‑278857

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAA  >  NC_002947_CJ‑RC/278615‑278889
                                                                                                                                  |                                                                                                                                                 
GACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCC                                                                                                                               <  NB501016:282:H2V7GAFX2:1:21306:19108:1768/151‑1 (MQ=60)
   AGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCA                                                                                                                            <  NB501016:282:H2V7GAFX2:4:21403:19176:10704/151‑1 (MQ=60)
     ACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAG                                                                                                                          <  NB501016:282:H2V7GAFX2:3:11605:3962:11379/151‑1 (MQ=60)
      CGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGA                                                                                                                         <  NB501016:282:H2V7GAFX2:2:11102:21999:20080/151‑1 (MQ=60)
         GTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGTAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATC                                                                                                                       <  NB501016:282:H2V7GAFX2:4:21511:3742:18185/150‑1 (MQ=60)
           GCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTG                                                                                                                     <  NB501016:282:H2V7GAFX2:4:21506:9652:16710/150‑1 (MQ=60)
                         CCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGC                                                                                                       <  NB501016:282:H2V7GAFX2:1:11208:2692:17558/150‑1 (MQ=60)
                                    ACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCG                                                                                            <  NB501016:282:H2V7GAFX2:3:21410:6927:4690/150‑1 (MQ=60)
                                                TATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTT                                                                                                  <  NB501016:282:H2V7GAFX2:2:11306:2065:6012/132‑1 (MQ=60)
                                                                                       GAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCAGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCG                                        <  NB501016:282:H2V7GAFX2:4:11607:21481:5455/151‑1 (MQ=60)
                                                                                                                         CGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCT                                                                                                                      <  NB501016:282:H2V7GAFX2:1:21309:7226:15203/39‑1 (MQ=20)
                                                                                                                              gtaAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAA  >  NB501016:282:H2V7GAFX2:2:21211:15226:11813/4‑150 (MQ=60)
                                                                                                                                  |                                                                                                                                                 
GACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAA  >  NC_002947_CJ‑RC/278615‑278889

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