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 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,841,200 +G intergenic (‑2736/‑2680) h16_A1684 ← / → ycaC conserved hypothetical protein/YcaC related amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,841,1971.G100.0% 69.5 / NA 19intergenic (‑2733/‑2683)h16_A1684/ycaCconserved hypothetical protein/YcaC related amidohydrolase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/11);  total (8/11)

CAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCG  >  AM260479/1841098‑1841319
                                                                                                    |                                                                                                                          
cAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTc                                                                                           >  2:778466/1‑134 (MQ=255)
                tgatCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCgctggc                                                                           >  1:376562/1‑134 (MQ=255)
                     cggccggTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                       <  1:509535/133‑1 (MQ=255)
                            tACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGAcc                                                               <  2:233774/134‑1 (MQ=255)
                                          tCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggtg                                                                      <  1:276001/113‑1 (MQ=255)
                                          tCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggtg                                                                      >  2:276001/1‑113 (MQ=255)
                                            gTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCACAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCt                                               <  2:957322/134‑1 (MQ=255)
                                                ccGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCg                                           <  2:75932/134‑1 (MQ=255)
                                                       aGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCg                                    >  1:106952/1‑134 (MQ=255)
                                                            gccgccTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCAc                               <  2:71881/134‑1 (MQ=255)
                                                               gccTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAg                            <  2:106952/134‑1 (MQ=255)
                                                                    cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCgc                                                                               <  2:561671/78‑1 (MQ=255)
                                                                    cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCgc                                                                               >  1:561671/1‑78 (MQ=255)
                                                                    cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGc                                          >  2:666438/1‑115 (MQ=255)
                                                                    cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGc                                          <  1:666438/115‑1 (MQ=255)
                                                                    cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGc                       <  1:266895/134‑1 (MQ=255)
                                                                                         cAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCAc                               <  2:585485/105‑1 (MQ=255)
                                                                                         cAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCAc                               >  1:585485/1‑105 (MQ=255)
                                                                                         cAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCcgcg  >  1:996631/1‑134 (MQ=255)
                                                                                                    |                                                                                                                          
CAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCG  >  AM260479/1841098‑1841319

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCTGGGCCTTGGCGATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCGAGTAGT  >  AM260479/1841054‑1841345
                                                                                                                                                   |                                                                                                                                                 
TCTGGGCCTTGGCGATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGG                                                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2105:17522:40405/148‑1 (MQ=60)
                                   ttgttttttttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagTGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAG                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2108:11036:57317/62‑1 (MQ=60)
                                             AGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCT                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:1209:6379:12856/1‑150 (MQ=60)
                                                                gagacaGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCC                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1203:8294:67149/144‑1 (MQ=60)
                                                                CCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATC                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1216:14704:71037/149‑1 (MQ=60)
                                                                                    CATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGG                                                             >  GWNJ‑0478:712:GW2002102894th:2:1108:3497:15404/1‑150 (MQ=60)
                                                                                            gacagGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGG                                                     <  GWNJ‑0478:712:GW2002102894th:2:2110:12698:4274/145‑1 (MQ=60)
                                                                                                 GTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCtgtctcttatacacatctccgagcccacgagacggagctac                                                >  GWNJ‑0478:712:GW2002102894th:2:2105:17542:36586/1‑109 (MQ=60)
                                                                                                               GCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCT                                  <  GWNJ‑0478:712:GW2002102894th:2:1203:14740:2379/150‑1 (MQ=60)
                                                                                                                      GCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGctgtctcttatacac                           >  GWNJ‑0478:712:GW2002102894th:2:2104:13883:94665/1‑135 (MQ=60)
                                                                                                                      GCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGT                           >  GWNJ‑0478:712:GW2002102894th:2:2215:11806:27738/1‑150 (MQ=60)
                                                                                                                                     CAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCG            >  GWNJ‑0478:712:GW2002102894th:2:1214:1455:57574/1‑150 (MQ=60)
                                                                                                                                           GTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGctgtctcttatacacatctccgagcccacgagacgg      >  GWNJ‑0478:712:GW2002102894th:2:2202:13087:14238/1‑114 (MQ=60)
                                                                                                                                               GGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaa  >  GWNJ‑0478:712:GW2002102894th:2:1210:9861:28015/1‑71 (MQ=60)
                                                                                                                                                   |                                                                                                                                                 
TCTGGGCCTTGGCGATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCGAGTAGT  >  AM260479/1841054‑1841345

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