Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,831,004 C→T 6.1% G59G (GGC→GGT U712_09105 → Uncharacterized protein ymzC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,831,0040CT 6.1% 207.3 / 3.2 66G59G (GGC→GGTU712_09105Uncharacterized protein ymzC
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base T (2/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CACTTGTATGGATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGATGCCAGCCGCT  >  CP006881/1830868‑1831142
                                                                                                                                        |                                                                                                                                          
cACTTGTATGGATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGt                                                                                                                                          >  8:87925/1‑139 (MQ=255)
cACTTGTATGGATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGt                                                                                                                                          >  2:380848/1‑139 (MQ=255)
     gTATGGATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGACCATAACCATTTTTGCGTTATgg                                                                                                                                     >  6:285732/1‑139 (MQ=255)
     gTATGGATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATg                                                                                                                                      <  6:141930/138‑1 (MQ=255)
     gTATGGATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATg                                                                                                                                      <  2:32200/138‑1 (MQ=255)
          gATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGAc                                                                                                                                >  5:224533/1‑139 (MQ=255)
           aTAGCTGTCATTTTACTGTTCGGGGCGTGCGGGAATCAAGAAACCATTATTGAAACCGACAACGCCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACGTGAACATAACCATTTTGGCGTTATGGAGGACg                                                                                                                               >  8:89042/1‑139 (MQ=255)
             aGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGAc                                                                                                                                >  3:284160/1‑136 (MQ=255)
               cTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGTGTTATggag                                                                                                                                   >  7:34461/1‑131 (MQ=255)
               cTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGTGTTATggag                                                                                                                                   <  8:34461/131‑1 (MQ=255)
               cTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCta                                                                                                                           >  6:115829/1‑139 (MQ=255)
               cTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAAAGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCta                                                                                                                           >  2:104115/1‑139 (MQ=255)
                    tttttACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATc                                                                                                                       >  1:25456/1‑138 (MQ=255)
                            attCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTa                                                                                                               <  8:93660/137‑1 (MQ=255)
                            attCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCAACCTCGTTTCCACTTGAACATAACCATTTTTGCGTTATGGAGGACGGCTATATCAAAATTTa                                                                                                               >  7:93660/2‑138 (MQ=255)
                             ttCGGGGCGTGCGGGAATCAAGATACCATTATTGAAATAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTa                                                                                                               <  5:7759/137‑1 (MQ=255)
                             ttCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAAtt                                                                                                                 <  4:342354/135‑1 (MQ=255)
                             ttCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAAtt                                                                                                                 <  4:258275/135‑1 (MQ=255)
                             ttCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTa                                                                                                               <  8:138527/137‑1 (MQ=255)
                             ttCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTa                                                                                                               <  1:95550/137‑1 (MQ=255)
                             ttCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTAt                                                                                                              <  6:117505/138‑1 (MQ=255)
                              tCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGa                                                                                                            >  6:303026/1‑139 (MQ=255)
                                     gTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGtat                                                                                                        <  4:284160/136‑1 (MQ=255)
                                       gCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACgg                                                                                                                              <  7:326865/112‑1 (MQ=255)
                                       gCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACgg                                                                                                                              >  8:326865/1‑112 (MQ=255)
                                       gCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAAt                                                                                                     >  6:395980/1‑137 (MQ=255)
                                                aataTACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCg                                                                                             <  3:369250/133‑1 (MQ=255)
                                                aaGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCaa                                                                                          <  1:348725/139‑1 (MQ=255)
                                                aaGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCaa                                                                                          <  2:67414/139‑1 (MQ=255)
                                                aaGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCaa                                                                                          <  1:380848/139‑1 (MQ=255)
                                                     aCCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCaa                                                                                          <  6:33722/134‑1 (MQ=255)
                                                     aCCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCaa                                                                                          >  5:33722/1‑134 (MQ=255)
                                                       cATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAg                                                                                      <  3:150213/136‑1 (MQ=255)
                                                            ttGAAACAGACAACGGCAATTCAGACTATGAAATGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGc                                                                               <  2:30293/138‑1 (MQ=255)
                                                                 aCAGACAACGGCAATTCAGACTATGAACCGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTgaag                                                                          >  2:391242/1‑138 (MQ=255)
                                                                    gACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGTGTTATGGAGGACGGCTATATCAAAATTTAt                                                                                                              <  8:380476/99‑1 (MQ=255)
                                                                    gACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGTGTTATGGAGGACGGCTATATCAAAATTTAt                                                                                                              >  7:380476/1‑99 (MQ=255)
                                                                     aCAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGaagaaag                                                                      >  1:320906/1‑138 (MQ=255)
                                                                               ttCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGgat                                                           >  1:305099/1‑139 (MQ=255)
                                                                                    aCTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGgatgat                                                        >  7:80544/1‑137 (MQ=255)
                                                                                       aTGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGgatgat                                                        >  3:354483/1‑134 (MQ=255)
                                                                                             cGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGTCGTTATGGAGGCCGGATATATAAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGATGAATAAAGAATACGCGGATGATGAGCTTGAATa                                             <  2:21067/139‑1 (MQ=255)
                                                                                                cTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATaaaa                                          >  2:316957/1‑139 (MQ=255)
                                                                                                 tCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAAtt                                                                                                                 >  3:246483/1‑67 (MQ=255)
                                                                                                 tCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAAtt                                                                                                                 <  4:246483/67‑1 (MQ=255)
                                                                                                      ccACCTCGTTTCCACTTTAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTcaca                                    <  5:285732/139‑1 (MQ=255)
                                                                                                         ccTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAg                                  >  2:27584/1‑138 (MQ=255)
                                                                                                         ccTCGTTTCAACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGAGGATGATGAGCTTGAATAAAATTCACAAg                                  >  6:174532/1‑138 (MQ=255)
                                                                                                          cTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGTGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCa                                >  2:278105/1‑139 (MQ=255)
                                                                                                          cTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTAt                                                                                                              <  1:163167/61‑1 (MQ=255)
                                                                                                          cTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTAt                                                                                                              <  1:125386/61‑1 (MQ=255)
                                                                                                          cTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTAt                                                                                                              >  2:125386/1‑61 (MQ=255)
                                                                                                          cTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTAt                                                                                                              >  2:163167/1‑61 (MQ=255)
                                                                                                          cTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCa                                >  7:305416/1‑139 (MQ=255)
                                                                                                             gTTTCCACTTGAACGTAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCAtt                              <  7:98874/138‑1 (MQ=255)
                                                                                                               ttCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCg                            >  7:16334/1‑138 (MQ=255)
                                                                                                                    cTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGCGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTa                      >  7:291337/1‑139 (MQ=255)
                                                                                                                       gAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACa                    <  2:320906/138‑1 (MQ=255)
                                                                                                                       gAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACa                    >  1:385045/1‑138 (MQ=255)
                                                                                                                          cATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACa                    >  4:402156/1‑135 (MQ=255)
                                                                                                                             aaCCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGa             >  8:102176/1‑139 (MQ=255)
                                                                                                                             aaCCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGa             >  6:377963/1‑139 (MQ=255)
                                                                                                                                 aTTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGATGcc         <  2:305099/139‑1 (MQ=255)
                                                                                                                                  ttttGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGATg           >  4:217505/1‑136 (MQ=255)
                                                                                                                                  ttttGGAGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCTAATTTACAAAGCGGATGCCa        >  7:74566/1‑139 (MQ=255)
                                                                                                                                   tttGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGATGCCAg       >  8:390577/1‑139 (MQ=255)
                                                                                                                                        cGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGATGCCAGCCGct  <  7:87925/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                          
CACTTGTATGGATAGCTGTCTTTTTACTGTTCGGGGCGTGCGGGAATCAAGATACCATTATTGAAACAGACAACGGCAATTCAGACTATGAAACGCCTCAGCCCACCTCGTTTCCACTTGAACATAACCATTTTGGCGTTATGGAGGACGGCTATATCAAAATTTATGAGTATAATGAGTCCCGCAATGAGGTAAAGCTGAAGAAAGAATACGCGGATGATGAGCTTGAATAAAATTCACAAGCATTCGAATTTACAAAGCGGATGCCAGCCGCT  >  CP006881/1830868‑1831142

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