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
A13 F1 I1 R1 778 68.7 3228504 93.4% 3015422 109.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,786,729 G→A pseudogene (333/2253 nt) ygaQ → uncharacterized protein, putative enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,786,7290GA96.5% 173.1 / ‑1.9 58pseudogene (333/2253 nt)ygaQuncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (30/26);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

AGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCT  >  NC_000913/2786594‑2786858
                                                                                                                                       |                                                                                                                                 
aGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGtt                                                                                                                                >  2:1087534/1‑139 (MQ=255)
 gAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGttt                                                                                                                               >  2:246417/1‑139 (MQ=255)
            taATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGAt                                                                                                                    >  2:1345949/1‑139 (MQ=255)
                  tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACgag                                                                                                                      >  2:907915/1‑131 (MQ=255)
                  tCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACgag                                                                                                                      <  1:907915/131‑1 (MQ=255)
                      ttAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGa                                                                                                         >  2:935533/1‑139 (MQ=255)
                        aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       <  2:930299/139‑1 (MQ=255)
                        aaaaaaCGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       <  1:1345949/139‑1 (MQ=255)
                                aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGAt                                                                                               >  2:742838/1‑139 (MQ=255)
                                aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGAt                                                                                               >  1:1123650/1‑139 (MQ=255)
                                aTTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGAt                                                                                               >  1:435557/1‑139 (MQ=255)
                                     aaaaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAg                                                                                          <  1:935533/139‑1 (MQ=255)
                                       aaTTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCg                                                                                        <  2:435557/139‑1 (MQ=255)
                                                       gATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATgg                                                                        >  2:199388/1‑139 (MQ=255)
                                                             aCTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGc                                                                  >  1:1565377/1‑139 (MQ=255)
                                                                tACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGa                                                               <  1:1087534/139‑1 (MQ=255)
                                                                    gagaGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACAt                                                           <  1:246417/139‑1 (MQ=255)
                                                                          tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  1:742838/139‑1 (MQ=255)
                                                                          tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  2:1123650/139‑1 (MQ=255)
                                                                          tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  2:468482/139‑1 (MQ=255)
                                                                          tGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  1:199388/139‑1 (MQ=255)
                                                                            aaTAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCca                                                   >  1:410899/1‑139 (MQ=255)
                                                                                   cagcaCAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATtata                                            >  1:1315422/1‑139 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       >  1:193389/1‑77 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGaaa                                                                                                       <  2:193389/77‑1 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                >  2:838420/1‑100 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTAt                                                                                <  1:838420/100‑1 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                                                   >  1:1449893/1‑113 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTg                                                                   <  2:1449893/113‑1 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCt                                         <  1:1023702/139‑1 (MQ=255)
                                                                                      cacaGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGAT‑TAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCt                                         >  2:1323016/1‑139 (MQ=255)
                                                                                            gTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtgt                                                                                    >  2:452528/1‑90 (MQ=255)
                                                                                            gTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGgtgt                                                                                    <  1:452528/90‑1 (MQ=255)
                                                                                                   tACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTccc                            >  1:979173/1‑139 (MQ=255)
                                                                                                     cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaaa                                                                                                                 >  1:510575/1‑51 (MQ=255)
                                                                                                     cAGCGCCAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATaaa                                                                                                                 <  2:510575/53‑3 (MQ=255)
                                                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     >  2:1494161/1‑107 (MQ=255)
                                                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  1:1494161/107‑1 (MQ=255)
                                                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                     >  2:1019216/1‑139 (MQ=255)
                                                                                                          ccAGACAGGATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGa                     >  1:1445902/1‑139 (MQ=255)
                                                                                                                  gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     <  2:1331718/99‑1 (MQ=255)
                                                                                                                  gATAGGAGAAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCAcc                                                     >  1:1331718/1‑99 (MQ=255)
                                                                                                                       gagaAAGAAAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCa        <  2:167872/139‑1 (MQ=255)
                                                                                                                        agaaagaaAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat       <  2:1315422/139‑1 (MQ=255)
                                                                                                                        agaaagaaAAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCat       <  2:314500/139‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                                                             <  1:1163632/81‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAAc                                                             >  2:1163632/1‑81 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCa                                                       >  2:228531/1‑87 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCa                                                       <  1:228531/87‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCATATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc   >  2:595158/1‑139 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                         <  1:1247997/117‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                         >  1:1550023/1‑117 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                         <  2:1550023/117‑1 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCaa                         >  2:1247997/1‑117 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc   >  2:680074/1‑139 (MQ=255)
                                                                                                                            agaaaAACATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGc   >  1:1261060/1‑139 (MQ=255)
                                                                                                                                 aaCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCt  <  1:1220285/135‑1 (MQ=255)
                                                                                                                                 aaCATAAGTTTATCAACGAGATA‑AATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCt  >  2:1220285/1‑135 (MQ=255)
                                                                                                                                       |                                                                                                                                 
AGAAAAAGTTATTAATGATCTGTTAAAAAACGATTATAAAATTTCCCTTTCTCCAGATACAACTTACCGAGAGTTGAATAAAGCAGCACAGCGTAGCATTACAGCGCCAGACAGGATAGGAGAAAGAAAAACATAGGTTTATCAACGAGATACAATGATTGAAAGAGGTGATAACAGCGGTGTTTATCAGTATGGTCGTGCTGAACATTTCACCCACATTATATCTGACAAACCTTCCCCAAAAGATAAATATGTTGCATATGCT  >  NC_000913/2786594‑2786858

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

GATK/CNVnator alignment

N/A