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 F28 I1 R1 155 60.3 2039550 99.3% 2025273 138.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,040,562 T→C V27A (GTC→GCC)  msrQ → periplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,5620TC98.0% 165.8 / ‑4.9 49V27A (GTC→GCC) msrQperiplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (22/26);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACGCAAAACAGCCGTTATT  >  NC_000913/2040427‑2040700
                                                                                                                                       |                                                                                                                                          
aCGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTgg                                                                                                                                        >  1:733772/1‑140 (MQ=255)
aCGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTgg                                                                                                                                        >  1:201892/1‑140 (MQ=255)
aCGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTgg                                                                                                                                        >  1:654530/1‑140 (MQ=255)
aCGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTgg                                                                                                                                        >  2:222056/1‑140 (MQ=255)
aCGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTgg                                                                                                                                        >  2:825859/1‑140 (MQ=255)
aCGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTgg                                                                                                                                        >  2:586438/1‑140 (MQ=255)
      cccAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGatc                                                                                                                                  >  1:159811/1‑140 (MQ=255)
        cAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGatcaa                                                                                                                                <  2:97056/140‑1 (MQ=255)
            tGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCAc                                                                                                                            <  2:201892/140‑1 (MQ=255)
             ggCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACg                                                                                                                           <  1:222056/140‑1 (MQ=255)
                aTCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTg                                                                                                                        >  2:396037/1‑140 (MQ=255)
                    cTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACt                                                                                                                    <  2:451577/140‑1 (MQ=255)
                     tGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTg                                                                                                                   >  2:364001/1‑140 (MQ=255)
                         tCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTg                                                                                                               >  2:834556/1‑140 (MQ=255)
                           gTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGcc                                                                                                             <  1:834556/140‑1 (MQ=255)
                               ccTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATc                                                                                                         <  1:508481/140‑1 (MQ=255)
                                 tGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATc                                                                                                         >  2:953729/1‑138 (MQ=255)
                                 tGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATc                                                                                                         <  1:953729/138‑1 (MQ=255)
                                      ttGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGaa                                                                                                  <  1:396037/140‑1 (MQ=255)
                                        gCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAg                                                                                                >  1:629016/1‑140 (MQ=255)
                                           ggAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGata                                                                                             <  2:654648/140‑1 (MQ=255)
                                             agaATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATAtt                                                                                           <  2:591195/140‑1 (MQ=255)
                                              gaATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGCCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTc                                                                                          <  2:763113/140‑1 (MQ=255)
                                               aaTTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCa                                                                                         <  1:825859/140‑1 (MQ=255)
                                                           cGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTgg                                                                             >  2:878554/1‑140 (MQ=255)
                                                            gTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGt                                                                            >  2:19245/1‑140 (MQ=255)
                                                              cTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCg                                                                          <  1:724948/140‑1 (MQ=255)
                                                                   aGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactg                                                                     <  2:267235/140‑1 (MQ=255)
                                                                    gCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgc                                                                    <  1:878554/140‑1 (MQ=255)
                                                                    gCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgc                                                                    <  1:364001/140‑1 (MQ=255)
                                                                               gACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAtttt                                                         >  1:718740/1‑140 (MQ=255)
                                                                                       tGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCg                                                 >  1:3473/1‑140 (MQ=255)
                                                                                       tGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCg                                                 <  2:629016/140‑1 (MQ=255)
                                                                                           aGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCt                                             <  1:388217/140‑1 (MQ=255)
                                                                                                  cTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAAt                                      <  1:586438/140‑1 (MQ=255)
                                                                                                     cATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAc                                   >  1:29129/1‑140 (MQ=255)
                                                                                                       tCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAccc                                 >  2:476333/1‑140 (MQ=255)
                                                                                                               gATTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCa                         >  2:736943/1‑140 (MQ=255)
                                                                                                                aTTGTTGCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCAc                        >  2:44693/1‑140 (MQ=255)
                                                                                                                 ttgttgCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACg                       <  1:19245/140‑1 (MQ=255)
                                                                                                                 ttgttgCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACg                       <  2:159811/140‑1 (MQ=255)
                                                                                                                 ttgttgCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACg                       <  1:755626/140‑1 (MQ=255)
                                                                                                                      gCCGTTTCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACg                  <  2:3473/140‑1 (MQ=255)
                                                                                                                          tttCTCTGGCTGGCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACGCaaa              >  1:761310/1‑140 (MQ=255)
                                                                                                                              tctGGCTGGCCTTGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGTACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACGCAAAACAg          >  2:892616/1‑140 (MQ=255)
                                                                                                                                 ggctggCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACGCAAAACAGCCg       <  1:265698/140‑1 (MQ=255)
                                                                                                                                    tggCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACGCAAAACAGCCGTTa    <  1:44693/140‑1 (MQ=255)
                                                                                                                                    tggCCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACGCAAAACAGCCGTTa    <  1:51409/140‑1 (MQ=255)
                                                                                                                                      gTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGCCACGCTACGCAAAACAGCCGTTAtt  >  2:352751/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                          
ACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCACGCTACGCAAAACAGCCGTTATT  >  NC_000913/2040427‑2040700

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

N/A