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 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,882,413 A→C V349G (GTT→GGT)  cysI ← assimilatory sulfite reductase (NADPH) hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,882,4130AC93.3% 90.7 / ‑1.5 30V349G (GTT→GGT) cysIassimilatory sulfite reductase (NADPH) hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (12/16);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT  >  NZ_CP009273/2882324‑2882493
                                                                                         |                                                                                
ggCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTaa                                                                                  >  1:649079/1‑90 (MQ=255)
  ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                <  2:183725/90‑1 (MQ=255)
  ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                <  2:23503/90‑1 (MQ=255)
  ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                <  1:711425/90‑1 (MQ=255)
  ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                <  1:663795/90‑1 (MQ=255)
  ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                >  1:632392/1‑90 (MQ=255)
  ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                >  1:363299/1‑90 (MQ=255)
  ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                >  1:289513/1‑90 (MQ=255)
   cGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAcccc                                                                               <  2:399468/90‑1 (MQ=255)
     gTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAg                                                                             >  2:726410/1‑90 (MQ=255)
          cAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAAt                                                                        >  1:135112/1‑90 (MQ=255)
            gCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                <  2:678529/80‑1 (MQ=255)
            gCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                >  1:678529/1‑80 (MQ=255)
               gACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGAt                                                                   >  2:65294/1‑90 (MQ=255)
                   cgcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc                                                               <  1:83035/90‑1 (MQ=255)
                         aTAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTcc                                                         <  1:164987/90‑1 (MQ=255)
                           aaTCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCgg                                                       <  1:451018/90‑1 (MQ=255)
                              cAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCACGTCCGGTGa                                                    <  1:635241/90‑1 (MQ=255)
                                aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  <  1:701848/90‑1 (MQ=255)
                                aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  <  2:180756/90‑1 (MQ=255)
                                aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  <  2:553371/90‑1 (MQ=255)
                                aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCACCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  <  1:270222/90‑1 (MQ=255)
                                 ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACt                                                 >  2:592213/1‑90 (MQ=255)
                                   aTGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  <  2:836331/87‑1 (MQ=255)
                                   aTGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  >  1:836331/1‑87 (MQ=255)
                                     gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                             >  2:401515/1‑90 (MQ=255)
                                      cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                            >  1:758937/1‑90 (MQ=255)
                                            ttttCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACgg                                      >  2:300862/1‑90 (MQ=255)
                                               tCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                   <  2:441051/90‑1 (MQ=255)
                                                  aTAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGt                                >  1:646299/1‑90 (MQ=255)
                                                    aaaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTc                              <  2:135112/90‑1 (MQ=255)
                                                    aaaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTc                              <  2:646299/90‑1 (MQ=255)
                                                      aCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCaa                            <  1:421490/90‑1 (MQ=255)
                                                         gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAc                                        >  2:484060/1‑75 (MQ=255)
                                                         gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAc                                        <  1:484060/75‑1 (MQ=255)
                                                          cGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc                                                               <  1:741204/51‑1 (MQ=255)
                                                          cGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc                                                               >  2:741204/1‑51 (MQ=255)
                                                              aGGTGCCAGTTATCATCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATc                    <  2:506562/90‑1 (MQ=255)
                                                                        tATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc          <  2:632392/90‑1 (MQ=255)
                                                                                aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt  >  1:182511/1‑90 (MQ=255)
                                                                                      ttACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc          <  2:788438/76‑1 (MQ=255)
                                                                                         |                                                                                
GGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT  >  NZ_CP009273/2882324‑2882493

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
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAA  >  NZ_CP009273/2882314‑2882512
                                                                                                   |                                                                                                   
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAA                                                                                                     >  SRR3722079.655785/1‑100 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                   >  SRR3722079.291906/1‑100 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                   >  SRR3722079.366142/1‑100 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                   >  SRR3722079.638862/1‑100 (MQ=60)
   TCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCC                                                                                                  <  SRR3722079.373865/100‑1 (MQ=60)
   TCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCC                                                                                                  <  SRR3722079.347237/100‑1 (MQ=60)
      AGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGC                                                                                               <  SRR3722079.376827/100‑1 (MQ=60)
       GCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCC                                                                                              <  SRR3722079.62652/100‑1 (MQ=60)
          GGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAAT                                                                                           >  SRR3722079.136253/1‑100 (MQ=60)
            CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATAC                                                                                         <  SRR3722079.718984/100‑1 (MQ=60)
            CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATAC                                                                                         >  SRR3722079.685651/1‑100 (MQ=60)
            CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATAC                                                                                         <  SRR3722079.670700/100‑1 (MQ=60)
                             gtgggctcggagatgtgtataagagacagCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTG                                                                        <  SRR3722079.749131/71‑1 (MQ=60)
                             CGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTG                                                                        <  SRR3722079.83704/100‑1 (MQ=60)
                                   ATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCA                                                                  >  SRR3722079.845283/1‑100 (MQ=60)
                                   ATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCA                                                                  <  SRR3722079.166429/100‑1 (MQ=60)
                                     AATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATA                                                                <  SRR3722079.454965/100‑1 (MQ=60)
                                      ATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATAT                                                               >  SRR3722079.767055/1‑100 (MQ=60)
                                        CAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCACGTCCGGTGAACTCATATGG                                                             <  SRR3722079.641763/100‑1 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                           <  SRR3722079.647548/100‑1 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                           <  SRR3722079.709285/100‑1 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCACCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                           <  SRR3722079.272482/100‑1 (MQ=60)
                                                  ACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGT                                                   >  SRR3722079.652982/1‑100 (MQ=60)
                                                    gacagTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTC                                                 <  SRR3722079.488284/95‑1 (MQ=60)
                                                                ACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCC                                     <  SRR3722079.425122/100‑1 (MQ=60)
                                                                                GTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT                     >  SRR3722079.184107/1‑100 (MQ=60)
                                                                                      ATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCtgt               >  SRR3722079.796886/1‑97 (MQ=60)
                                                                                       TCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT              >  SRR3722079.346405/1‑100 (MQ=60)
                                                                                       TCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT              >  SRR3722079.213261/1‑100 (MQ=60)
                                                                                            GCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAAC         >  SRR3722079.510189/1‑100 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG        >  SRR3722079.305678/1‑100 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG        >  SRR3722079.522648/1‑100 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG        >  SRR3722079.474922/1‑100 (MQ=60)
                                                                                                  ACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCA   >  SRR3722079.312473/1‑100 (MQ=60)
                                                                                                   CCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAA  >  SRR3722079.227986/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAA  >  NZ_CP009273/2882314‑2882512

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