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 3,333,763 +G coding (5/615 nt) serB2 → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,333,7601.G100.0% 74.1 / NA 20M1R (ATG→AGG) †serB2phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (10/10);  total (10/10)

CCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGC  >  AM260479/3333633‑3333888
                                                                                                                                |                                                                                                                                
ccAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATggggcc                                                                                                                            >  1:109997/1‑131 (MQ=255)
      cGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTcc                                                                                                                       <  1:748610/134‑1 (MQ=255)
      cGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTc                                                                                                                        <  2:46787/133‑1 (MQ=255)
                gcctgcctgcctGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTcc                                                                                                                       >  2:550945/1‑124 (MQ=255)
                gcctgcctgcctGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTcc                                                                                                                       <  1:43928/124‑1 (MQ=255)
                gcctgcctgcctGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTcc                                                                                                                       <  1:550945/124‑1 (MQ=255)
                gcctgcctgcctGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTcc                                                                                                                       >  2:43928/1‑124 (MQ=255)
                              ccgcccgACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTcgacg                                                                                               <  1:39244/134‑1 (MQ=255)
                              ccgcccgACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTcgac                                                                                                >  1:353841/1‑133 (MQ=255)
                                                            gccgccCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAg                                                                  <  2:353841/133‑1 (MQ=255)
                                                                 ccATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAgaagaa                                                                         >  1:966221/1‑121 (MQ=255)
                                                                 ccATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAgaagaa                                                                         <  2:966221/121‑1 (MQ=255)
                                                                                              cacaCCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGatat                               >  1:727737/1‑134 (MQ=255)
                                                                                                cacCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATc                              >  1:23566/1‑133 (MQ=255)
                                                                                                           cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTg                  >  2:354172/1‑134 (MQ=255)
                                                                                                           cgaccgacAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTg                  >  2:44582/1‑134 (MQ=255)
                                                                                                                 acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTa                                                             >  2:387668/1‑85 (MQ=255)
                                                                                                                 acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTa                                                             >  2:467560/1‑85 (MQ=255)
                                                                                                                 acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTa                                                             <  1:467560/85‑1 (MQ=255)
                                                                                                                 acAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTa                                                             <  1:387668/85‑1 (MQ=255)
                                                                                                                           tGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGgcgc  <  2:109997/134‑1 (MQ=255)
                                                                                                                                |                                                                                                                                
CCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGC  >  AM260479/3333633‑3333888

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
GCGCCCATACCGCACCCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGCCGCTGGCGGCCAACCCGCGCGACC  >  AM260479/3333618‑3333912
                                                                                                                                                  |                                                                                                                                                     
GCGCCCATACCGCACCCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCG                                                                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:1108:1914:43646/1‑150 (MQ=60)
        ACCGCACCCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGgcctgtctctt                                                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2205:16032:32733/1‑138 (MQ=60)
           gtttgttagctttcttttgtttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTC                                                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2216:18311:69410/54‑1 (MQ=60)
                    GCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2114:12468:31109/150‑1 (MQ=60)
                    GCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1102:17079:3675/150‑1 (MQ=60)
                    GCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2103:16509:97495/150‑1 (MQ=60)
                              CGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACG                                                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:1206:19775:48530/1‑149 (MQ=60)
                                            CCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGC                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1105:9391:72399/150‑1 (MQ=60)
                                                                 AGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACG                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2212:15671:96999/1‑150 (MQ=60)
                                                                              cgtctcgtgggctcggagatgtgtataagagacagCCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGC                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1215:19470:2046/115‑1 (MQ=60)
                                                                              cgtctcgtgggctcggagatgtgtataagagacagCCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGC                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1209:9962:56818/115‑1 (MQ=60)
                                                                                              CTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATC                                                      >  GWNJ‑0478:712:GW2002102894th:2:2115:20368:17829/1‑150 (MQ=60)
                                                                                                CTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCC                                                     >  GWNJ‑0478:712:GW2002102894th:2:1103:10072:68324/1‑149 (MQ=60)
                                                                                                                                      ACCATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGCCGCTGGCGGCC               >  GWNJ‑0478:712:GW2002102894th:2:1212:17187:32199/1‑149 (MQ=60)
                                                                                                                                         ATGAATGGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGCCGCTGGCGGCCAACC           <  GWNJ‑0478:712:GW2002102894th:2:2208:5350:12698/150‑1 (MQ=60)
                                                                                                                                                  GCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGCCGCTGGCGGCCAACCCGCGCGACC  >  GWNJ‑0478:712:GW2002102894th:2:2110:5968:95729/1‑150 (MQ=60)
                                                                                                                                                  GCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGCCGCTGGCGGCCAACCCGCGCGACC  <  GWNJ‑0478:712:GW2002102894th:2:1116:2540:35220/150‑1 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
GCGCCCATACCGCACCCAATGCGGCCCCGGCGCCTGCCTGCCTGCCCGCCCGACCTCTCGCCTTAAGGTAAAATCGCCGCCCATGAATCTGGCACTCTTTGACCTCGACCACACCCTGATCCCGACCGACAGCGACCATGAATGGGCCGCTTCCTGGTCCGCCTGGGCGTCGTCGACGAGGAAACCTACCGGCAGAAGAACGACGAGTTCTACGGCCACTACAAAGCCGGCACCCTCGATATCCAGGCCTTCCTGCGCTTTGCGCTGGCGCCGCTGGCGGCCAACCCGCGCGACC  >  AM260479/3333618‑3333912

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