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 155 25.8 2285122 88.0% 2010907 104.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 1,192,562 +C coding (163/594 nt) cccA → cytochrome c, mono‑or diheme variant

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604801,192,5601.C100.0% 95.0 / NA 26V54A (GTC→GCC) cccAcytochrome c, mono‑or diheme variant
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (13/13);  total (13/13)

GCCTCCGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCC  >  AM260480/1192430‑1192671
                                                                                                                                   |                                                                                                               
gcCTCCGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTccc                                                                                                              >  2:180626/1‑133 (MQ=255)
       ggCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGc                                                                                                        <  2:529153/134‑1 (MQ=255)
                 ggctggctGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCccggcc                                                                                              <  1:180626/134‑1 (MQ=255)
                 ggctggctGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCccggcc                                                                                              <  1:207101/134‑1 (MQ=255)
                            aTGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCt                                                                                   >  1:360590/1‑134 (MQ=255)
                                    tGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCg                                                                                       <  2:540716/122‑1 (MQ=255)
                                    tGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCg                                                                                       >  1:540716/1‑122 (MQ=255)
                                        gAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGa                                                                       >  1:752690/1‑134 (MQ=255)
                                         aaaCAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGAc                                                                      <  2:752690/134‑1 (MQ=255)
                                                      aacCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAgg                                                         <  2:9699/134‑1 (MQ=255)
                                                                      cTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGc                                         >  1:414636/1‑134 (MQ=255)
                                                                      cTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGc                                         >  1:174810/1‑134 (MQ=255)
                                                                         tAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCg                                        <  2:174810/132‑1 (MQ=255)
                                                                                cATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGAc                                                                      >  1:672755/1‑95 (MQ=255)
                                                                                cATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGAc                                                                      <  2:672755/95‑1 (MQ=255)
                                                                                     aTGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGccc                                                             >  2:1127775/1‑99 (MQ=255)
                                                                                     aTGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGccc                                                             <  1:1127775/99‑1 (MQ=255)
                                                                                       gAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATa                                                <  1:468764/110‑1 (MQ=255)
                                                                                       gAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATa                                                >  2:468764/1‑110 (MQ=255)
                                                                                       gAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTgg                           <  2:347867/131‑1 (MQ=255)
                                                                                       gAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTgg                           >  1:347867/1‑131 (MQ=255)
                                                                                       gAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGc                         >  1:532156/1‑133 (MQ=255)
                                                                                       gAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGc                         >  2:538677/1‑133 (MQ=255)
                                                                                       gAAAACGCAGCGCCGCCTGCTTGCCTTACTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGc                         >  1:634320/1‑133 (MQ=255)
                                                                                         aaaCGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGcgac                      <  2:360590/134‑1 (MQ=255)
                                                                                              cAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATg                 >  1:280558/1‑134 (MQ=255)
                                                                                                              ccttgcTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATcc  <  2:532156/133‑1 (MQ=255)
                                                                                                                                   |                                                                                                               
GCCTCCGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCC  >  AM260480/1192430‑1192671

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ 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
TGCTTTCAGGCGCCTCCGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGCACAAGGAGGCTCTGCGG  >  AM260480/1192419‑1192708
                                                                                                                                                |                                                                                                                                                  
cgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACC                                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2215:19592:69424/102‑1 (MQ=60)
                        GTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTG                                                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:1204:13591:87434/1‑150 (MQ=60)
                           CGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCC                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1110:14477:69759/150‑1 (MQ=60)
                           CGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCC                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1112:10139:8568/150‑1 (MQ=60)
                                GGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGG                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1214:17317:30023/1‑150 (MQ=60)
                                    GGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGAC                                                                                                           >  GWNJ‑0478:712:GW2002102894th:2:2110:10226:73352/1‑150 (MQ=60)
                                                            ggagatgtgtataagagacagCTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACC                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2216:4157:27901/129‑1 (MQ=60)
                                                                  ACTAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCG                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1208:3381:93551/1‑150 (MQ=60)
                                                                  ACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCG                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1110:16902:38905/1‑150 (MQ=60)
                                                                          aagagacagGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCAC                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1210:11439:79068/140‑1 (MQ=60)
                                                                            ggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAAC                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2214:11070:52424/97‑1 (MQ=60)
                                                                            GTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCctgtctcttatacacatctccgagc                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2106:13081:22712/1‑125 (MQ=60)
                                                                                   GTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCG                                                             >  GWNJ‑0478:712:GW2002102894th:2:1204:18958:10056/1‑149 (MQ=60)
                                                                                   GTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCG                                                             >  GWNJ‑0478:712:GW2002102894th:2:1215:5639:89622/1‑149 (MQ=60)
                                                                                   GTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTACTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCG                                                             >  GWNJ‑0478:712:GW2002102894th:2:2105:10328:89547/1‑149 (MQ=60)
                                                                                          ACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGG                                                     >  GWNJ‑0478:712:GW2002102894th:2:1116:9443:87908/1‑150 (MQ=60)
                                                                                                                            TGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACctgtctcttatacacatctccgagcccacgagacactcgctaatctcg                   >  GWNJ‑0478:712:GW2002102894th:2:1212:9184:97691/1‑102 (MQ=60)
                                                                                                                                  CGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCtgtctcttatacacatctcccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaac              >  GWNJ‑0478:712:GW2002102894th:2:1111:10922:100393/1‑71 (MQ=60)
                                                                                                                                            gtcCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGCACAAGGAGGCTCTGCG   <  GWNJ‑0478:712:GW2002102894th:2:1214:20855:20729/147‑1 (MQ=60)
                                                                                                                                             tcCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaacaagcctcgtagaaa   >  GWNJ‑0478:712:GW2002102894th:2:1214:13732:96604/3‑55 (MQ=60)
                                                                                                                                             tcCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgctt  >  GWNJ‑0478:712:GW2002102894th:2:2206:8281:48672/3‑85 (MQ=60)
                                                                                                                                             tcCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTgtctcttatacacatctccgagcccacgagacact  >  GWNJ‑0478:712:GW2002102894th:2:2101:15744:44995/3‑115 (MQ=60)
                                                                                                                                                |                                                                                                                                                  
TGCTTTCAGGCGCCTCCGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGCACAAGGAGGCTCTGCGG  >  AM260480/1192419‑1192708

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