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 F142 I1 R1 99 64.8 3148104 96.3% 3031624 99.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,436,937 A→T D305E (GAT→GAA)  pdxB ← erythronate‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,436,9370AT95.6% 235.4 / 1.5 68D305E (GAT→GAA) pdxBerythronate‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (31/34);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTAT  >  NC_000913/2436841‑2437031
                                                                                                |                                                                                              
aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             >  2:1254650/1‑100 (MQ=255)
aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             <  1:1190297/100‑1 (MQ=255)
aTCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCg                                                                                             >  2:1190290/1‑100 (MQ=255)
 tCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGa                                                                                            >  2:476909/1‑100 (MQ=255)
   gAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                          <  1:1393210/100‑1 (MQ=255)
   gAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGc                                                                                          >  1:862288/1‑100 (MQ=255)
     aCTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCgg                                                                                        >  1:269989/1‑100 (MQ=255)
      cTCACCCTGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCggg                                                                                       <  2:188764/100‑1 (MQ=255)
      cTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCggg                                                                                       <  1:355189/100‑1 (MQ=255)
      cTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCggg                                                                                       <  2:1521553/100‑1 (MQ=255)
       tCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                      >  1:143252/1‑100 (MQ=255)
        cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                     >  1:706866/1‑100 (MQ=255)
        cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                     >  2:1467483/1‑100 (MQ=255)
        cACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGcc                                                                                     >  2:405448/1‑100 (MQ=255)
           ccGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGc                                                                                      <  2:1485088/96‑1 (MQ=255)
            cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGTGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGTGCCATGc                                                                                 <  2:1084048/100‑1 (MQ=255)
            cGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGc                                                                                 >  1:904680/1‑100 (MQ=255)
                 tCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGt                                                                            <  2:1089108/100‑1 (MQ=255)
                 tCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGt                                                                            <  2:862281/100‑1 (MQ=255)
                  cccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTa                                                                           >  2:489551/1‑100 (MQ=255)
                   ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                          >  1:8729/1‑100 (MQ=255)
                   ccGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTaa                                                                          >  2:59119/1‑100 (MQ=255)
                         aCTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAc                                                                    <  2:1427706/100‑1 (MQ=255)
                         aCTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAc                                                                    <  1:1286865/100‑1 (MQ=255)
                          cTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAcc                                                                   <  2:269988/100‑1 (MQ=255)
                          cTTTACACAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGAcc                                                                   <  1:622843/100‑1 (MQ=255)
                             tACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCaaa                                                                >  2:1339563/1‑100 (MQ=255)
                               cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                              >  1:1353908/1‑100 (MQ=255)
                               cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                              >  2:1100419/1‑100 (MQ=255)
                               cGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAAct                                                              >  2:1315853/1‑100 (MQ=255)
                                   gCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTgg                                                          <  1:1443279/100‑1 (MQ=255)
                                   gCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTgg                                                          <  1:414874/100‑1 (MQ=255)
                                     ggTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcg                                                        <  1:1449012/100‑1 (MQ=255)
                                     ggTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcg                                                        >  1:1078597/1‑100 (MQ=255)
                                      gTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGgcgc                                                       >  2:1050535/1‑100 (MQ=255)
                                       tGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgca                                                      <  2:380732/100‑1 (MQ=255)
                                        gCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcag                                                     <  2:934790/100‑1 (MQ=255)
                                            cATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCgcaggcag                                                 >  2:364220/1‑100 (MQ=255)
                                               cgcgGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTaa                                              >  2:1226278/1‑100 (MQ=255)
                                                gcggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAAt                                             <  2:1498503/100‑1 (MQ=255)
                                                 cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  1:1254658/100‑1 (MQ=255)
                                                 cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  1:462361/100‑1 (MQ=255)
                                                 cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATg                                            <  1:1172463/100‑1 (MQ=255)
                                                 cggcgCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAAt                                             <  1:213067/99‑1 (MQ=255)
                                                   gcgcACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTa                                          <  1:1483770/100‑1 (MQ=255)
                                                         atcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgc                                    >  1:1221435/1‑100 (MQ=255)
                                                          tcatACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAgcg                                   >  1:551510/1‑100 (MQ=255)
                                                            atACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCa                                 >  1:1445925/1‑100 (MQ=255)
                                                               cacCAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACg                              >  1:638918/1‑100 (MQ=255)
                                                                  cAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                           >  1:840811/1‑100 (MQ=255)
                                                                  cAAATGCACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGc                           <  2:461200/100‑1 (MQ=255)
                                                                        cACCAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCa                     <  2:551507/100‑1 (MQ=255)
                                                                        cACCAGCCTTTACAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCa                     >  2:824519/1‑100 (MQ=255)
                                                                          ccAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATg                   >  1:63057/1‑100 (MQ=255)
                                                                          ccAGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATg                   <  2:930403/100‑1 (MQ=255)
                                                                            aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                 >  1:338467/1‑100 (MQ=255)
                                                                            aGCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGcc                 <  1:1226285/100‑1 (MQ=255)
                                                                             gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                <  2:889684/100‑1 (MQ=255)
                                                                             gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                <  2:1321873/100‑1 (MQ=255)
                                                                             gCCTTTTCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGccc                >  1:65539/1‑100 (MQ=255)
                                                                                 tttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATa            >  1:106367/1‑100 (MQ=255)
                                                                                 tttCAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATa            <  1:186193/100‑1 (MQ=255)
                                                                                    cAGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAAc         >  2:1161001/1‑100 (MQ=255)
                                                                                     aGCGTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACt        <  2:1161357/100‑1 (MQ=255)
                                                                                        gTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGc     <  2:1080016/100‑1 (MQ=255)
                                                                                        gTCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGc     <  2:904672/100‑1 (MQ=255)
                                                                                         tCGGTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCt    <  1:59119/100‑1 (MQ=255)
                                                                                           ggTTGTTCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCtat  >  2:801957/1‑100 (MQ=255)
                                                                                                |                                                                                              
ATCGAACTCACCCGGTATCCCGGCGACTTTACGCAGCGGTGCGTCATCGCGGCGCACATCATACACCAAATGCACCAGCCTTTTCAGCGTCGGTTGATCGAGCGGGCCATGCAGCGTAATGCGACCAAACTCTGGCGCAGGCAGTAATGTATCCAGCGCAACGTGCTGTTCATGCCCAATAAACTTGCTAT  >  NC_000913/2436841‑2437031

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

GATK/CNVnator alignment

N/A