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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,494,539 T→C F222S (TTC→TCC)  nanX → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,494,5390TC96.2% 77.6 / ‑3.5 26F222S (TTC→TCC) nanXMFS transporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (10/15);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

ATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG  >  NZ_CP009273/4494471‑4494625
                                                                    |                                                                                      
atatAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTAtt                                                                   <  1:203124/90‑1 (MQ=255)
          aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                         >  2:645218/1‑90 (MQ=255)
           aaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACt                                                        <  2:199126/90‑1 (MQ=255)
                      ttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAc                                             <  2:684944/90‑1 (MQ=255)
                      ttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAc                                             <  2:140130/90‑1 (MQ=255)
                      ttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAc                                             <  1:623558/90‑1 (MQ=255)
                          ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                         >  1:194382/1‑90 (MQ=255)
                              ctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGactact                                     <  2:578598/90‑1 (MQ=255)
                              ctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGactact                                     >  1:627913/1‑90 (MQ=255)
                               ttCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                    <  2:557761/90‑1 (MQ=255)
                                   gAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActtcctt                                <  1:537822/90‑1 (MQ=255)
                                      aaCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTtccttcct                             >  1:656760/1‑90 (MQ=255)
                                       aCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTa                            <  1:710837/90‑1 (MQ=255)
                                          acaTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCt                         >  1:114065/1‑90 (MQ=255)
                                                     tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTTCTTCCTACCTGGCAGATAATg              <  2:557174/90‑1 (MQ=255)
                                                     tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  1:319902/90‑1 (MQ=255)
                                                     tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg              <  2:166143/90‑1 (MQ=255)
                                                       tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa            <  1:673292/90‑1 (MQ=255)
                                                        ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg           >  1:79122/1‑90 (MQ=255)
                                                          aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt         <  1:302209/90‑1 (MQ=255)
                                                              tCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATa     >  1:506748/1‑90 (MQ=255)
                                                              tCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACgg                                           >  2:419694/1‑52 (MQ=255)
                                                              tCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACgg                                           <  1:419694/52‑1 (MQ=255)
                                                              tCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATa     >  1:229489/1‑90 (MQ=255)
                                                                gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  2:342490/1‑90 (MQ=255)
                                                                gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa   >  1:5968/1‑90 (MQ=255)
                                                                 tttcCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAg  <  1:422420/90‑1 (MQ=255)
                                                                    |                                                                                      
ATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG  >  NZ_CP009273/4494471‑4494625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
CTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  >  NZ_CP009273/4494440‑4494635
                                                                                                   |                                                                                                
CTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTC                                                                                                  >  SRR3722112.487236/1‑100 (MQ=60)
                               ATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAC                                                                   <  SRR3722112.205592/100‑1 (MQ=60)
                                               ACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGAC                                                   >  SRR3722112.196758/1‑100 (MQ=60)
                                                   TTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACT                                               >  SRR3722112.637613/1‑100 (MQ=60)
                                                     TTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTC                                             <  SRR3722112.633190/100‑1 (MQ=60)
                                                       ggcagcgtcagatgtgtataagagacagATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCT                                           <  SRR3722112.425254/72‑1 (MQ=60)
                                                           GTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCT                                       >  SRR3722112.667008/1‑100 (MQ=60)
                                                               TCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCT                                   >  SRR3722112.115480/1‑100 (MQ=60)
                                                                  GAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGC                                <  SRR3722112.545846/100‑1 (MQ=60)
                                                                      ACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGAT                            <  SRR3722112.721951/100‑1 (MQ=60)
                                                                             CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                     >  SRR3722112.80086/1‑100 (MQ=60)
                                                                                   ATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATA               >  SRR3722112.514126/1‑100 (MQ=60)
                                                                                   ATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATA               >  SRR3722112.232247/1‑100 (MQ=60)
                                                                                    TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC              <  SRR3722112.323642/100‑1 (MQ=60)
                                                                                     CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA             >  SRR3722112.6043/1‑100 (MQ=60)
                                                                                      TCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG            <  SRR3722112.683804/100‑1 (MQ=60)
                                                                                       CTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGT           >  SRR3722112.514769/1‑100 (MQ=60)
                                                                                       CTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGT           >  SRR3722112.611016/1‑100 (MQ=60)
                                                                                         ATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACctgtc         >  SRR3722112.674667/1‑95 (MQ=60)
                                                                                         ATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGG         <  SRR3722112.305687/100‑1 (MQ=60)
                                                                                                TTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  <  SRR3722112.428042/100‑1 (MQ=60)
                                                                                                   |                                                                                                
CTCCAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  >  NZ_CP009273/4494440‑4494635

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