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
A0 F0 I1 R1 149 26.7 2337422 87.5% 2045244 105.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,787,923 +C intergenic (+165/‑55) cafA2 → / → moaA1 Ribonuclease G and E/molybdopterin biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,787,9201.C100.0% 42.3 / NA 13intergenic (+162/‑58)cafA2/moaA1Ribonuclease G and E/molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (7/6);  total (7/6)

GGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGC  >  AM260479/2787793‑2788038
                                                                                                                                |                                                                                                                      
ggACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGcccct                                                                                                                   >  1:105262/1‑131 (MQ=255)
                    tttttCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGc                                                                                               >  1:450516/1‑134 (MQ=255)
                                                                               aCCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGcccac                                                                        >  1:818585/1‑98 (MQ=255)
                                                                               aCCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGcccac                                                                        <  2:818585/98‑1 (MQ=255)
                                                                                   ccGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGAcc                                <  2:930613/134‑1 (MQ=255)
                                                                                          ggCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGa                         >  2:989782/1‑134 (MQ=255)
                                                                                          ggCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGa                         >  2:989783/1‑134 (MQ=255)
                                                                                                         agAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGcccacc                                                                       >  1:230651/1‑73 (MQ=255)
                                                                                                         agAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGcccacc                                                                       <  2:230651/73‑1 (MQ=255)
                                                                                                         agAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGa                                                           <  1:715173/85‑1 (MQ=255)
                                                                                                         agAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGa                                                           >  2:715173/1‑85 (MQ=255)
                                                                                                                 tCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAg                                                                  >  1:611985/1‑70 (MQ=255)
                                                                                                                 tCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAg                                                                  <  2:611985/70‑1 (MQ=255)
                                                                                                                 tCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATGCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGc  <  2:3268/134‑1 (MQ=255)
                                                                                                                                |                                                                                                                      
GGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGC  >  AM260479/2787793‑2788038

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCCAAAGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGC  >  AM260479/2787778‑2788038
                                                                                                                                                  |                                                                                                                   
GCCAAAGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCctg                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:1106:8437:56428/1‑146 (MQ=60)
  gtgtctcgtgggctcggagatgtgtataagagacagTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATT                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2215:6702:89108/114‑1 (MQ=60)
  CAAAGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATT                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2203:2743:87695/150‑1 (MQ=60)
     AGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGG                                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2208:18312:26279/1‑150 (MQ=60)
                    GCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCC                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1211:16955:95212/1‑150 (MQ=60)
                                                ttcctggagccttgcgggtttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagcCCCTGATTTGGCCACGCAGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCAC                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1216:13206:33807/52‑1 (MQ=60)
                                                  gttgttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGC                                                                <  GWNJ‑0478:712:GW2002102894th:2:2204:4133:38215/69‑1 (MQ=60)
                                                                       gtctcgtgggctcggagatgtgtataagagacagGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTC                                            <  GWNJ‑0478:712:GW2002102894th:2:2109:12488:40993/115‑1 (MQ=60)
                                                                               GGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCctgtctcttatacacatctccg                                   >  GWNJ‑0478:712:GW2002102894th:2:2115:10541:62460/1‑128 (MQ=60)
                                                                                                         GGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGACGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGctgtctcttatacacatctccgagcccacgagacactcgctaatctc         >  GWNJ‑0478:712:GW2002102894th:2:1115:13510:35729/1‑103 (MQ=60)
                                                                                                                 TGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgt  >  GWNJ‑0478:712:GW2002102894th:2:2102:7224:40186/1‑100 (MQ=60)
                                                                                                                                                  |                                                                                                                   
GCCAAAGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGC  >  AM260479/2787778‑2788038

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