Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,903,717 C→T 6.9% L212L (CTC→CTT)  manY → mannose‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,903,7170CT 6.9% 176.8 / 2.8 58L212L (CTC→CTT) manYmannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (32/22);  new base T (2/2);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

GCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTAACAACGAT  >  NC_000913/1903575‑1903861
                                                                                                                                              |                                                                                                                                                
gCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCtt                                                                                                                                             >  8:90326/1‑148 (MQ=255)
         aGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTGCTTCTACCTCGGCTTCGTAACCgc                                                                                                                                    <  2:353561/148‑1 (MQ=255)
          gAACATGCTGAATGCGATTCCGGAAGTGGTGCCCATTGGTCTGAATATCTCCGGTGGCATGATCGTGGTGGTTGGATATGCGATGTTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCTTAAcct                                                                                                                                     >  1:251861/1‑145 (MQ=255)
          gAACATGCTGAATGCGATTCCGGAAGTGGTGACGAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTTGGCtt                                                                                                                                             <  6:54120/138‑1 (MQ=255)
          gAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTTGGCtt                                                                                                                                             >  5:54120/1‑138 (MQ=255)
          gAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTGGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGGCGTTCTTCTACCTCGGCTTCGTAACCg                                                                                                                                     >  6:266339/1‑146 (MQ=255)
            aCATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCgcag                                                                                                                                  >  7:215556/1‑147 (MQ=255)
                     aTGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCaa                                                                                                                        >  7:233159/1‑148 (MQ=255)
                      tGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAAc                                                                                                                       >  7:37291/1‑148 (MQ=255)
                       gCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACt                                                                                                                      >  8:77630/1‑148 (MQ=255)
                               ggAAGTGGTGACCAATGGTCTGAATATCGCGGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTTGCTACCTGATGCCGTTCTTCTAACTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTg                                                                                                              <  2:216464/148‑1 (MQ=255)
                               ggAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTg                                                                                                              >  7:282237/1‑148 (MQ=255)
                                 aaGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGt                                                                                                            >  6:24710/1‑148 (MQ=255)
                                     ggtgACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTaa                                                                                                                  >  2:44642/1‑138 (MQ=255)
                                     ggtgACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTaa                                                                                                                  <  1:44642/138‑1 (MQ=255)
                                     ggtgACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGc                                                                                                         >  7:259214/1‑147 (MQ=255)
                                        gACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTaa                                                                                                                  >  7:320499/1‑135 (MQ=255)
                                        gACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTaa                                                                                                                  <  8:320499/135‑1 (MQ=255)
                                         aCCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTgg                                                                                                    >  8:178488/1‑148 (MQ=255)
                                          ccAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTaa                                                                                                                  >  6:434/1‑133 (MQ=255)
                                          ccAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTaa                                                                                                                  <  5:434/133‑1 (MQ=255)
                                           cAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGgt                                                                                                  >  3:194723/1‑148 (MQ=255)
                                            aaTGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTTGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGgtg                                                                                                 >  7:261248/1‑148 (MQ=255)
                                              tGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGtgtg                                                                                               <  4:378105/148‑1 (MQ=255)
                                                 tCTGAATATCGCCGGTGGCATGATCGTGGGGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGAtt                                                                                            <  6:179658/148‑1 (MQ=255)
                                                   tGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTgg                                                                                          >  1:90508/1‑148 (MQ=255)
                                                     aaTATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTa                                                                                        <  8:233159/148‑1 (MQ=255)
                                                           gCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGtt                                                                                   <  7:77630/147‑1 (MQ=255)
                                                                        tCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGt                                                                                         >  6:206960/1‑128 (MQ=255)
                                                                        tCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGt                                                                                         <  5:206960/128‑1 (MQ=255)
                                                                                          cGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGaa                                                   <  4:194723/148‑1 (MQ=255)
                                                                                          cGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGaa                                                   <  5:127362/148‑1 (MQ=255)
                                                                                            aTGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGaaa                                                  >  1:288325/1‑147 (MQ=255)
                                                                                             tGGATATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTACTCCCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGAGTGATTGGTACTGTTATGGCAATGCACTACATCCAACTTAGCCCGAAATa                                                >  4:152109/1‑148 (MQ=255)
                                                                                                  aTCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCCACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTg                                                                                                     <  6:143054/90‑1 (MQ=255)
                                                                                                  aTCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCCACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTg                                                                                                     >  5:143054/1‑90 (MQ=255)
                                                                                                      aCATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTTGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGt                                       <  8:261248/148‑1 (MQ=255)
                                                                                                      aCATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGt                                       <  5:24710/148‑1 (MQ=255)
                                                                                                       cATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATAAAACcgag                                        >  6:344192/1‑144 (MQ=255)
                                                                                                              cGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGGTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTg                               <  7:314759/148‑1 (MQ=255)
                                                                                                                    ggCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCt                                                                                                               <  4:44911/62‑1 (MQ=255)
                                                                                                                    ggCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCt                                                                                                               >  3:44911/1‑62 (MQ=255)
                                                                                                                     gCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGc                        >  4:296794/1‑148 (MQ=255)
                                                                                                                      cTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCCCCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCt                       >  7:312805/1‑148 (MQ=255)
                                                                                                                              tGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACcg                                          <  6:33168/121‑1 (MQ=255)
                                                                                                                              tGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACcg                                          <  6:27740/121‑1 (MQ=255)
                                                                                                                              tGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACcg                                          >  5:27740/1‑121 (MQ=255)
                                                                                                                              tGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACcg                                          >  5:33168/1‑121 (MQ=255)
                                                                                                                                ccGTTATTCTACCTCTGCTTCGTATCAGGAGCATTCACCAACTTTATCCTGGTTGCTGTGGGTGTGATGGGTAGGGTTATGGCCGTGCTCTACAGCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCt              <  2:251861/147‑1 (MQ=255)
                                                                                                                                  gTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACATTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTAAGCCCGAAATACAACCGCGTAGCCGCTGCGCCTGCTCAGGCAgaagg            >  7:164185/1‑143 (MQ=255)
                                                                                                                                   ttcttcTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGt           >  5:248926/1‑147 (MQ=255)
                                                                                                                                   ttcttcTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGt           >  5:254595/1‑147 (MQ=255)
                                                                                                                                    tcttctACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTa          >  8:171287/1‑147 (MQ=255)
                                                                                                                                     cttctACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTaa         >  8:58510/1‑147 (MQ=255)
                                                                                                                                        ctACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTaaca       <  2:288325/146‑1 (MQ=255)
                                                                                                                                           ccTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTAACAACGa   <  7:178488/147‑1 (MQ=255)
                                                                                                                                           ccTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTAACAACGa   <  3:296794/147‑1 (MQ=255)
                                                                                                                                           ccTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTAACAACGAt  >  8:240509/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCATGATCGTGGTGGTTGGTTATGCGATGGTTATCAACATGATGCGTGCTGGCTACCTGATGCCGTTCTTCTACCTCGGCTTCGTAACCGCAGCATTCACCAACTTTAACCTGGTTGCTCTGGGTGTGATTGGTACTGTTATGGCAGTGCTCTACATCCAACTTAGCCCGAAATACAACCGCGTAGCCGGTGCGCCTGCTCAGGCAGCTGGTAACAACGAT  >  NC_000913/1903575‑1903861

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