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
A11 F28 I1 R1 228 65.7 2434974 91.7% 2232871 138.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,913,058 T→C N587D (AAC→GAC)  prc ← tail‑specific protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,913,0580TC97.6% 136.9 / ‑4.9 41N587D (AAC→GAC) prctail‑specific protease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (24/16);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAGCGCG  >  NC_000913/1912923‑1913179
                                                                                                                                       |                                                                                                                         
tAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTgg                                                                                                                       >  2:1016624/1‑140 (MQ=255)
    gCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACt                                                                                                                   >  2:1042268/1‑140 (MQ=255)
       ttcgttcAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACt                                                                                                                   <  2:596361/137‑1 (MQ=255)
       ttcgttcAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACt                                                                                                                   >  1:596361/1‑137 (MQ=255)
        tcgttcAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAgg                                                                                                               <  1:1004648/140‑1 (MQ=255)
                              tCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTa                                                                                         >  1:236628/1‑140 (MQ=255)
                               cGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTAt                                                                                        <  2:999072/140‑1 (MQ=255)
                                    catcTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTcc                                                                                   <  1:1016624/140‑1 (MQ=255)
                                      tcTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCtt                                                                                 >  1:409990/1‑140 (MQ=255)
                                          cATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGTTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTcagc                                                                             <  2:389049/140‑1 (MQ=255)
                                               ttCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTc                                                                        >  1:803615/1‑140 (MQ=255)
                                               ttCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTc                                                                        <  2:236628/140‑1 (MQ=255)
                                                tctcTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTcc                                                                       <  1:1072155/140‑1 (MQ=255)
                                                    tttCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGtt                                                                   >  1:1025162/1‑140 (MQ=255)
                                                    tttCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGtt                                                                   >  1:180689/1‑140 (MQ=255)
                                                     ttCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTc                                                                  >  1:627870/1‑140 (MQ=255)
                                                                  cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAGGGCCGTTaaa                                                     >  2:904418/1‑140 (MQ=255)
                                                                  cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTaaa                                                     >  2:863104/1‑140 (MQ=255)
                                                                  cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTaaa                                                     >  2:379519/1‑140 (MQ=255)
                                                                  cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTaaa                                                     >  1:836806/1‑140 (MQ=255)
                                                                        tCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCt                                               <  1:863104/140‑1 (MQ=255)
                                                                          agagAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGa                                             >  2:124352/1‑140 (MQ=255)
                                                                           gagaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGAt                                            <  1:379519/140‑1 (MQ=255)
                                                                            agaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAAtc                                                   >  2:107279/1‑132 (MQ=255)
                                                                            agaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAAtc                                                   <  1:107279/132‑1 (MQ=255)
                                                                            agaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGAtt                                           <  1:124352/140‑1 (MQ=255)
                                                                                 cGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTcaca                                      >  1:674722/1‑140 (MQ=255)
                                                                                        gcgcTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTcg                               <  2:409990/140‑1 (MQ=255)
                                                                                              gTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCAt                         <  2:1025162/140‑1 (MQ=255)
                                                                                              gTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCAt                         <  2:438385/140‑1 (MQ=255)
                                                                                               tCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATc                        <  1:1042268/140‑1 (MQ=255)
                                                                                                    cATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGc                   <  1:1124712/140‑1 (MQ=255)
                                                                                                      tAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTa                 >  1:564694/1‑140 (MQ=255)
                                                                                                           ttGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCa            >  1:337884/1‑140 (MQ=255)
                                                                                                              aaGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACgg         >  2:1122511/1‑140 (MQ=255)
                                                                                                              aaGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACgg         >  2:669208/1‑140 (MQ=255)
                                                                                                                 cgcgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAg      <  2:586528/140‑1 (MQ=255)
                                                                                                                   cgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAgcg    >  1:867328/1‑140 (MQ=255)
                                                                                                                   cgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAgcg    >  2:707024/1‑140 (MQ=255)
                                                                                                                   cgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAgcg    >  2:1098517/1‑140 (MQ=255)
                                                                                                                     cgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAGagcg  >  2:181189/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                         
TAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAGCGCG  >  NC_000913/1912923‑1913179

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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