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 2,670,359 T→C G454G (GGA→GGG yphE ← sugar ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,670,3590TC97.9% 148.7 / ‑3.7 47G454G (GGA→GGGyphEsugar ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (31/15);  total (32/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAAC  >  NZ_CP009273/2670283‑2670443
                                                                            |                                                                                    
cgtgcTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         <  1:299223/90‑1 (MQ=255)
  tgcTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagccagc                                                                       <  2:171245/90‑1 (MQ=255)
          taaCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGAc                                                               <  1:191910/90‑1 (MQ=255)
           aaCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACa                                                              >  2:488513/1‑90 (MQ=255)
            aCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACaa                                                             <  2:1338812/90‑1 (MQ=255)
             cAGGATGCGGTCACACACCAGCGGCGACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt                                                            <  1:1424061/90‑1 (MQ=255)
              aGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATa                                                           >  1:9731/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:200360/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         >  1:138703/1‑74 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         <  2:138703/74‑1 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:657214/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         <  2:339425/90‑1 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:1206538/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:1027500/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:1087651/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:903828/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:830694/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:560811/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:430071/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:138138/1‑90 (MQ=255)
                  tGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGt                                                       <  1:657214/90‑1 (MQ=255)
                       tCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATc                                                  <  2:16750/90‑1 (MQ=255)
                         acacacCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATctg                                                >  1:1305914/1‑90 (MQ=255)
                              ccAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCCCGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGtt                                           <  1:1328240/90‑1 (MQ=255)
                              ccAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGtt                                           >  2:333988/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt                                                            <  1:200016/71‑1 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt                                                            >  2:200016/1‑71 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  1:425710/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  1:1352262/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  2:75466/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  1:48643/1‑90 (MQ=255)
                                    gCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCtt                                     <  2:1244808/90‑1 (MQ=255)
                                    gCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCtt                                     >  1:521195/1‑90 (MQ=255)
                                      aaCTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCg                                   >  2:918350/1‑90 (MQ=255)
                                          cctccACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGatat                               >  2:1213033/1‑90 (MQ=255)
                                               aCTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         >  1:287290/1‑43 (MQ=255)
                                               aCTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         <  2:287290/43‑1 (MQ=255)
                                                  tCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGAcgccg                       <  2:152946/90‑1 (MQ=255)
                                                          gaTAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTgg               >  2:745344/1‑90 (MQ=255)
                                                               acacCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgt          >  1:1219895/1‑90 (MQ=255)
                                                                 accacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgtcg        >  2:1445683/1‑90 (MQ=255)
                                                                 accacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgtcg        >  2:626532/1‑90 (MQ=255)
                                                                    acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGc     >  2:1376901/1‑90 (MQ=255)
                                                                    acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGc     >  1:882047/1‑90 (MQ=255)
                                                                    acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCATCGATATCGACGCCGCGCGTTGGCTCGTCGAGc     >  2:1195516/1‑90 (MQ=255)
                                                                     cGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGca    >  2:1297124/1‑90 (MQ=255)
                                                                       cTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGcaac  <  2:1305914/90‑1 (MQ=255)
                                                                            |                                                                                    
CGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAAC  >  NZ_CP009273/2670283‑2670443

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ 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
TGAAACTCCTGCGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCTGGCTG  >  NZ_CP009273/2670261‑2670455
                                                                                                  |                                                                                                
TGAAACTCCTGCGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCC                                                                                                 >  SRR3722078.294887/1‑100 (MQ=60)
         TGCGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAG                                                                                        <  SRR3722078.234558/100‑1 (MQ=60)
                      CGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGAC                                                                           <  SRR3722078.302883/100‑1 (MQ=60)
                          CTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATA                                                                       >  SRR3722078.9835/1‑100 (MQ=60)
                            GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACctgtct                                                                     >  SRR3722078.140316/1‑94 (MQ=60)
                            GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG                                                                     >  SRR3722078.435250/1‑100 (MQ=60)
                            GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG                                                                     >  SRR3722078.1104137/1‑100 (MQ=60)
                            GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG                                                                     >  SRR3722078.915948/1‑100 (MQ=60)
                            GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG                                                                     >  SRR3722078.841731/1‑100 (MQ=60)
                            GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG                                                                     >  SRR3722078.567403/1‑100 (MQ=60)
                            GTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACG                                                                     >  SRR3722078.139739/1‑100 (MQ=60)
                                TAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAA                                                                 <  SRR3722078.194193/100‑1 (MQ=60)
                                   aagagacaGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATC                                                              <  SRR3722078.202418/92‑1 (MQ=60)
                                   CAGGATGCGGTCACACACCAGCGGCGACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATC                                                              <  SRR3722078.1447194/100‑1 (MQ=60)
                                     GGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTG                                                            >  SRR3722078.1326749/1‑100 (MQ=60)
                                        TGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTG                                                         <  SRR3722078.664744/100‑1 (MQ=60)
                                            GTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTG                                                     >  SRR3722078.430838/1‑100 (MQ=60)
                                            GTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTG                                                     >  SRR3722078.1373971/1‑100 (MQ=60)
                                            GTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTG                                                     >  SRR3722078.49212/1‑100 (MQ=60)
                                                CACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTT                                                 >  SRR3722078.527358/1‑100 (MQ=60)
                                                    CCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCCCGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGAT                                             <  SRR3722078.1349491/100‑1 (MQ=60)
                                                           CAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACctgtctcttatacacatctgacgctgccgacgatggc                                      >  SRR3722078.290796/1‑63 (MQ=60)
                                                                           CTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGT                      >  SRR3722078.1239094/1‑100 (MQ=60)
                                                                                GATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGC                 >  SRR3722078.893900/1‑100 (MQ=60)
                                                                                           CGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCTG      >  SRR3722078.655399/1‑100 (MQ=60)
                                                                                               TTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCTGGCTG  >  SRR3722078.766486/1‑100 (MQ=60)
                                                                                                  |                                                                                                
TGAAACTCCTGCGAGAACGTGCCGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAACAAAATCTGGCTG  >  NZ_CP009273/2670261‑2670455

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