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
A18 F53 I0 R2 221 66.0 2680176 94.9% 2543487 126.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,980 T→C 100% K63K (AAA→AAG U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,9800TC92.7% 107.9 / ‑4.8 33K63K (AAA→AAGU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (4/26);  minor base G (3/0);  total (7/26)
Fisher's exact test for biased strand distribution p-value = 6.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCG  >  CP006881/3724851‑3725110
                                                                                                                                 |                                                                                                                                  
cTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAACGCGCGGCGGCTTTCGCTTTCTGCTTttcttc                                                                                                                            <  6:158514/138‑1 (MQ=255)
cTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTttcttc                                                                                                                            <  4:11853/138‑1 (MQ=255)
       cTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTca                                                                                                                    <  2:181213/139‑1 (MQ=255)
       cTGGGATTTTACCTACTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACTTCTTCTTTTATATAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTca                                                                                                                    <  1:109814/139‑1 (MQ=255)
          ggATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAg                                                                                                                 <  1:53854/139‑1 (MQ=255)
           gATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGc                                                                                                                <  4:28819/139‑1 (MQ=255)
           gATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGc                                                                                                                <  4:24512/139‑1 (MQ=255)
                      ctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTGTGGCGTTCGTGTGTTAGTCCTCGTTGACAGCCCTCCGCCat                                                                                                      >  8:79951/1‑137 (MQ=255)
                      ctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGGTTTGTGCTTTTATTCCTCATTCACAGCCCTCCGCCat                                                                                                      >  7:65004/1‑137 (MQ=255)
                              cGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAAc                                                                                             <  5:99240/139‑1 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                                 <  1:194577/139‑2 (MQ=255)
                                             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCtgcg                                                                              <  6:172027/139‑1 (MQ=255)
                                             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTAGCTTTCTGGTTTTATTACACATTACCAGCCCAACGCCAATGATAAAACCAAACGAGTCCagcc                                                                              >  6:64368/1‑138 (MQ=255)
                                               ggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTATGCTTTTATTACTCATTCACAGCCATCCGCCAATCATAAAACGAAACCAGTCCTGCGGc                                                                            >  8:52751/1‑139 (MQ=255)
                                                     ctctctCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTATGGTTTTATTACTCATTAACACCCAACCGCCCTTCATAAAAAGAAACCATTCCTGCGGCCAGccc                                                                      >  8:162333/1‑137 (MQ=255)
                                                          tcCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagc                                                                   <  3:115503/137‑1 (MQ=255)
                                                            cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAcagccg                                                                 <  3:3581/137‑1 (MQ=255)
                                                                  tGGCGATTCCTTGATGTCTTGCAAACATCTTATTTTATAAAGCGCGGCGGCTTTCTCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCTTAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                          <  5:11982/138‑1 (MQ=255)
                                                                  tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                          <  7:52751/138‑1 (MQ=255)
                                                                  tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                          <  5:287798/138‑1 (MQ=255)
                                                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAACCGACACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                        <  1:274333/138‑1 (MQ=255)
                                                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                        <  5:208236/138‑1 (MQ=255)
                                                                    gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                        <  5:203995/138‑1 (MQ=255)
                                                                         tCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATg                                                  >  7:81861/1‑139 (MQ=255)
                                                                          ccTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATg                                                  <  2:77232/138‑1 (MQ=255)
                                                                              gatggCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                              <  5:128284/135‑1 (MQ=255)
                                                                              gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                              <  3:150048/138‑1 (MQ=255)
                                                                               aaGGCTGGCAAACATCTTCTTTTATAAAGCGCGCCGGCTTTCGCTTTCTGCTTTTCTTCCTCAATCCCAGCCATCCGCCATTCATAAAACGACACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                            <  2:77764/139‑1 (MQ=255)
                                                                                    ttGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTATTCCTGATTCACAGCCATCCGGCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGGCGCCCTTCGTTTGTTTATACGCtatacc                                       >  3:167996/1‑134 (MQ=255)
                                                                                                    ttATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCgg                          <  7:162333/136‑1 (MQ=255)
                                                                                                       taACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                    <  3:202022/139‑1 (MQ=255)
                                                                                                                    cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc       <  8:111415/139‑1 (MQ=255)
                                                                                                                         gCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCg  <  8:65004/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCG  >  CP006881/3724851‑3725110

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