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 2,102,646 T→C G140G (GGA→GGG wzx ← O16 family O‑antigen flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,102,6460TC96.8% 95.4 / ‑3.4 31G140G (GGA→GGGwzxO16 family O‑antigen flippase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/16);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTA  >  NZ_CP009273/2102561‑2102727
                                                                                     |                                                                                 
ttGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCt                                                                               >  1:521684/1‑90 (MQ=255)
  gtgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGa                                                                             >  1:163874/1‑90 (MQ=255)
   tgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGaa                                                                            <  1:190059/90‑1 (MQ=255)
     ttCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAAc                                                                          <  1:712829/90‑1 (MQ=255)
     ttCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAAc                                                                          <  2:163874/90‑1 (MQ=255)
        acaAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                       <  2:207663/90‑1 (MQ=255)
                 aatGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAggg                                                              >  1:333607/1‑90 (MQ=255)
                  atGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGt                                                             <  2:33782/90‑1 (MQ=255)
                         agGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTaa                                                      <  1:722686/90‑1 (MQ=255)
                            gATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAAtt                                                   >  1:173025/1‑90 (MQ=255)
                                    ccAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTa                                           >  2:124492/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  1:446291/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  1:449490/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  2:292914/1‑90 (MQ=255)
                                     cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                          >  2:549280/1‑90 (MQ=255)
                                      aaGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAAt                                         >  1:63812/1‑90 (MQ=255)
                                        gCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGc                                       <  1:331312/90‑1 (MQ=255)
                                                  aTACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAAc                             <  2:278702/90‑1 (MQ=255)
                                                   tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                            <  1:1351/90‑1 (MQ=255)
                                                   tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                            <  2:417156/90‑1 (MQ=255)
                                                   tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                            <  1:124492/90‑1 (MQ=255)
                                                       gctTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCa                        >  1:172799/2‑90 (MQ=255)
                                                                 aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg              >  1:457287/1‑90 (MQ=255)
                                                                 aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg              >  2:365809/1‑90 (MQ=255)
                                                                  gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg             >  1:540426/1‑90 (MQ=255)
                                                                  gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg             <  2:449490/90‑1 (MQ=255)
                                                                  gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg             <  2:691473/90‑1 (MQ=255)
                                                                       ttcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGc        <  2:716271/90‑1 (MQ=255)
                                                                        tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt       <  1:317182/90‑1 (MQ=255)
                                                                        tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt       <  1:176428/90‑1 (MQ=255)
                                                                          atttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTg     >  2:267976/1‑90 (MQ=255)
                                                                             ttttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTta  <  2:157064/90‑1 (MQ=255)
                                                                             ttttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTta  <  1:479469/90‑1 (MQ=255)
                                                                                     |                                                                                 
TTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTA  >  NZ_CP009273/2102561‑2102727

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAA  >  NZ_CP009273/2102551‑2102737
                                                                                               |                                                                                           
ATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCT                                                                                         >  SRR3722112.529355/1‑100 (MQ=60)
  ATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGA                                                                                       >  SRR3722112.165893/1‑100 (MQ=60)
             TGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGAT                                                                            <  SRR3722112.192389/100‑1 (MQ=60)
               TTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAG                                                                          <  SRR3722112.723970/100‑1 (MQ=60)
                 CACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGG                                                                        >  SRR3722112.337655/1‑100 (MQ=60)
                            ATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATT                                                             >  SRR3722112.175149/1‑100 (MQ=60)
                                   AGGGATAATAGCCAAGCGGGATAAAATACTCCTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACT                                                      <  SRR3722112.733966/100‑1 (MQ=60)
                                     GGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAA                                                    >  SRR3722112.452502/1‑100 (MQ=60)
                                     GGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAA                                                    >  SRR3722112.455762/1‑100 (MQ=60)
                                      GATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAAT                                                   >  SRR3722112.64584/1‑100 (MQ=60)
                                                  GCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAAC                                       <  SRR3722112.335310/100‑1 (MQ=60)
                                                       ATAAAATACTGCTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCA                                  >  SRR3722112.174923/1‑100 (MQ=60)
                                                             TACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAAT                            <  SRR3722112.1367/100‑1 (MQ=60)
                                                             TACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAAT                            <  SRR3722112.126032/100‑1 (MQ=60)
                                                                 ACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACG                        >  SRR3722112.463723/1‑100 (MQ=60)
                                                                  CTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGG                       >  SRR3722112.548507/1‑100 (MQ=60)
                                                                                  TCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC       <  SRR3722112.320863/100‑1 (MQ=60)
                                                                                  TCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC       <  SRR3722112.178592/100‑1 (MQ=60)
                                                                                   CATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCA      >  SRR3722112.401735/1‑100 (MQ=60)
                                                                                       TTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAA  <  SRR3722112.486338/100‑1 (MQ=60)
                                                                                               |                                                                                           
ATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAA  >  NZ_CP009273/2102551‑2102737

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: