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 2,739,096 T→C E61G (GAA→GGA)  rimM ← ribosome maturation factor RimM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,739,0960TC97.5% 130.3 / ‑3.5 40E61G (GAA→GGA) rimMribosome maturation factor RimM
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (20/19);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

ACAATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGGAA  >  NZ_CP009273/2739011‑2739179
                                                                                     |                                                                                   
acaATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCttccagcttcccag                                                                                 >  1:1187928/1‑90 (MQ=255)
 caATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCttccagcttcccagc                                                                                <  2:595916/90‑1 (MQ=255)
   aTTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTg                                                                              >  1:1151423/1‑90 (MQ=255)
               gTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCa                                                                  <  1:639113/90‑1 (MQ=255)
                 tCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACt                                                                >  1:1131439/1‑90 (MQ=255)
                     cgcATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGAcc                                                            <  1:594085/90‑1 (MQ=255)
                     cgcATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGAcc                                                            <  1:461370/90‑1 (MQ=255)
                       cATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCg                                                          <  1:926373/90‑1 (MQ=255)
                       cATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCg                                                          <  1:569043/90‑1 (MQ=255)
                       cATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCg                                                          <  1:345114/90‑1 (MQ=255)
                        aTCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGc                                                         >  1:278796/1‑90 (MQ=255)
                        aTCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGc                                                         >  2:391018/1‑90 (MQ=255)
                            cGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTc                                                     <  2:539937/90‑1 (MQ=255)
                                 gTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                                <  1:411158/90‑1 (MQ=255)
                                 gTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                                <  1:53660/90‑1 (MQ=255)
                                 gTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                                <  1:829191/90‑1 (MQ=255)
                                 gTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAt                                                >  2:1076627/1‑90 (MQ=255)
                                   cAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATaa                                              >  1:154117/1‑90 (MQ=255)
                                     aCGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAAc                                            <  2:1013072/90‑1 (MQ=255)
                                       gCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCa                                          <  1:98712/90‑1 (MQ=255)
                                        ccTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAg                                         <  1:959729/90‑1 (MQ=255)
                                               gCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGAt                                  >  2:431749/1‑90 (MQ=255)
                                                  tgatgaTCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTc                                                     <  2:565543/68‑1 (MQ=255)
                                                  tgatgaTCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTc                                                     >  1:565543/1‑68 (MQ=255)
                                                  tgatgaTCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGt                               >  2:280731/1‑90 (MQ=255)
                                                  tgatgaTCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGt                               >  2:433283/1‑90 (MQ=255)
                                                  tgatgaTCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGt                               >  2:856129/1‑90 (MQ=255)
                                                  tgatgaTCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGt                               >  1:683870/1‑90 (MQ=255)
                                                                 gATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCgg                >  2:138669/1‑90 (MQ=255)
                                                                 gATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCgg                >  1:609929/1‑90 (MQ=255)
                                                                 gATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCgg                >  1:1033567/1‑90 (MQ=255)
                                                                 gATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCgg                <  1:306557/90‑1 (MQ=255)
                                                                 gATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCgg                >  1:538577/1‑90 (MQ=255)
                                                                   ttGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCg              <  1:418669/90‑1 (MQ=255)
                                                                     gtggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGAccc                                                           >  2:847055/1‑43 (MQ=255)
                                                                     gtggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGAccc                                                           <  1:847055/43‑1 (MQ=255)
                                                                      tggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCa                                          <  1:813766/59‑1 (MQ=255)
                                                                      tggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCa                                          >  2:813766/1‑59 (MQ=255)
                                                                      tggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCt           <  2:300237/90‑1 (MQ=255)
                                                                      tggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCt           >  1:862133/1‑90 (MQ=255)
                                                                      tggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCt           >  1:575019/1‑90 (MQ=255)
                                                                       ggtgCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCtt          >  2:478286/1‑90 (MQ=255)
                                                                               cagcttcccagctGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGgaa  >  1:644434/1‑90 (MQ=255)
                                                                               cagcttcccagctGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGgaa  <  1:695520/90‑1 (MQ=255)
                                                                                     |                                                                                   
ACAATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGGAA  >  NZ_CP009273/2739011‑2739179

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGACAATTTCACAATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGGAAGAAAACACTC  >  NZ_CP009273/2739001‑2739189
                                                                                               |                                                                                             
CGACAATTTCACAATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAG                                                                                           >  SRR3722089.1205999/1‑100 (MQ=60)
   CAATTTCACAATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTG                                                                                        >  SRR3722089.1168952/1‑100 (MQ=60)
                 GTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACT                                                                          >  SRR3722089.1148645/1‑100 (MQ=60)
                        GGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGC                                                                   >  SRR3722089.282355/1‑100 (MQ=60)
                         GTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCC                                                                  <  SRR3722089.647273/100‑1 (MQ=60)
                          actaagttggtggccggcagcgtcagatgtgtaaaagagacagTGTCCTGTTTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCT                                                                 <  SRR3722089.537197/57‑1 (MQ=60)
                               CGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGG                                                            <  SRR3722089.601323/100‑1 (MQ=60)
                               CGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGG                                                            <  SRR3722089.466878/100‑1 (MQ=60)
                                atggtcgtcggcagcgtcagatgtgtataagagacagTGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGA                                                           <  SRR3722089.859104/63‑1 (MQ=60)
                                 CATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAT                                                          <  SRR3722089.349409/100‑1 (MQ=60)
                                 CATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAT                                                          <  SRR3722089.939890/100‑1 (MQ=60)
                                 CATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGAT                                                          <  SRR3722089.575765/100‑1 (MQ=60)
                                   TCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAA                                                        >  SRR3722089.156029/1‑100 (MQ=60)
                                           GTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGC                                                <  SRR3722089.416127/100‑1 (MQ=60)
                                           GTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGC                                                <  SRR3722089.840902/100‑1 (MQ=60)
                                           GTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGC                                                <  SRR3722089.54323/100‑1 (MQ=60)
                                                 tcagatgtgtataagagacagGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAG                                          <  SRR3722089.825150/79‑1 (MQ=60)
                                                 GCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAG                                          <  SRR3722089.99942/100‑1 (MQ=60)
                                                  CCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCctgtctcttata                                         >  SRR3722089.572194/1‑88 (MQ=60)
                                                  CCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGT                                         >  SRR3722089.692756/1‑100 (MQ=60)
                                                  CCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGT                                         <  SRR3722089.973861/100‑1 (MQ=60)
                                                                 ATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGG                          >  SRR3722089.617520/1‑100 (MQ=60)
                                                                 ATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGG                          >  SRR3722089.544848/1‑100 (MQ=60)
                                                                 ATCATGTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGG                          >  SRR3722089.1049020/1‑100 (MQ=60)
                                                                      GTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCT                     >  SRR3722089.581835/1‑100 (MQ=60)
                                                                      GTCCTGATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCT                     >  SRR3722089.874468/1‑100 (MQ=60)
                                                                           GATTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGT                <  SRR3722089.310448/100‑1 (MQ=60)
                                                                             TTGTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGG              <  SRR3722089.423727/100‑1 (MQ=60)
                                                                               GTGGTGCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGGAA            >  SRR3722089.652675/1‑100 (MQ=60)
                                                                                    GCTTCCAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTGTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGGAAGAAAA       >  SRR3722089.499515/1‑100 (MQ=60)
                                                                                         CAGCTTCCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGGAAGAAAACACTC  <  SRR3722089.704587/100‑1 (MQ=60)
                                                                                               |                                                                                             
CGACAATTTCACAATTCGTCAGCAGGTTCGCCGCATCACGATCGTCAACGCCTTTCAGCTTGATGATCATGTCCTGATTGTGGTGCTTCCAGCTTTCCAGCTGGACTTGCTGCCACTGACCCGCCTTCTGGATAAACCAGGGCTGATAGTCAAAAATGCTTTCGGCGTCTTCGGTGGAAGAAAACACTC  >  NZ_CP009273/2739001‑2739189

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