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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,076,727 (T)5→4 intergenic (+458/‑85) putP → / → BW25113_RS05295 sodium/proline symporter PutP/FTR1 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,076,7230T.96.5% 109.6 / ‑2.9 29intergenic (+454/‑89)putP/BW25113_RS05295sodium/proline symporter PutP/FTR1 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (15/13);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTACCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAGGGTGGGT  >  NZ_CP009273/1076635‑1076808
                                                                                        |                                                                                     
tttCTACCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGtt                                                                                      >  2:321869/1‑90 (MQ=255)
  tCTACCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGtttt                                                                                    >  2:592372/1‑90 (MQ=255)
      cccTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGttttcttt                                                                                <  1:801561/90‑5 (MQ=255)
       ccTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTg                                                                              >  2:450461/1‑90 (MQ=255)
        cTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACAGTCCGCTTGTGTTTTCTTTGc                                                                             >  2:896736/1‑90 (MQ=255)
             ataatCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTaa                                                                        <  2:129430/90‑1 (MQ=255)
             ataatCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTAAAAACATTCCGCTTGTGTTTTCTTTGCCGTaa                                                                        <  1:528142/90‑1 (MQ=255)
               aatCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATg                                                                      >  2:120698/1‑90 (MQ=255)
               aatCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATg                                                                      >  2:254282/1‑90 (MQ=255)
                    tttATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAAt                                                                 >  1:406056/1‑90 (MQ=255)
                           cccTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATc                                                          <  1:418424/90‑1 (MQ=255)
                             cTTGGCATTAGCTCTCTTTGTTTACATTCCAACTTCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATTAc                                                        >  2:7347/1‑90 (MQ=255)
                             cTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCAc                                                        >  1:595399/1‑90 (MQ=255)
                              ttGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACt                                                       >  2:118324/1‑90 (MQ=255)
                              ttGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACt                                                       <  1:268193/90‑1 (MQ=255)
                              ttGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGAGAATCGCTATCACt                                                       >  1:718081/1‑90 (MQ=255)
                               tGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTg                                                      <  1:370065/90‑1 (MQ=255)
                                     tACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGAttt                                                <  2:678091/90‑1 (MQ=255)
                                                 tttACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACtt                                            >  1:764782/1‑82 (MQ=255)
                                                 tttACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACtt                                            <  2:764782/82‑1 (MQ=255)
                                                      aTTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCAt                                  >  1:611303/1‑87 (MQ=255)
                                                      aTTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCAt                                  <  2:611303/87‑1 (MQ=255)
                                                          cAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGa                           >  2:156717/1‑90 (MQ=255)
                                                            aCATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACt                         <  1:450461/90‑1 (MQ=255)
                                                                 atTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGa                                                   <  1:884214/59‑1 (MQ=255)
                                                                 atTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGa                                                   >  2:884214/1‑59 (MQ=255)
                                                                 atTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGaa                    >  2:859764/1‑90 (MQ=255)
                                                                    ttATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGaaaaa                 <  2:344526/90‑1 (MQ=255)
                                                                            aTTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAg         >  1:707343/1‑90 (MQ=255)
                                                                               ccGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAgggt      <  1:441503/90‑1 (MQ=255)
                                                                               ccGCTTGTTTTTTCTTTGACCTAATGATAATCGCTATCACTGCGATTTGCTTTTTTTTGCATAGATTGACTCAGAAAAACGTTTAAgggt      <  2:326577/90‑1 (MQ=255)
                                                                                   ttGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAgggtgggt  <  1:838824/90‑1 (MQ=255)
                                                                                        |                                                                                     
TTTCTACCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAGGGTGGGT  >  NZ_CP009273/1076635‑1076808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 30 ≤ 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
TCCCATCTTTCTACCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAGGGTGGGTGGCATGTTTGTTC  >  NZ_CP009273/1076628‑1076821
                                                                                                   |                                                                                              
acgtctcgtgggctcggagatgtgtataagagacagCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTC                                                                                                <  SRR3722072.431051/64‑1 (MQ=60)
             CCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGA                                                                                  <  SRR3722072.815563/100‑1 (MQ=60)
                 GGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAAT                                                                              >  SRR3722072.413390/1‑100 (MQ=60)
                    ATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTAAAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGC                                                                           <  SRR3722072.537655/100‑1 (MQ=60)
                          GTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCAC                                                                     >  SRR3722072.606018/1‑100 (MQ=60)
                           TTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGAGAATCGCTATCACT                                                                    >  SRR3722072.730859/1‑100 (MQ=60)
                                  CCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTT                                                             <  SRR3722072.426010/100‑1 (MQ=60)
                                     TTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACT                                                          <  SRR3722072.272756/100‑1 (MQ=60)
                                      TGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTT                                                         <  SRR3722072.376641/100‑1 (MQ=60)
                                         ggagatgtgtataagagacagTTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTC                                                      <  SRR3722072.899434/79‑1 (MQ=60)
                                              CCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGC                                                 >  SRR3722072.778260/1‑100 (MQ=60)
                                                   CTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGA                                            >  SRR3722072.622202/1‑100 (MQ=60)
                                                                   ACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACG                            <  SRR3722072.458709/100‑1 (MQ=60)
                                                                         TTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAG                      >  SRR3722072.719955/1‑100 (MQ=60)
                                                                                      CCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAGGGTGGGTGGCATG         <  SRR3722072.449575/100‑1 (MQ=60)
                                                                                          TTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAGGGTGGGTGGCATGTTTG     <  SRR3722072.853360/100‑1 (MQ=60)
                                                                                              tgTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAGGGTGGGTGGCATGTTTGTTC  <  SRR3722072.466902/98‑1 (MQ=60)
                                                                                                   |                                                                                              
TCCCATCTTTCTACCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAACGTTTAAGGGTGGGTGGCATGTTTGTTC  >  NZ_CP009273/1076628‑1076821

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