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 I2 R1 882 254.4 12957588 98.0% 12698436 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,941,738 2 bp→GG noncoding (413‑414/2906 nt) rrlC → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,941,7380AG100.0% 128.9 / NA 65noncoding (413/2906 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/26);  total (39/26)
*NC_0009133,941,7390CG90.6% 93.2 / 8.6 64noncoding (414/2906 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base G (34/24);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGAGGGAC‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAA  >  NC_000913/3941648‑3941797
                                                                                            |||                                                            
cAGGGTGACAGCCCCGTACTCAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGc‑‑c                                                           >  2:6273733/1‑93 (MQ=255)
         aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCa                                                   >  1:3116169/1‑93 (MQ=11)
          gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCaa                                                  >  1:5712588/1‑93 (MQ=14)
          gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGC‑‑CATCCTCCAAg                                                 >  2:1588095/1‑93 (MQ=255)
           ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGG‑‑CATCCTCCAAgg                                                >  1:1445065/1‑93 (MQ=255)
            cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAgg                                                >  2:2155459/1‑93 (MQ=14)
            cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CAGCCTCCAAgg                                                >  2:311875/1‑93 (MQ=11)
             cccaacacAAAAATGCACATATTGTGAGCCCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGGGCATCCTCCAAGGc                                               <  2:6075164/89‑1 (MQ=14)
             ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CCTCCTCCAAGGc                                               >  2:1473406/1‑93 (MQ=11)
             ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAggt                                               >  1:3913488/1‑92 (MQ=14)
             ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGc                                               >  2:931274/1‑93 (MQ=14)
              cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCt                                              >  2:4858622/1‑93 (MQ=14)
              cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCt                                              >  1:2460808/1‑93 (MQ=14)
              cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCt                                              >  1:645120/1‑93 (MQ=14)
               gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGGG‑‑CATCCTCCAAGGCt                                              >  1:2038145/1‑93 (MQ=14)
               gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTa                                             >  2:2122032/1‑93 (MQ=14)
                tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTaa                                            >  2:4749652/1‑93 (MQ=14)
                tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTaa                                            >  1:5211382/1‑93 (MQ=14)
                tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGC‑‑CATCCTCCAAGGCTaaa                                           >  1:839205/1‑93 (MQ=255)
                 acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑GATCCTCCAAGGCTaaa                                           >  1:1154549/1‑93 (MQ=14)
                  cgcaAAAATGCGCATATTGTGAGCTCGATGGGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAAt                                          <  2:5922525/91‑1 (MQ=14)
                  cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAAt                                          <  1:1031130/93‑1 (MQ=11)
                  cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAAt                                          <  2:4814714/93‑1 (MQ=11)
                  cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAAt                                          >  2:3400246/1‑93 (MQ=14)
                  cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAAt                                          >  1:866125/1‑93 (MQ=14)
                   acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATa                                         <  2:4090950/93‑1 (MQ=11)
                   acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATa                                         >  1:4006111/1‑93 (MQ=14)
                   acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATa                                         <  1:4749652/93‑1 (MQ=11)
                   aaaaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATa                                         <  1:4400753/91‑1 (MQ=11)
                    caAAAATGCACATATTGTGAGCTCGATGAGTAGGGGGGGACCCGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGGG‑‑CATCCTCCAAGGCTAAATAc                                        <  1:83801/93‑1 (MQ=11)
                    caAAAATGCACATATTGTGAGCTCGATGAGTAGGGGGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGGGGCATCCTCCAAGGCTAAATa                                         <  1:2580697/93‑1 (MQ=14)
                    caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATAc                                        >  1:3702383/1‑93 (MQ=14)
                    caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATAc                                        <  1:1430736/93‑1 (MQ=11)
                    caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGGG‑‑CATCCTCCAAGGCTAAATAc                                        <  2:3142469/93‑1 (MQ=14)
                     aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACt                                       >  1:2787499/1‑93 (MQ=14)
                     aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCCAAATACt                                       >  2:942568/1‑93 (MQ=14)
                       aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTcc                                     <  1:5279350/93‑1 (MQ=11)
                       aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTcc                                     <  2:3210058/93‑1 (MQ=11)
                       aaaTGCACATATGGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCct                                    <  2:376092/93‑1 (MQ=14)
                        aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCct                                    <  1:6309823/93‑1 (MQ=11)
                        aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCct                                    >  2:69703/1‑93 (MQ=14)
                        aaTGCACATATTGGGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTTTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCct                                    <  1:1126002/93‑1 (MQ=11)
                         aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctg                                   >  2:1066217/1‑93 (MQ=14)
                         aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctg                                   >  1:1749477/1‑93 (MQ=14)
                          tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctga                                  >  2:1963151/1‑93 (MQ=14)
                          tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCCAAATACTCctga                                  >  1:163937/1‑93 (MQ=14)
                           gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgac                                 <  1:5881249/93‑1 (MQ=11)
                           gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgac                                 <  1:1059644/93‑1 (MQ=11)
                           gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgac                                 >  2:2055587/1‑93 (MQ=14)
                           gCAAATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGGGCATCCTCCAAGGCTAAATACTCctgac                                 <  1:3411698/93‑1 (MQ=11)
                            cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgact                                >  1:4647437/1‑93 (MQ=14)
                            cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgact                                <  1:1646951/93‑1 (MQ=14)
                             acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgactg                               >  2:2803057/1‑93 (MQ=14)
                             acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgactg                               >  2:231021/1‑93 (MQ=14)
                             acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgactg                               >  1:1616296/1‑93 (MQ=14)
                             acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgactg                               >  2:5906840/1‑93 (MQ=14)
                             acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgactg                               >  2:5413094/1‑93 (MQ=14)
                              caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑GATCCCCCAAGGCTAAATACTCccgactga                              >  2:380014/1‑93 (MQ=11)
                              caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCctgactga                              >  2:804026/1‑93 (MQ=11)
                               atatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGG‑‑GATCCTCCAAGGCTAAATACTCCTGACTCAcc                            >  2:5649042/1‑93 (MQ=14)
                                tatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑GATCCTCCAAGGCTAAATACTCCTGACTGAcc                            >  1:310321/1‑93 (MQ=11)
                                tatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCCTGACTGAcc                            <  1:5802105/93‑1 (MQ=11)
                                tatTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGC‑‑CATCCTCCAAGGCTTAATACTCCTGACTGACCg                           >  1:3999151/1‑93 (MQ=255)
                                 gttGTGAGCTCGAGGAGTGGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCg                           <  1:4797222/92‑1 (MQ=11)
                                 atTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCg                           >  1:4313942/1‑93 (MQ=11)
                                  ttGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGa                          >  1:1974852/1‑93 (MQ=11)
                                  ttGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGa                          <  1:3443664/93‑1 (MQ=11)
                                  ttGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGa                          >  1:410944/1‑93 (MQ=11)
                                     tgAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGG‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGt                      <  1:3111072/93‑1 (MQ=255)
                                                      gCGGGACACGTGGTCTCCTGTCTGAATAT‑‑GGGGGGGGC‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTggag     >  2:4683291/1‑90 (MQ=11)
                                                        gggACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGC‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGa   <  1:4610222/93‑1 (MQ=16)
                                                         ggACACGTGGTATCCTGTCTGAATAT‑‑GGGGGGGGC‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGaa  <  2:4932633/93‑1 (MQ=16)
                                                                                            |||                                                            
CAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATAT‑‑GGGAGGGAC‑‑CATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAA  >  NC_000913/3941648‑3941797

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

GATK/CNVnator alignment

N/A