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 F1 I0 R1 4577 57.4 3436498 77.7% 2670158 63.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,803,427 C→T 54.1% S80S (AGC→AGT)  dnaB → replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,803,4270CT54.1% 9.1 / 84.2 99S80S (AGC→AGT) dnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base T (17/36);  total (35/64)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GACGATTTTTACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTG  >  minE/2803359‑2803491
                                                                    |                                                                
gACGATTTTTACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGc                                                                  >  1:941319/1‑69 (MQ=255)
       tttACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCg                                                           >  1:1145690/1‑69 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:3311690/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2088361/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2081698/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2258977/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:1986282/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2321793/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2325661/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:2747154/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:296094/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:1339209/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:1147586/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:754854/69‑1 (MQ=255)
              ccGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGAt                                                    <  1:575617/69‑1 (MQ=255)
                   cacaCCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTActc                                               >  1:551182/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCtt                                             >  1:1244934/1‑68 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:2706027/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1020280/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:2078331/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1184790/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:614211/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1814971/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:248691/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:2627235/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1300710/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:3052250/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:31330/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1357075/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:401958/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:1338827/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:3396803/1‑69 (MQ=255)
                      acCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTg                                            >  1:3428960/1‑69 (MQ=255)
                       ccGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTg                                            >  1:3034850/1‑68 (MQ=255)
                       ccGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTg                                            >  1:2065100/1‑68 (MQ=255)
                         gTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:960513/1‑69 (MQ=255)
                          tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:2027805/69‑1 (MQ=255)
                          tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:3112402/69‑1 (MQ=255)
                          tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:2466577/69‑1 (MQ=255)
                          tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:885139/69‑1 (MQ=255)
                          tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGa                                        <  1:2684307/69‑1 (MQ=255)
                          tCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTAGCGGa                                        <  1:2254475/69‑1 (MQ=255)
                              atCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATc                                     <  1:2884277/68‑1 (MQ=255)
                              atCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCg                                    >  1:3251129/1‑69 (MQ=255)
                               tCTTTACTTAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:390423/1‑69 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:3348489/1‑63 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:186613/1‑63 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:3142452/1‑63 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:356790/1‑63 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCgg                                         >  1:657031/1‑63 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCg                                    >  1:1267321/1‑68 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:2263738/1‑69 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:2673569/1‑69 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:1169247/1‑69 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:245036/1‑69 (MQ=255)
                               tCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGc                                   >  1:2451892/1‑69 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTTATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:3269905/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2553773/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:337460/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:132900/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:3421255/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:315000/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2363066/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1375952/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1423375/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2934170/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:650432/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2864139/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2857353/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:279841/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2771554/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1533163/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1566619/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1264739/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1215531/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1251465/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:2633081/68‑1 (MQ=255)
                                          aTGGCGCGTTTGCAGGAAAGCGGTAGTCCAATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1426088/68‑1 (MQ=255)
                                           tGGCGCGTTTGCAGGAAAGCGGTAGTCCTATCGATCTGATTACCCTTGCGGAATCGCTGGAACGCCa                         <  1:1280100/67‑1 (MQ=255)
                                           tGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTgg                                <  1:2631724/60‑1 (MQ=255)
                                                   ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:76406/69‑1 (MQ=255)
                                                   ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1012141/69‑1 (MQ=255)
                                                   ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:3190555/69‑1 (MQ=255)
                                                   ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1633636/69‑1 (MQ=255)
                                                   ttGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:148451/69‑1 (MQ=255)
                                                    tGCAGGAAAGCGGTAGCCCTATCGATGTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:2731301/68‑1 (MQ=255)
                                                    tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1372530/68‑1 (MQ=255)
                                                    tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:2211413/68‑1 (MQ=255)
                                                    tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:2684050/68‑1 (MQ=255)
                                                    tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:3052348/68‑1 (MQ=255)
                                                    tGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1324618/68‑1 (MQ=255)
                                                       aGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt            <  1:2779217/68‑1 (MQ=255)
                                                       aGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt            <  1:2600957/68‑1 (MQ=255)
                                                        ggAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGAt            <  1:2605038/67‑1 (MQ=255)
                                                               ggTAGGCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTc               <  1:1972487/57‑1 (MQ=255)
                                                                gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:37620/69‑1 (MQ=255)
                                                                gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:2689008/69‑1 (MQ=255)
                                                                gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:1609581/69‑1 (MQ=255)
                                                                gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:2524284/69‑1 (MQ=255)
                                                                gTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCggtg  <  1:3357789/69‑1 (MQ=255)
                                                                    |                                                                
GACGATTTTTACACCCGCCCACACCGTCATATCTTTACTGAAATGGCGCGTTTGCAGGAAAGCGGTAGCCCTATCGATCTGATTACTCTTGCGGAATCGCTGGAACGCCAGGGGCAACTCGATAGCGTCGGTG  >  minE/2803359‑2803491

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

GATK/CNVnator alignment

N/A