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
A6 F21 I0 R1 2169 63.0 3952800 76.2% 3012033 62.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,155,607 T→A 7.6% intergenic (+118/+68) dacB → / ← obgE D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,155,6070TA 7.6% 181.4 / 6.5 66intergenic (+118/+68)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base T (40/21);  new base A (3/2);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATT  >  minE/2155538‑2155676
                                                                     |                                                                      
gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACata                                                                       <  1:1502857/71‑1 (MQ=39)
 tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                                      >  1:3163074/1‑71 (MQ=39)
   cgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                                                                    >  1:592977/1‑71 (MQ=39)
          tttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtc                                                              <  1:112347/70‑1 (MQ=255)
            tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTc‑tct                                                          >  1:160016/1‑71 (MQ=255)
                ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTa                                                      <  1:757168/71‑1 (MQ=255)
                  ggATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAg                                                     <  1:2820287/70‑1 (MQ=255)
                   gATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                                                    >  1:378758/1‑71 (MQ=11)
                     tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCt                                                 >  1:2547970/1‑71 (MQ=255)
                      gcggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCt                                                 >  1:1473641/1‑70 (MQ=255)
                        ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGa                                               >  1:1435700/1‑70 (MQ=255)
                        ggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGAt                                              <  1:2525497/71‑1 (MQ=255)
                        ggcgTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                                               >  1:3172492/1‑71 (MQ=14)
                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                                                                    <  1:204757/48‑1 (MQ=255)
                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATaa                                            <  1:1881545/71‑1 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGc                                                   >  1:3349055/1‑60 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCg                                         >  1:1649692/1‑70 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCg                                         >  1:2966096/1‑70 (MQ=255)
                              aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCg                                         >  1:1818822/1‑70 (MQ=255)
                               aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCCACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGtc                                       >  1:900512/1‑70 (MQ=255)
                               aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCg                                         >  1:505943/1‑69 (MQ=255)
                                  ccTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAgcg                                    >  1:1178601/1‑71 (MQ=255)
                                   cTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCaa                                                                 >  1:634349/1‑42 (MQ=40)
                                   cTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCgca                                       <  1:1723047/68‑3 (MQ=11)
                                     tATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCAt                                 >  1:3104888/1‑71 (MQ=255)
                                      aTCCGGCCTACAAAATCGTGCAAATTCAACATATTGCTATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCAt                                 >  1:828670/1‑70 (MQ=255)
                                         cGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                             <  1:669121/71‑1 (MQ=255)
                                           gCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                           >  1:1363920/1‑71 (MQ=255)
                                            ccTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                            >  1:627270/1‑69 (MQ=255)
                                            ccTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                           >  1:1067304/1‑70 (MQ=255)
                                            ccTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                           >  1:1387815/1‑70 (MQ=255)
                                             cTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTg                                                        <  1:2097864/40‑1 (MQ=255)
                                             cTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaa                          >  1:1006971/1‑69 (MQ=2)
                                                cAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGt                                                       >  1:2639210/1‑38 (MQ=255)
                                                cAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGAt                                              >  1:15816/1‑47 (MQ=255)
                                                cAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt                      >  1:380661/1‑71 (MQ=255)
                                                cAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGcaat                         <  1:2393310/69‑4 (MQ=0)
                                                 aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATa                                             <  1:3633402/47‑1 (MQ=255)
                                                 aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCg                                         <  1:3844684/51‑1 (MQ=255)
                                                 aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCa                               <  1:962781/61‑1 (MQ=255)
                                                 aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCa                               <  1:3614853/61‑1 (MQ=255)
                                                 aaaaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCa                               <  1:1851692/61‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                   >  1:983841/1‑71 (MQ=255)
                                                   aaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                   >  1:227844/1‑71 (MQ=255)
                                                   aaTCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc                   >  1:2928902/1‑71 (MQ=255)
                                                     tCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt                 <  1:3388900/71‑1 (MQ=255)
                                                     tCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt                 >  1:2420039/1‑71 (MQ=255)
                                                     tCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGt                 >  1:725990/1‑71 (MQ=255)
                                                      cGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGtt                >  1:1361377/1‑71 (MQ=255)
                                                      cGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGtt                >  1:2834110/1‑71 (MQ=255)
                                                         gCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTat             >  1:36949/1‑71 (MQ=255)
                                                           aaaTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTa                                       <  1:1476321/43‑1 (MQ=255)
                                                           aaaTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGttt               <  1:265077/67‑1 (MQ=255)
                                                            aaTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAtt                       >  1:2471532/1‑58 (MQ=255)
                                                            aaTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGttt               <  1:2116251/66‑1 (MQ=255)
                                                            aaTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat          >  1:354385/1‑71 (MQ=255)
                                                             aTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCg                                         <  1:3824181/39‑1 (MQ=255)
                                                             aTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat          >  1:2639318/1‑70 (MQ=255)
                                                             aTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat          >  1:643573/1‑70 (MQ=255)
                                                             aTTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat          >  1:641304/1‑70 (MQ=255)
                                                              ttCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGTCGTTTAtcatca        <  1:3049201/71‑1 (MQ=255)
                                                               tCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGc                                          >  1:2589245/1‑36 (MQ=255)
                                                               tCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTa                                       >  1:1665264/1‑39 (MQ=255)
                                                                 aaCATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGt      <  1:3788766/70‑1 (MQ=255)
                                                                  aCATATTGCAATTC‑TCTAGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCACTg     >  1:2637782/1‑70 (MQ=255)
                                                                   cATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatcat          <  1:1192717/64‑1 (MQ=255)
                                                                   cATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAt   >  1:1179388/1‑71 (MQ=255)
                                                                   cATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAt   >  1:3320652/1‑71 (MQ=255)
                                                                    atatTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAtt  <  1:2267692/71‑1 (MQ=255)
                                                                    atatTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGAtt  <  1:522826/71‑1 (MQ=255)
                                                                     |                                                                      
GTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTC‑TCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATCATCAGTGATT  >  minE/2155538‑2155676

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

GATK/CNVnator alignment

N/A