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
A1 F2 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,671,899 T→G T88P (ACG→CCG)  pntA ← Re/Si‑specific NAD(P)(+) transhydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,671,8990TG90.0% 86.7 / 2.7 30T88P (ACG→CCG) pntARe/Si‑specific NAD(P)(+) transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (13/14);  total (14/16)
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.36e-01

CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC  >  NZ_CP009273/1671816‑1671981
                                                                                   |                                                                                  
cATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTccc                                                                              >  2:672183/1‑90 (MQ=255)
       aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAGTTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCACAGGATTc                                                                       >  2:198096/1‑90 (MQ=255)
       aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGTGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTc                                                                       <  1:758312/90‑1 (MQ=255)
       aTCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTc                                                                       >  1:28665/1‑90 (MQ=255)
                 cgttacgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCaa                                                             <  2:66373/90‑1 (MQ=255)
                   ttacgttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAAtt                                                           >  1:314509/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTtcat                                                       <  1:437078/90‑1 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                       <  2:611882/90‑1 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                       >  2:14790/1‑90 (MQ=255)
                       gttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcat                                                       <  2:13108/90‑1 (MQ=255)
                        ttCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTtcatc                                                      <  1:14790/90‑1 (MQ=255)
                           cGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCatcatc                                                   >  2:488264/1‑90 (MQ=255)
                                 tttttGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACgg                                             <  1:729047/90‑1 (MQ=255)
                                 tttttGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACgg                                             <  1:443148/90‑1 (MQ=255)
                                       cATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCAtt                                       <  1:138315/90‑1 (MQ=255)
                                        aTTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTg                                      <  1:89974/90‑1 (MQ=255)
                                         ttAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGa                                     >  2:570175/1‑90 (MQ=255)
                                          tAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAc                                    >  1:573631/1‑90 (MQ=255)
                                           aaTTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAcc                                   >  2:814096/1‑90 (MQ=255)
                                             ttCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCtt                                 <  2:60170/90‑1 (MQ=255)
                                                      cTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAt                        >  1:755313/1‑90 (MQ=255)
                                                      cTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAt                        >  2:780980/1‑90 (MQ=255)
                                                             ggCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc                 <  1:563366/90‑1 (MQ=255)
                                                             ggCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGAc                 <  2:755313/90‑1 (MQ=255)
                                                                  gagaaaaCTCACCAGCGGTGTCCCAGGATTCAGTACCGCAATTTCATCATCTACCGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa            <  2:522458/87‑1 (MQ=255)
                                                                  gATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCa            <  2:682072/90‑1 (MQ=255)
                                                                   aTAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAg           <  1:420090/90‑1 (MQ=255)
                                                                       aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc       >  2:556363/1‑90 (MQ=255)
                                                                       aaCTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGc       >  1:731588/1‑90 (MQ=255)
                                                                            aCCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTc  >  1:297183/1‑90 (MQ=255)
                                                                                   |                                                                                  
CATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC  >  NZ_CP009273/1671816‑1671981

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTTCTACAATTTCAGC  >  NZ_CP009273/1671801‑1671998
                                                                                                  |                                                                                                   
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGT                                                                                                    >  SRR3722079.254653/1‑100 (MQ=60)
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGT                                                                                                    >  SRR3722079.354428/1‑100 (MQ=60)
         AGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGA                                                                                           <  SRR3722079.656550/100‑1 (MQ=60)
            GTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTC                                                                                        >  SRR3722079.28878/1‑100 (MQ=60)
                      ATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGTGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAA                                                                              <  SRR3722079.766425/100‑1 (MQ=60)
                        CACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATT                                                                            >  SRR3722079.317050/1‑100 (MQ=60)
                                      GTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGG                                                              <  SRR3722079.440871/100‑1 (MQ=60)
                                       TTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGC                                                             <  SRR3722079.14896/100‑1 (MQ=60)
                                               GTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGAC                                                     >  SRR3722079.579282/1‑100 (MQ=60)
                                                TTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACC                                                    <  SRR3722079.447014/100‑1 (MQ=60)
                                                TTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACC                                                    <  SRR3722079.736839/100‑1 (MQ=60)
                                                      CATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGA                                              <  SRR3722079.139486/100‑1 (MQ=60)
                                                       ATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAA                                             <  SRR3722079.90686/100‑1 (MQ=60)
                                                           ATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGAT                                         >  SRR3722079.763388/1‑100 (MQ=60)
                                                                            GGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGC                        <  SRR3722079.568852/100‑1 (MQ=60)
                                                                            GGCCAGATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGC                        >  SRR3722079.739417/1‑100 (MQ=60)
                                                                                 GATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTC                   >  SRR3722079.299619/1‑100 (MQ=60)
                                                                                  ATAAAACTCACCAGCGGTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCC                  <  SRR3722079.423708/100‑1 (MQ=60)
                                                                                                  GTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTTCTACAATTTCAGC  >  SRR3722079.504313/1‑100 (MQ=60)
                                                                                                  |                                                                                                   
ACGCGGCACAGAGTCCATCGCCATCACGGTCACGTTACGTTCCGCAAGTTTTTGCATTAATTCCGGATTCTGCGCAGGCCAGATAAAACTCACCAGCGTTGTCCCAGGATTCAGTAACGCAATTTCATCATCTAACGGCGCATTGACCTTCAGAATGATCTCTGACTGCCAGACGCTATTCCCTTCTACAATTTCAGC  >  NZ_CP009273/1671801‑1671998

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: