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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,064,564 A→G Y185H (TAT→CAT)  yihS ← sulfoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,064,5640AG91.5% 124.3 / 4.8 47Y185H (TAT→CAT) yihSsulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (23/20);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCG  >  NZ_CP009273/4064479‑4064640
                                                                                     |                                                                            
tGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACa                                                                          <  2:425670/90‑1 (MQ=255)
  attattTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAAt                                                                        <  1:316133/90‑1 (MQ=255)
       ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                   <  1:944811/90‑1 (MQ=255)
       ttCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGa                                                                   <  1:467872/90‑1 (MQ=255)
        tCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaa                                                                  >  2:1191731/1‑90 (MQ=255)
        tCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaa                                                                  >  1:275843/1‑90 (MQ=255)
         cTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaaa                                                                 >  2:1136206/1‑90 (MQ=255)
         cTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGaaa                                                                 <  2:961630/90‑1 (MQ=255)
           ggCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGc                                                               >  2:28511/1‑90 (MQ=255)
                  tcgtGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAcc                                                        <  2:1289626/90‑1 (MQ=255)
                     tGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACcgcg                                                     <  1:1564664/90‑1 (MQ=255)
                      ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                    <  2:1478271/90‑1 (MQ=255)
                      ggATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGt                                                    <  1:994626/90‑1 (MQ=255)
                        aTAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGc                                                  <  2:684930/90‑1 (MQ=255)
                           aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                               >  2:266916/1‑90 (MQ=255)
                           aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                               >  2:532210/1‑90 (MQ=255)
                           aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCata                                               >  1:1388989/1‑90 (MQ=255)
                           aTCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCATGCata                                               >  2:938108/1‑90 (MQ=255)
                              aCGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTg                                            >  1:85163/1‑90 (MQ=255)
                                ggAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGc                                          >  1:123411/1‑90 (MQ=255)
                                ggAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGc                                          >  2:1290283/1‑90 (MQ=255)
                                  aaGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCAt                                        >  1:1502671/1‑90 (MQ=255)
                                     ccACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGCGCATATTGGCATTgc                                     >  2:1113170/1‑90 (MQ=255)
                                      cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                    <  1:532210/90‑1 (MQ=255)
                                      cACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgcc                                    >  1:501830/1‑90 (MQ=255)
                                        cgcgAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgccgc                                  <  1:91285/90‑1 (MQ=255)
                                           gAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACa                                                                          <  2:684564/47‑1 (MQ=255)
                                           gAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACa                                                                          >  1:684564/1‑47 (MQ=255)
                                           gAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTgccgcccc                               <  2:201048/90‑3 (MQ=255)
                                           gAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCcgccgc                               >  2:1106745/1‑90 (MQ=255)
                                           gAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCcgccgc                               >  2:26689/1‑90 (MQ=255)
                                           gAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCcgccgc                               <  2:411594/90‑1 (MQ=255)
                                             aTCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCgg                             >  1:830729/1‑90 (MQ=255)
                                              tCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGt                            <  2:916263/90‑1 (MQ=255)
                                              tCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGt                            >  2:24539/1‑90 (MQ=255)
                                              tCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGt                            <  2:1388989/90‑1 (MQ=255)
                                                               tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg           >  2:1320182/1‑90 (MQ=255)
                                                               tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg           >  1:1063935/1‑90 (MQ=255)
                                                               tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg           >  1:1485174/1‑90 (MQ=255)
                                                               tttttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTg           >  1:300035/1‑90 (MQ=255)
                                                                 tttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGaa         <  2:1502671/90‑1 (MQ=255)
                                                                  ttGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAg        <  2:241800/90‑1 (MQ=255)
                                                                      cGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCtt    <  1:1290283/90‑1 (MQ=255)
                                                                       gtgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTc   <  1:53/90‑1 (MQ=255)
                                                                        tgagtgaCGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCg  >  2:874913/1‑90 (MQ=255)
                                                                         gagtgaCGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCt     >  2:1354749/1‑86 (MQ=255)
                                                                         gagtgaCGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCt     <  1:1354749/86‑1 (MQ=255)
                                                                                     |                                                                            
TGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCG  >  NZ_CP009273/4064479‑4064640

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
TTAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC  >  NZ_CP009273/4064467‑4064658
                                                                                                 |                                                                                              
TTAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAA                                                                                              <  SRR3722078.439543/100‑1 (MQ=60)
          AATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAA                                                                                    >  SRR3722078.279207/1‑100 (MQ=60)
              ATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCT                                                                                <  SRR3722078.320001/100‑1 (MQ=60)
                   TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722078.958005/100‑1 (MQ=60)
                   TTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCAC                                                                           <  SRR3722078.473444/100‑1 (MQ=60)
                             GTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATA                                                                 >  SRR3722078.1411427/1‑100 (MQ=60)
                                GTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTG                                                              >  SRR3722078.86164/1‑100 (MQ=60)
                                 TGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGG                                                             <  SRR3722078.1590133/100‑1 (MQ=60)
                                  GGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGC                                                            <  SRR3722078.1009096/100‑1 (MQ=60)
                                  GGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGC                                                            >  SRR3722078.124847/1‑100 (MQ=60)
                                    ATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCAT                                                          >  SRR3722078.1527178/1‑100 (MQ=60)
                                        TCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCC                                                      >  SRR3722078.507770/1‑100 (MQ=60)
                                             GAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTgtctcttatacacatctgacgctgccgacga                                                 >  SRR3722078.692578/1‑69 (MQ=60)
                                               AGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGG                                               >  SRR3722078.841767/1‑100 (MQ=60)
                                                  CACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTAC                                            <  SRR3722078.538502/100‑1 (MQ=60)
                                                    CGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTC                                          <  SRR3722078.92356/100‑1 (MQ=60)
                                                                 ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             >  SRR3722078.303702/1‑100 (MQ=60)
                                                                 ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             >  SRR3722078.1509387/1‑100 (MQ=60)
                                                                 ATCCAGCCATTTTTTGTCGTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTG                             >  SRR3722078.1079960/1‑100 (MQ=60)
                                                                                 TCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGG             <  SRR3722078.1376495/100‑1 (MQ=60)
                                                                                  CGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGA            <  SRR3722078.1310821/100‑1 (MQ=60)
                                                                                   GTGAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGAT           <  SRR3722078.54/100‑1 (MQ=60)
                                                                                     GAGTGACGTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTC         >  SRR3722078.1146820/1‑100 (MQ=60)
                                                                                            GTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC  >  SRR3722078.42767/1‑100 (MQ=60)
                                                                                            GTCATGAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC  >  SRR3722078.1532224/1‑100 (MQ=60)
                                                                                                 |                                                                                              
TTAACGCGATAATGATTATTTCTGGCGACGTCGTGGATAATCACGGAAGCCACGCGAATCGCGCGATCCAGCCATTTTTTGTCGTGAGTGACGTCATAAACAATCAAGAAAGCTTCCACCGCGTGCATATTGGCATTGCCGCCGCGGTACTCTTCGGTTTTGCTGAAGGCTTCGTCCCAGGATTCCAGGCAC  >  NZ_CP009273/4064467‑4064658

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