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 539,975 A→G D258D (GAT→GAC allC ← allantoate deiminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273539,9750AG100.0% 85.7 / NA 25D258D (GAT→GACallCallantoate deiminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/10);  total (15/10)

AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGT  >  NZ_CP009273/539890‑540061
                                                                                     |                                                                                      
aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                    >  2:501398/1‑90 (MQ=255)
          gAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCtttt                                                                          >  1:297445/1‑90 (MQ=255)
            aCGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCg                                                                        <  2:221039/90‑1 (MQ=255)
               tGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCtt                                                                     >  1:793062/1‑90 (MQ=255)
                         cGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgaccga                                                           >  2:552951/1‑90 (MQ=255)
                         cGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgaccga                                                           >  1:882711/1‑90 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAAt                                               >  2:373048/1‑90 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAAt                                               >  2:211735/1‑90 (MQ=255)
                                        cGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                            >  2:386361/1‑90 (MQ=255)
                                           aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                         >  2:467753/1‑90 (MQ=255)
                                           aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                         >  1:231435/1‑90 (MQ=255)
                                           aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                         >  2:596675/1‑90 (MQ=255)
                                             tCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                       <  1:368621/90‑1 (MQ=255)
                                               gggCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaaa                                     <  1:501398/90‑1 (MQ=255)
                                                  cgcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCg                                  <  2:412895/90‑1 (MQ=255)
                                                   gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaa                                      <  1:321338/85‑1 (MQ=255)
                                                   gcgGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaa                                      >  2:321338/1‑85 (MQ=255)
                                                       ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  2:753376/90‑1 (MQ=255)
                                                       ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  1:399179/90‑1 (MQ=255)
                                                              ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCa                      >  2:628382/1‑90 (MQ=255)
                                                               tGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCac                     <  1:247816/90‑1 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                    >  2:27367/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             <  2:177150/65‑1 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                             >  1:177150/1‑65 (MQ=255)
                                                                                  cGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCgt  <  1:163613/90‑1 (MQ=255)
                                                                                     |                                                                                      
AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGT  >  NZ_CP009273/539890‑540061

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ 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
GACAATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGCCTGCATGG  >  NZ_CP009273/539887‑540074
                                                                                        |                                                                                                   
gagatgtgtataagagacaGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCT                                                                                          <  SRR3722072.523676/81‑1 (MQ=60)
   AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTT                                                                                       >  SRR3722072.302600/1‑100 (MQ=60)
        ATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTT                                                                                  >  SRR3722072.806936/1‑100 (MQ=60)
                  TGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGA                                                                        >  SRR3722072.897910/1‑100 (MQ=60)
                                    CATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACT                                                      >  SRR3722072.235310/1‑100 (MQ=60)
                                                TCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAC                                          <  SRR3722072.375173/100‑1 (MQ=60)
                                                 CGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACT                                         <  SRR3722072.326954/100‑1 (MQ=60)
                                                  GGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTG                                        <  SRR3722072.510447/100‑1 (MQ=60)
                                                          CTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGA                                <  SRR3722072.406382/100‑1 (MQ=60)
                                                                  TGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACC                        <  SRR3722072.252000/100‑1 (MQ=60)
                                                                         GGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACtgtctcttatacac                 >  SRR3722072.180217/1‑86 (MQ=60)
                                                                              GAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGT            >  SRR3722072.592964/1‑100 (MQ=60)
                                                                                 CCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGCC         <  SRR3722072.917671/100‑1 (MQ=60)
                                                                                     CGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGCCTGCA     <  SRR3722072.166458/100‑1 (MQ=60)
                                                                                      GGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGCCTGCAT    <  SRR3722072.355734/100‑1 (MQ=60)
                                                                                        GTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGCCTGCATGG  <  SRR3722072.215340/100‑1 (MQ=60)
                                                                                        |                                                                                                   
GACAATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGCCTGCATGG  >  NZ_CP009273/539887‑540074

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