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 3,786,467 T→C I193V (ATT→GTT)  yibB ← protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,4670TC92.9% 76.9 / ‑0.5 28I193V (ATT→GTT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (18/8);  total (19/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

AATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGA  >  NZ_CP009273/3786394‑3786553
                                                                         |                                                                                      
aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                        >  1:13383/1‑90 (MQ=255)
aaTTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAAcc                                                                        >  1:44197/1‑90 (MQ=255)
        gAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGttgttg                                                                <  2:13383/90‑1 (MQ=255)
          aTTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTgag                                                              >  1:339930/1‑90 (MQ=255)
                    tttATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAAcc                                                                        <  2:674119/70‑1 (MQ=255)
                    tttATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAAcc                                                                        >  1:674119/1‑70 (MQ=255)
                     ttATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTa                                                   >  1:663692/1‑90 (MQ=255)
                         gCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGcc                                               <  2:94882/90‑1 (MQ=255)
                            gtgAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCtt                                            >  2:31665/1‑90 (MQ=255)
                                 ccAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCttttttt                                       >  2:259735/1‑90 (MQ=255)
                                   aacaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtat                                     <  2:83162/90‑1 (MQ=255)
                                   aacaaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTtat                                     <  2:107342/90‑1 (MQ=255)
                                     caaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAg                                   <  1:134756/90‑1 (MQ=255)
                                     caaTTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAg                                   <  2:184068/90‑1 (MQ=255)
                                              ccccAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGAt                          <  1:194298/90‑1 (MQ=255)
                                                         tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt               >  2:12681/1‑90 (MQ=255)
                                                         tataCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACtt               >  2:182482/1‑90 (MQ=255)
                                                           taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  1:355451/1‑90 (MQ=255)
                                                           taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:312119/1‑90 (MQ=255)
                                                           taCATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtc             >  2:677533/1‑90 (MQ=255)
                                                              aTGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatc          >  2:154904/1‑90 (MQ=255)
                                                               tGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTtcatca         >  1:522070/1‑90 (MQ=255)
                                                                gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGGTGCATCTTACTTTCATCaa        >  2:600797/1‑90 (MQ=255)
                                                                gATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCaa        >  1:239236/1‑90 (MQ=255)
                                                                  ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  1:268917/1‑90 (MQ=255)
                                                                  ttGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAt      >  1:13929/1‑90 (MQ=255)
                                                                    gCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATgg    >  2:130285/1‑90 (MQ=255)
                                                                      cAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGa  <  1:391469/90‑1 (MQ=255)
                                                                         |                                                                                      
AATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGA  >  NZ_CP009273/3786394‑3786553

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
GTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAAT  >  NZ_CP009273/3786368‑3786563
                                                                                                   |                                                                                                
GTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAAC                                                                                                  <  SRR3722112.661249/100‑1 (MQ=60)
                TTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACC                                                                                  >  SRR3722112.13553/1‑100 (MQ=60)
                TTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACC                                                                                  >  SRR3722112.44731/1‑100 (MQ=60)
                          AATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAG                                                                        >  SRR3722112.344121/1‑100 (MQ=60)
                           ccaccgagatctacacattgacattcgtcggcagcgtcagatgtgtataagagacagATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGA                                                                       <  SRR3722112.363761/43‑1 (MQ=60)
                                    ATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGctgtctctta                                                              >  SRR3722112.684656/1‑90 (MQ=60)
                                     TTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTA                                                             >  SRR3722112.674046/1‑100 (MQ=60)
                                                               CAATTGCGCCCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATG                                   <  SRR3722112.136434/100‑1 (MQ=60)
                                                                        CCCCAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACT                          <  SRR3722112.196672/100‑1 (MQ=60)
                                                                           CAATAATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTC                       >  SRR3722112.359971/1‑100 (MQ=60)
                                                                               AATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCA                   >  SRR3722112.529748/1‑100 (MQ=60)
                                                                                ATATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAA                  >  SRR3722112.242083/1‑100 (MQ=60)
                                                                                  ATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAT                >  SRR3722112.14109/1‑100 (MQ=60)
                                                                                  ATATACATGATTGCCAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAAT                >  SRR3722112.272070/1‑100 (MQ=60)
                                                                                                CAACCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAAT  <  SRR3722112.396593/100‑1 (MQ=60)
                                                                                                   |                                                                                                
GTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAACAATTGCGCCCCCAATAATATATACATGATTGCCAATCATGAAATCAAAAACCTGTTGTTGAGATGTGACAGTAAGCCCTTTTTTTATAGTGAAAAGATGCATCTTACTTTCATCAAATGGGAAATCCCAAAT  >  NZ_CP009273/3786368‑3786563

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: