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 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% 93.7 / NA 25V54A (GTC→GCC) cccAcytochrome c, mono‑or diheme variant
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (12/13);  total (12/13)

CGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATG  >  AM260480/1192435‑1192688
                                                                                                                              |                                                                                                                                
cGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGAcc                                                                                                                           >  2:543250/1‑134 (MQ=255)
            ggctggctGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCccggcc                                                                                                               <  1:543250/134‑1 (MQ=255)
            ggctggctGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCccggcc                                                                                                               <  1:910027/134‑1 (MQ=255)
             gctggctgGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCccggccg                                                                                                              <  2:981369/134‑1 (MQ=255)
                ggctggTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCt                                                                                                           >  1:433633/1‑134 (MQ=255)
                                    aaaCAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGAc                                                                                       >  2:985941/1‑134 (MQ=255)
                                    aaaCAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGAc                                                                                       >  1:355639/1‑134 (MQ=255)
                                     aaCAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACg                                                                                      <  2:433633/134‑1 (MQ=255)
                                             aaaCAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGccc                                                                              <  2:227416/134‑1 (MQ=255)
                                              aacaacCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTcc                                                                                                                    >  1:1017394/1‑95 (MQ=255)
                                              aacaacCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTcc                                                                                                                    <  2:1017394/95‑1 (MQ=255)
                                                               cACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCa                                                              >  1:915546/1‑132 (MQ=255)
                                                                                aTGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGc                                                                                                                         <  1:1158411/56‑1 (MQ=255)
                                                                                aTGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGc                                                                                                                         >  2:1158411/1‑56 (MQ=255)
                                                                                aTGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTggg                                           >  2:854098/1‑134 (MQ=255)
                                                                                      aCGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGc                                                                         >  1:1068338/1‑98 (MQ=255)
                                                                                      aCGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGc                                                                         <  2:1068338/98‑1 (MQ=255)
                                                                                         cAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATg                                  <  1:854098/134‑1 (MQ=255)
                                                                                             gccgccTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGgcgc                              >  1:229509/1‑134 (MQ=255)
                                                                                               cgccTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGccc                            >  1:223072/1‑134 (MQ=255)
                                                                                                gccTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCg                           <  1:538071/134‑1 (MQ=255)
                                                                                                                cgcgCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACaa                      <  2:1133511/123‑1 (MQ=255)
                                                                                                                cgcgCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACaa                      >  1:1133511/1‑123 (MQ=255)
                                                                                                                   gcTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCg        <  1:2399/134‑1 (MQ=255)
                                                                                                                         cgcgTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATg  <  1:422733/134‑1 (MQ=255)
                                                                                                                              |                                                                                                                                
CGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATG  >  AM260480/1192435‑1192688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ 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
CAGGCGCCTCCGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGCACAAGGAGGCTCTGCGG  >  AM260480/1192425‑1192708
                                                                                                                                          |                                                                                                                                                  
cgtgggctcggagatgtgtataagagacagGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCC                                                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2112:3560:61100/117‑1 (MQ=60)
           agcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTC                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:2216:18727:55547/86‑1 (MQ=60)
           GGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTC                                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1211:10577:8047/1‑150 (MQ=60)
                     CGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCC                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2202:9958:54407/150‑1 (MQ=60)
                     CGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCC                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1214:7052:52115/150‑1 (MQ=60)
                       ttttgtttttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCC                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2210:17593:76022/62‑1 (MQ=60)
                               GTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACG                                                                                                          >  GWNJ‑0478:712:GW2002102894th:2:1204:5979:54369/1‑150 (MQ=60)
                                         TGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCctgtctcttatacacatctcgagcc                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2208:6142:74913/1‑125 (MQ=60)
                                                          CAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAA                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2202:12573:88593/1‑148 (MQ=60)
                                                                           CTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCctgaccctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaa                                                              >  GWNJ‑0478:712:GW2002102894th:2:2103:8323:24212/1‑63 (MQ=60)
                                                                                 GAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCctgtctcttatacacatctccg                                                        >  GWNJ‑0478:712:GW2002102894th:2:2212:21229:41562/1‑128 (MQ=60)
                                                                                        CCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCC                                                 >  GWNJ‑0478:712:GW2002102894th:2:1115:2129:30347/1‑150 (MQ=60)
                                                                                          ATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCG                                               >  GWNJ‑0478:712:GW2002102894th:2:1112:17775:97864/1‑150 (MQ=60)
                                                                                             attagcgagtgtctcgtgggctcggagatgtgtataagagacagCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGAC                                            <  GWNJ‑0478:712:GW2002102894th:2:2202:11256:12054/106‑1 (MQ=60)
                                                                                                  GCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCC                                       <  GWNJ‑0478:712:GW2002102894th:2:2116:18964:28252/150‑1 (MQ=60)
                                                                                                         CGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGT                                <  GWNJ‑0478:712:GW2002102894th:2:1214:5079:27643/150‑1 (MQ=60)
                                                                                                           CCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGT                                >  GWNJ‑0478:712:GW2002102894th:2:2215:17567:47926/1‑148 (MQ=60)
                                                                                                            CTGCTTGCCTTGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaagactccaaacat                             >  GWNJ‑0478:712:GW2002102894th:2:1115:3139:28199/1‑58 (MQ=60)
                                                                                                                      TGCTCGCGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGC                   >  GWNJ‑0478:712:GW2002102894th:2:1106:16427:73176/1‑150 (MQ=60)
                                                                                                                            CGCTTTGCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGCACAAGG             <  GWNJ‑0478:712:GW2002102894th:2:1101:8769:17308/150‑1 (MQ=60)
                                                                                                                                  GCGCGTCCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGCACAAGGAGGCTC       <  GWNJ‑0478:712:GW2002102894th:2:1209:14162:46394/150‑1 (MQ=60)
                                                                                                                                       tcCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTgtctcttatacacatctccgagcccacgagacact  >  GWNJ‑0478:712:GW2002102894th:2:2101:18627:41624/3‑115 (MQ=60)
                                                                                                                                          |                                                                                                                                                  
CAGGCGCCTCCGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGCGTCCTGACCGCCCTCCCGGCCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATGCGCACAAGGAGGCTCTGCGG  >  AM260480/1192425‑1192708

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