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 F1 I1 R1 152 24.9 2083062 91.6% 1908084 109.2

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% 57.5 / NA 16intergenic (+162/‑58)cafA2/moaA1Ribonuclease G and E/molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (10/6);  total (10/6)

CACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCAG  >  AM260479/2787800‑2788042
                                                                                                                         |                                                                                                                          
cACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGc                                                                                                                >  2:546783/1‑134 (MQ=255)
              ttttCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAgg                                                                                                    <  2:174726/132‑1 (MQ=255)
                         aGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGAc                                                                                       <  1:546783/134‑1 (MQ=255)
                            cggccggATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGcc                                                                                    >  1:892340/1‑134 (MQ=255)
                                   aTGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGccca                                                                             <  2:892340/134‑1 (MQ=255)
                                                                              gCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCa                                                   >  1:409438/1‑117 (MQ=255)
                                                                              gCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCa                                                   <  2:409438/117‑1 (MQ=255)
                                                                                 tCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGCTTTTCGACCTgcg                                >  1:868783/1‑133 (MQ=255)
                                                                                 tCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTgcg                                <  2:868783/133‑1 (MQ=255)
                                                                                       aGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCAt                         >  2:753582/1‑134 (MQ=255)
                                                                                        gAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCAt                         <  1:753582/133‑1 (MQ=255)
                                                                                                  agAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTc              >  1:240337/1‑134 (MQ=255)
                                                                                                   gAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTc              >  2:920187/1‑133 (MQ=255)
                                                                                                      ggCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGAt           >  1:735831/1‑133 (MQ=255)
                                                                                                       gCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGa         >  1:449953/1‑134 (MQ=255)
                                                                                                               cccTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCAg  >  1:613083/1‑133 (MQ=255)
                                                                                                                         |                                                                                                                          
CACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCAG  >  AM260479/2787800‑2788042

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCAAAGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCAGCATTCCCGGCAAGCTGGCT  >  AM260479/2787779‑2788061
                                                                                                                                                 |                                                                                                                                          
CCAAAGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGAT                                                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:2113:5208:98311/150‑1 (MQ=60)
                                  TTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCG                                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:2207:11733:59604/1‑150 (MQ=60)
                                             CAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCAC                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2103:20964:68208/150‑1 (MQ=60)
                                                                                    GGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCctg                                                    >  GWNJ‑0478:712:GW2002102894th:2:1211:16710:80163/1‑147 (MQ=60)
                                                                                       TGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGCTTTTCGACCTGCGCG                                                 >  GWNJ‑0478:712:GW2002102894th:2:2206:11431:9503/1‑150 (MQ=60)
                                                                                                        GGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCG                                >  GWNJ‑0478:712:GW2002102894th:2:1114:9059:94167/1‑150 (MQ=60)
                                                                                                            AGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATG                             <  GWNJ‑0478:712:GW2002102894th:2:2114:9957:58686/149‑1 (MQ=60)
                                                                                                            AGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATG                             >  GWNJ‑0478:712:GW2002102894th:2:2115:9137:61278/1‑149 (MQ=60)
                                                                                                             GAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGAC                           >  GWNJ‑0478:712:GW2002102894th:2:1212:16260:3504/1‑150 (MQ=60)
                                                                                                                     GCAGAAGGCATCCGGCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCAGC                    >  GWNJ‑0478:712:GW2002102894th:2:2106:12146:26457/1‑149 (MQ=60)
                                                                                                                                   GCCCTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCAGCATTCCCtgtctctt      >  GWNJ‑0478:712:GW2002102894th:2:1104:18694:58580/1‑141 (MQ=60)
                                                                                                                                      CTCGAGCGCCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCctgtctcttatacacatctcc  >  GWNJ‑0478:712:GW2002102894th:2:2108:2295:12918/1‑129 (MQ=60)
                                                                                                                                                 |                                                                                                                                          
CCAAAGCAAGAAGGGGACGGCCACGCCGTCTCCTTTTTTCATGGGCAGCCGGCCGGATGAATCCCCGCCATCTCCGGGGGATGCGGGTGCCCCACCCCCGCCTCGGCTAGAATGGCAGCAGAAGGCATCCGGCCCTCGAGCGCCCTGATTTGGCCACGCCGCCAGGCCGCTCCCCGGACGCCGTGCCCACCCACAGGCCATGACCGAATCCGTCATTCCGATTTTCGACCTGCGCGACCATCGCTACCGCTCGATGACGCCCAGCATTCCCGGCAAGCTGGCT  >  AM260479/2787779‑2788061

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