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 1,289,523 A→G 61.0% E384E (GAA→GAG)  amn → AMP nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,289,5230AG61.0% 112.9 / 133.4 197E384E (GAA→GAG) amnAMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base A (51/26);  new base G (19/101);  total (70/127)
Fisher's exact test for biased strand distribution p-value = 1.02e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 197 total aligned reads displayed.

GAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGAAAG  >  minE/1289461‑1289588
                                                              |                                                                 
gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:3156835/1‑69 (MQ=255)
gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:3332996/1‑69 (MQ=255)
gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:640528/1‑69 (MQ=255)
gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:3053933/1‑69 (MQ=255)
gAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGt                                                             >  1:1953187/1‑69 (MQ=255)
   gAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTaa                                                               >  1:1802926/1‑63 (MQ=255)
   gAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTAc                                                               >  1:26095/1‑64 (MQ=255)
     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1012454/69‑1 (MQ=255)
     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:2268575/69‑1 (MQ=255)
     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1247531/69‑1 (MQ=255)
     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1321907/69‑1 (MQ=255)
     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:1770244/69‑1 (MQ=255)
     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:3435301/69‑1 (MQ=255)
     aGTCAAACAGCGGCTGCGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTc                                                        <  1:2319304/69‑1 (MQ=255)
     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:1524640/69‑1 (MQ=255)
     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:3073125/69‑1 (MQ=255)
     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:354940/69‑1 (MQ=255)
     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:2959379/69‑1 (MQ=255)
     aGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTc                                                        <  1:775260/69‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGTGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1737916/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCTGCTTCTGCACtt                                           <  1:3315229/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1908740/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2154727/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2017288/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2029621/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1636134/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:1622166/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:365932/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2506208/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:5318/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:682749/66‑1 (MQ=255)
                  ctgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACtt                                           <  1:2628394/66‑1 (MQ=255)
                   tgcGTACAGGTACTGTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTTCACtt                                           <  1:1075346/66‑1 (MQ=255)
                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:2981720/68‑1 (MQ=255)
                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:2736650/68‑1 (MQ=255)
                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:2102169/68‑1 (MQ=255)
                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:3222644/68‑1 (MQ=255)
                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:3231018/68‑1 (MQ=255)
                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:506752/68‑1 (MQ=255)
                       tacaggtactgTGGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGtt                                       <  1:711392/68‑1 (MQ=255)
                        actggtactgTGGTAACCACATATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:838903/69‑1 (MQ=255)
                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:1929461/69‑1 (MQ=255)
                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:2124452/69‑1 (MQ=255)
                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:497864/69‑1 (MQ=255)
                        actggtactgTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGtttt                                     <  1:2871956/69‑1 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCGCTTCGTTTTAACTTAAGcc                           >  1:3146857/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACtt                                >  1:2336130/1‑64 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:943183/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:2769584/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:759848/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:2999080/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:302948/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:2708359/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:1180815/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:1542916/1‑69 (MQ=255)
                                  tgGTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGcc                           >  1:157389/1‑69 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:2551453/68‑1 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:2520257/68‑1 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:2452039/68‑1 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:3126775/68‑1 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCg                          <  1:3250327/68‑1 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCgg                         >  1:2869034/1‑69 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCgg                         >  1:1456035/1‑69 (MQ=255)
                                    gTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCgg                         >  1:2386660/1‑69 (MQ=255)
                                         cacaGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTa                               >  1:1386825/1‑58 (MQ=255)
                                         cacaGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGt                    >  1:2816549/1‑69 (MQ=255)
                                               ttaCAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1841356/67‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCTTAGCAAt              <  1:336185/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2239282/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:10023/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:27777/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1282989/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1299632/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1375943/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:698149/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2903262/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2368235/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:438954/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1618785/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:1879842/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:3127571/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:3406297/69‑1 (MQ=255)
                                               tGACAGGAACTGGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2353981/69‑1 (MQ=255)
                                               tGACAGGAACTCGGAGTTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAt              <  1:2036692/69‑1 (MQ=255)
                                                gACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAtt             >  1:339913/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAAtt             >  1:475643/1‑66 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:1488633/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:3136985/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:903827/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:810327/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:1850029/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:1706085/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:521256/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:2752300/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:448008/1‑69 (MQ=255)
                                                   aGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGat          >  1:3385041/1‑69 (MQ=255)
                                                          gggAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGaaa   >  1:2658477/1‑69 (MQ=255)
                                                          gggAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGaaa   >  1:1657117/1‑69 (MQ=255)
                                                          gggAATTACGTTACTCAGCTTCTGCACTACGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGaaa   >  1:2926844/1‑69 (MQ=255)
                                                            gAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGAAAg  >  1:1609232/1‑68 (MQ=255)
                                                              |                                                                 
GAGGAAGTCAAACAGCGGCTACGTACTGGTACTGTGGTAACCACAGATGACAGGAACTGGGAATTACGTTACTCAGCTTCTGCACTTCGTTTTAACTTAAGCCGGGCCGTAGCAATTGATATGGAAAG  >  minE/1289461‑1289588

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

GATK/CNVnator alignment

N/A