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 1,021,663 (A)5→4 coding (1737/2055 nt) helD → DNA helicase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,021,6590A.98.0% 201.9 / ‑1.5 50coding (1733/2055 nt)helDDNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (33/16);  total (34/16)
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.00e-01

TGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC  >  NZ_CP009273/1021575‑1021744
                                                                                    |                                                                                     
tGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaagc                                                                                  <  2:152013/90‑3 (MQ=255)
      aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAac                                                                           >  2:1566458/1‑90 (MQ=255)
      aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAac                                                                           >  2:895404/1‑90 (MQ=255)
      aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAac                                                                           >  2:586738/1‑90 (MQ=255)
           tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCt                                                                      >  2:1477782/1‑90 (MQ=255)
           tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCt                                                                      >  2:1362501/1‑90 (MQ=255)
           tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCt                                                                      >  2:1008477/1‑90 (MQ=255)
                 tATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGa                                                                >  2:478247/1‑90 (MQ=255)
                  aTGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGaa                                                               >  2:961018/1‑90 (MQ=255)
                   tGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAg                                                              <  1:586738/90‑1 (MQ=255)
                   tGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAg                                                              <  1:1566458/90‑1 (MQ=255)
                   tGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAg                                                              >  1:688844/1‑90 (MQ=255)
                        aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGa                                                                >  2:1193956/1‑83 (MQ=255)
                        aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGa                                                                <  1:1193956/83‑1 (MQ=255)
                        aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCa                                                         >  1:292237/1‑90 (MQ=255)
                        aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCa                                                         >  2:1130960/1‑90 (MQ=255)
                            ggAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATc                                                     >  1:816957/1‑90 (MQ=255)
                             gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    <  1:143564/90‑1 (MQ=255)
                             gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    <  1:895404/90‑1 (MQ=255)
                               agagCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGAc                                                  <  2:79009/90‑1 (MQ=255)
                                  gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                               >  1:1135165/1‑90 (MQ=255)
                                  gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                               >  2:65947/1‑90 (MQ=255)
                                  gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                               >  2:1225484/1‑90 (MQ=255)
                                           gATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCAtt                                      >  1:627117/1‑90 (MQ=255)
                                           gATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCAtt                                      >  1:464631/1‑90 (MQ=255)
                                            aTCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTc                                     >  2:1279646/1‑90 (MQ=255)
                                               cTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATg                                  >  2:621231/1‑90 (MQ=255)
                                                tGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGc                                 >  1:726731/1‑90 (MQ=255)
                                                 ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGc                                 >  1:636536/1‑90 (MQ=255)
                                                 ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    >  1:1270303/1‑70 (MQ=255)
                                                 ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                                    <  2:1270303/70‑1 (MQ=255)
                                                 ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                                <  1:287591/90‑1 (MQ=255)
                                                 ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                                >  1:1485347/1‑90 (MQ=255)
                                                 ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                                >  2:608340/1‑90 (MQ=255)
                                                      gTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa                           >  2:543938/1‑90 (MQ=255)
                                                      gTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa                           <  2:1095405/90‑1 (MQ=255)
                                                        tcccATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAg                         <  1:278728/88‑1 (MQ=255)
                                                             tCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCaa                    <  2:726731/90‑1 (MQ=255)
                                                                  tGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGc               <  1:1052310/90‑1 (MQ=255)
                                                                  tGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGc               <  1:759476/90‑1 (MQ=255)
                                                                    aGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCgg             >  1:516166/1‑90 (MQ=255)
                                                                     ggCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGa            >  2:855309/1‑90 (MQ=255)
                                                                       ccTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAtt          >  2:1321692/1‑90 (MQ=255)
                                                                       ccTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAtt          >  1:354090/1‑90 (MQ=255)
                                                                       ccTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGAtt          <  2:938314/90‑1 (MQ=255)
                                                                        cTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa         >  2:577526/1‑90 (MQ=255)
                                                                        cTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa         >  1:936338/1‑90 (MQ=255)
                                                                          gCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACg       >  1:1134515/1‑90 (MQ=255)
                                                                              gCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAAc                   <  2:1303825/74‑1 (MQ=255)
                                                                              gCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAAc                   >  1:1303825/1‑74 (MQ=255)
                                                                               ccTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGtcatc  <  2:816957/90‑1 (MQ=255)
                                                                                    |                                                                                     
TGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC  >  NZ_CP009273/1021575‑1021744

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
TCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCAT  >  NZ_CP009273/1021584‑1021743
                                                                               |                                                                                
TCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAG                                                             >  SRR3722078.696981/1‑100 (MQ=60)
     GGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCA                                                        >  SRR3722078.295802/1‑100 (MQ=60)
        TATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAAT                                                     <  SRR3722078.1212656/100‑1 (MQ=60)
         ATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATC                                                    >  SRR3722078.827762/1‑100 (MQ=60)
          TGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCG                                                   <  SRR3722078.1591963/100‑1 (MQ=60)
          TGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCG                                                   <  SRR3722078.593571/100‑1 (MQ=60)
               AACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTT                                              >  SRR3722078.1152564/1‑100 (MQ=60)
                    GAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGAC                                         <  SRR3722078.907436/100‑1 (MQ=60)
                    GAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGAC                                         <  SRR3722078.145228/100‑1 (MQ=60)
                        AGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATT                                     >  SRR3722078.634332/1‑100 (MQ=60)
                        AGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATT                                     >  SRR3722078.470174/1‑100 (MQ=60)
                             ATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGC                                >  SRR3722078.735798/1‑100 (MQ=60)
                              TTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGC                                >  SRR3722078.643850/1‑100 (MQ=60)
                              TTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCtgtctctta                               >  SRR3722078.1290455/1‑91 (MQ=60)
                              TTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCC                               >  SRR3722078.1509564/1‑100 (MQ=60)
                                        GGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGC                     <  SRR3722078.291098/100‑1 (MQ=60)
                                               TCCCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGC              <  SRR3722078.282133/100‑1 (MQ=60)
                                                 CCATCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGG            >  SRR3722078.522273/1‑100 (MQ=60)
                                                    TCACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATT         >  SRR3722078.358409/1‑100 (MQ=60)
                                                     CACATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTA        >  SRR3722078.949312/1‑100 (MQ=60)
                                                       CATGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACG      >  SRR3722078.1151894/1‑100 (MQ=60)
                                                         TGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTC    <  SRR3722078.769228/100‑1 (MQ=60)
                                                         TGAGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTC    <  SRR3722078.1068099/100‑1 (MQ=60)
                                                           AGGCCTGCCAGCCTGGAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACtgtct  >  SRR3722078.1324612/1‑95 (MQ=60)
                                                                               |                                                                                
TCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCAT  >  NZ_CP009273/1021584‑1021743

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