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 F56 I1 R1 155 65.2 3164312 95.9% 3034575 99.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,910,662 A→T D288E (GAT→GAA)  pstS ← phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,910,6620AT97.1% 233.8 / ‑1.4 69D288E (GAT→GAA) pstSphosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (24/43);  total (26/43)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910570‑3910758
                                                                                            |                                                                                                
gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           <  2:1289852/100‑1 (MQ=255)
gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           <  2:183536/100‑1 (MQ=255)
gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           <  1:1425641/100‑1 (MQ=255)
 cGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgtg                                                                                          <  2:1553631/100‑1 (MQ=255)
   tAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGa                                                                                        <  2:968191/100‑1 (MQ=255)
     aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                      <  1:839081/100‑1 (MQ=255)
      tCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCa                                                                                     >  1:285672/1‑100 (MQ=255)
       ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    >  2:686904/1‑100 (MQ=255)
       ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    <  2:565930/100‑1 (MQ=255)
       ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    <  2:461769/100‑1 (MQ=255)
       ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    <  2:31463/100‑1 (MQ=255)
        cAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGa                                                                                   >  2:1197041/1‑100 (MQ=255)
          ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 <  1:1197046/100‑1 (MQ=255)
          ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 >  1:404762/1‑100 (MQ=255)
           gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATg                                                                                <  2:506587/100‑1 (MQ=255)
             cGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGaa                                                                              >  1:1246200/1‑100 (MQ=255)
             cGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGaa                                                                              >  2:1270793/1‑100 (MQ=255)
               ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  2:1331595/100‑1 (MQ=255)
               ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  1:456414/100‑1 (MQ=255)
               ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  2:844187/100‑1 (MQ=255)
               ctCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  1:804625/99‑1 (MQ=255)
                tCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACga                                                                           >  2:804620/1‑99 (MQ=255)
                tCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGt                                                                           >  2:1077486/1‑100 (MQ=255)
                 cGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTg                                                                          <  2:1187648/100‑1 (MQ=255)
                   ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                        <  1:1022215/100‑1 (MQ=255)
                      gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                     <  1:138974/100‑1 (MQ=255)
                         tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                  <  2:337554/100‑1 (MQ=255)
                            ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                               >  2:1209921/1‑100 (MQ=255)
                             cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              >  1:1546505/1‑100 (MQ=255)
                             cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              <  2:212507/100‑1 (MQ=255)
                               cGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                            <  1:12712/100‑1 (MQ=255)
                                 gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGcc                                                          <  1:1386925/100‑1 (MQ=255)
                                   tttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCAt                                                        <  1:447888/100‑1 (MQ=255)
                                      gTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                     <  2:1325063/100‑1 (MQ=255)
                                             cAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTACTAGGCCATGCATCTTCGc                                              >  2:667885/1‑100 (MQ=255)
                                             cAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGc                                              <  1:238664/100‑1 (MQ=255)
                                              aGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGcc                                             <  2:789413/100‑1 (MQ=255)
                                                  gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  1:621461/1‑100 (MQ=255)
                                                  gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  1:1389606/1‑100 (MQ=255)
                                                   aagtaTTTCAGCACTTCTGTGCGTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCttt                                          <  2:1339082/98‑1 (MQ=255)
                                                   aagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCt                                        >  2:1087297/1‑100 (MQ=255)
                                                   aagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCt                                        >  2:816008/1‑100 (MQ=255)
                                                      aaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggt                                     >  2:1103378/1‑100 (MQ=255)
                                                         ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                  <  2:464122/100‑1 (MQ=255)
                                                         ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAACAGGCCATGCATCTTCGCCTTTCtggttgg                                  <  1:21355/100‑1 (MQ=255)
                                                           cAGCACTTCTGTGCTTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAGCGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                <  2:1456173/100‑1 (MQ=255)
                                                           cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                >  2:1265999/1‑100 (MQ=255)
                                                           cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                >  2:252211/1‑100 (MQ=255)
                                                            aGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCa                               >  1:517657/1‑100 (MQ=255)
                                                             gCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCACCTTCGCCTTTCTGGTTGGTCAg                              <  2:918779/100‑1 (MQ=255)
                                                                cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGTATCTTCGCCTTTCTGGTTGGTCAGATc                           <  2:726157/100‑1 (MQ=255)
                                                                cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATc                           <  1:564109/100‑1 (MQ=255)
                                                                 ttCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATcc                          <  1:43970/100‑1 (MQ=255)
                                                                  tCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCt                         >  1:171004/1‑100 (MQ=255)
                                                                   cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        >  2:277967/1‑100 (MQ=255)
                                                                   cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        >  2:161974/1‑100 (MQ=255)
                                                                     gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      <  1:1266004/100‑1 (MQ=255)
                                                                      tgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGc                     <  1:686906/100‑1 (MQ=255)
                                                                       gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                    <  2:464448/100‑1 (MQ=255)
                                                                       gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                    <  1:63058/100‑1 (MQ=255)
                                                                              tCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGttt             <  1:1270798/100‑1 (MQ=255)
                                                                               cTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtttt            <  2:1202940/100‑1 (MQ=255)
                                                                                tagTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTg           >  1:1069641/3‑100 (MQ=255)
                                                                                 ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  1:1571824/100‑1 (MQ=255)
                                                                                 ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  1:816013/100‑1 (MQ=255)
                                                                                        tctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc   <  2:1246195/100‑1 (MQ=255)
                                                                                        tctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGCTCAGATCCTGAGCGAAGGTTTTGCTCCAGt    >  1:1330916/1‑99 (MQ=255)
                                                                                         ctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCt  >  2:425511/1‑100 (MQ=255)
                                                                                         ctGTTCTCTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCt  >  1:766133/1‑100 (MQ=255)
                                                                                            |                                                                                                
GCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910570‑3910758

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

GATK/CNVnator alignment

N/A