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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,153,816 T→C S84P (TCA→CCA)  btuB → TonB‑dependent vitamin B12 receptor BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,153,8160TC100.0% 147.1 / NA 45S84P (TCA→CCA) btuBTonB‑dependent vitamin B12 receptor BtuB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/22);  total (23/22)

CTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153734‑4153901
                                                                                  |                                                                                     
cTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGGTCAGGGCAGCTCCCACCTAt                                                                                >  1:761470/1‑90 (MQ=255)
     aGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGCAGCTCCCCTCTAtttttt                                                                           >  1:119932/1‑89 (MQ=255)
      gTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTAt                                                                          >  2:66749/1‑90 (MQ=255)
       tCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTAtttttttt                                                                         >  1:1067799/1‑90 (MQ=255)
        cGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTc                                                                        <  1:1699/90‑1 (MQ=255)
          aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                                      <  2:900243/90‑1 (MQ=255)
          aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                                      <  2:767994/90‑1 (MQ=255)
          aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                                      <  1:794494/90‑1 (MQ=255)
          aCCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTcgc                                                                      <  1:264723/90‑1 (MQ=255)
            cTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCgg                                                                    <  1:64040/90‑1 (MQ=255)
                gTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGGACa                                                                >  2:484953/1‑90 (MQ=255)
                 tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTGTTTTTATTCGCGGTACaa                                                               <  1:1134922/90‑1 (MQ=255)
                 tCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACaa                                                               <  1:541836/90‑1 (MQ=255)
                   aaTGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAAt                                                             <  1:3778/90‑1 (MQ=255)
                    atgatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATg                                                            >  2:75030/1‑90 (MQ=255)
                      gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  2:1182363/90‑1 (MQ=255)
                      gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  1:562269/90‑1 (MQ=255)
                      gatgTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGcc                                                          <  1:619452/90‑1 (MQ=255)
                         gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGCAGCTCCCATCTATTTTTATTTGCGGGACAAAAGCCAGt                                                       >  1:503973/1‑90 (MQ=255)
                         gtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGTACAAAAGCCCGt                                                       >  2:1229799/1‑90 (MQ=255)
                          tgctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGGACAAATGCCAGTc                                                      >  2:915234/1‑90 (MQ=255)
                          tgctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTCGCGGTACAAATGCCAGTc                                                      >  2:511179/1‑90 (MQ=255)
                           gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                                     <  1:718854/90‑1 (MQ=255)
                           gctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCa                                                     <  1:1121524/90‑1 (MQ=255)
                              gcgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTTGCGGGACAAAAGCCAGTCAtgt                                                  >  1:369287/1‑90 (MQ=255)
                                gcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATgtgt                                                <  1:930528/90‑1 (MQ=255)
                                   gTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTgg                                             >  2:477753/1‑90 (MQ=255)
                                   gTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCAGTCATGTGTTgg                                             >  1:105647/1‑90 (MQ=255)
                                             cGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGGACCAATGCCAGTCATGTGTTGGGGTTTATTGa                                   >  2:1159820/1‑90 (MQ=255)
                                             cGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCCGGCATGTGTTGGTGTTTATTGa                                   >  1:678402/1‑90 (MQ=255)
                                             cGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTTGCGGGACAAAAGCCAGTCATGTGTTTGGGTTAATTGa                                   >  1:243758/1‑90 (MQ=255)
                                             cGTCGATATCACCCAAAACGGCGGTTCACGGCCGCTCCCCTCCATTTTTTTTTGCGGGACAAAAGCCAGTCCTGGGTTTGTGTTTATTGa                                   >  1:633700/1‑90 (MQ=255)
                                              gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                                  <  2:1010686/90‑1 (MQ=255)
                                              gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                                  <  1:158047/90‑1 (MQ=255)
                                              gTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCAGCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGAt                                  <  2:909181/90‑1 (MQ=255)
                                                 gATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGGGTTAATTGATGGc                               >  2:784399/1‑90 (MQ=255)
                                                  atatCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTATTTTTTTTTGCGGGACAAATGCCAGTCATGTGTTGGGGGTAATTTATGGCg                              >  1:376486/1‑90 (MQ=255)
                                                            aaaCGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGTACAAAAGCCCGTCATGTGTTGGGGTTTATTTTTGGCGTACGCCTGaa                    >  1:632182/1‑90 (MQ=255)
                                                             aaCGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAAt                   <  2:214989/90‑1 (MQ=255)
                                                                  cggTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGc              >  2:1186607/1‑90 (MQ=255)
                                                                    gTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCgg            <  1:957373/90‑1 (MQ=255)
                                                                          gTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGa      <  2:915102/90‑1 (MQ=255)
                                                                            cAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGt    <  1:5230/90‑1 (MQ=255)
                                                                              gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg  >  1:286744/1‑90 (MQ=255)
                                                                              gCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTgg  >  1:1106321/1‑90 (MQ=255)
                                                                                  |                                                                                     
CTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG  >  NZ_CP009273/4153734‑4153901

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCG  >  NZ_CP009273/4153717‑4153909
                                                                                                   |                                                                                             
CGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCC                                                                                               >  SRR3722089.1097517/1‑100 (MQ=60)
       ATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGGTCAGGGCAGCTCCCACCTAT                                                                                        >  SRR3722089.771656/1‑100 (MQ=60)
         ATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTT                                                                                      <  SRR3722089.1049164/100‑1 (MQ=60)
          TCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTT                                                                                     <  SRR3722089.1158439/100‑1 (MQ=60)
            GACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGCAGCTCCCCTCTATTTTTT                                                                                   >  SRR3722089.121398/1‑100 (MQ=60)
              CCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTT                                                                                 >  SRR3722089.1083832/1‑100 (MQ=60)
                         CGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAA                                                                      <  SRR3722089.1715/100‑1 (MQ=60)
                           ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATG                                                                    <  SRR3722089.805441/100‑1 (MQ=60)
                           ACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATG                                                                    <  SRR3722089.268110/100‑1 (MQ=60)
                             CTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCC                                                                  <  SRR3722089.64821/100‑1 (MQ=60)
                                GGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTTAGGGCAGCTCCCATCTATTTTTATTTGCGGGACAAAAGCCAGT                                                               >  SRR3722089.509918/1‑100 (MQ=60)
                                  TCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTGTTTTTATTCGCGGTACAAATGCCAGTCA                                                             <  SRR3722089.1152184/100‑1 (MQ=60)
                                  TCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCA                                                             <  SRR3722089.548144/100‑1 (MQ=60)
                                    AATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATG                                                           <  SRR3722089.3815/100‑1 (MQ=60)
                                     ATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTTGCGGGACAAAAGCCAGTCATGT                                                          >  SRR3722089.373826/1‑100 (MQ=60)
                                       GATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGT                                                        <  SRR3722089.568849/100‑1 (MQ=60)
                                       GATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGT                                                        <  SRR3722089.627230/100‑1 (MQ=60)
                                          GTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCAGTCATGTGTTGG                                                     >  SRR3722089.106948/1‑100 (MQ=60)
                                            GCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTG                                                   <  SRR3722089.1138529/100‑1 (MQ=60)
                                            GCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTG                                                   <  SRR3722089.728198/100‑1 (MQ=60)
                                                 GCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAAT                                              <  SRR3722089.944119/100‑1 (MQ=60)
                                                    GTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCCTCTATTTTTTTTTGCGGGACAAAAGCCCGGCATGTGTTGGTGTTTATTGA                                           >  SRR3722089.687205/1‑100 (MQ=52)
                                                    GTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGGCAGCTCCCATCTATTTTTTTTTGCGGGACAAAAGCCAGTCATGTGTTTGGGTTAATTGA                                           >  SRR3722089.246867/1‑100 (MQ=60)
                                                         CCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCCTCTATTTTTTTTTGCGGGACAAATGCCAGTCATGTGTTGGGGGTAATTTATGGCG                                      >  SRR3722089.381090/1‑100 (MQ=60)
                                                               GTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCC                                <  SRR3722089.160004/100‑1 (MQ=60)
                                                                   ATATCACCCAAAACGGCGGTTCAGGTCAGCTCCCATCTATTTTTTTTCGCGGTACAAAAGCCCGTCATGTGTTGGGGTTTATTTTTGGCGTACGCCTGAA                            >  SRR3722089.640214/1‑100 (MQ=60)
                                                                                     GTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG          >  SRR3722089.290413/1‑100 (MQ=60)
                                                                                     GTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG          >  SRR3722089.1123052/1‑100 (MQ=60)
                                                                                     GTTCAGGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGG          <  SRR3722089.971458/100‑1 (MQ=60)
                                                                                          GGTCAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTG     >  SRR3722089.1004658/1‑100 (MQ=60)
                                                                                             CAGCTCCCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCG  <  SRR3722089.5279/100‑1 (MQ=60)
                                                                                                   |                                                                                             
CGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCG  >  NZ_CP009273/4153717‑4153909

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