Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,924,518 G→C 9.1% T99S (ACC→AGC) ‡ peaB ← quinohemoprotein amine dehydrogenase modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,924,5180GC 9.1% 170.9 / 6.2 55T99S (ACC→AGC) ‡peaBquinohemoprotein amine dehydrogenase modification protein
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base C (2/3);  total (25/30)
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.68e-01

CAGCATCTCCACCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAG  >  NC_002947/3924379‑3924651
                                                                                                                                           |                                                                                                                                      
cagcaTCTCCACCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCACCGCGCTcaa                                                                                                                               <  5:358980/149‑3 (MQ=255)
         cACCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCt                                                                                                                       >  6:253988/1‑147 (MQ=255)
         cACCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTACCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGc                                                                                                                         <  2:105632/145‑1 (MQ=255)
          aCCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCt                                                                                                                       <  8:170657/146‑1 (MQ=255)
          aCCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCt                                                                                                                       >  7:170657/1‑146 (MQ=255)
          aCCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAg                                                                                                                    <  4:297857/149‑1 (MQ=255)
           ccGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAgt                                                                                                                   <  7:86466/149‑2 (MQ=255)
                gCTTCGGCGGTGGTGGTACTCATCATCTTCCCGGCGGAAGGCTTGTACAGGACTTCCTTCTAGCAATAGGTGAAGCTGAGGTTGCAAACGGTGTTGACGTTCAGCACCACGCAGTTGAGCGCGGTGAGCTCGACTTGCCTGAGGGCGAt                                                                                                              >  8:99936/1‑149 (MQ=255)
                     ggcggTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCg                                                                                                         <  3:65816/149‑1 (MQ=255)
                     ggcggTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCg                                                                                                         <  3:289134/149‑1 (MQ=255)
                       aggtggttgtACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGGGCTCGACTTGCCTGAGGGCGATTTCCggg                                                                                                       <  7:99936/148‑1 (MQ=255)
                        ggtggtggtACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGAACTTGCCTGAGGGCGATTTCCggg                                                                                                       <  7:48300/149‑1 (MQ=255)
                        ggtggtggtACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCgggg                                                                                                      >  3:257069/1‑149 (MQ=255)
                                 aCTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAgcggcg                                                                                             >  6:163117/1‑149 (MQ=255)
                                   tCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTCGAGCGCACCGCGCTCAAcaccgtg                                                                                                                       >  7:307187/1‑115 (MQ=255)
                                   tCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTCGAGCGCACCGCGCTCAAcaccgtg                                                                                                                       <  8:307187/121‑7 (MQ=255)
                                    cATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCACCGCGCTCAAc                                                                                                                             >  6:331316/1‑114 (MQ=255)
                                    cATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCACCGCGCTCAAc                                                                                                                             <  5:331316/114‑1 (MQ=255)
                                         tcttGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCt                                                                                                                                  >  5:60817/1‑105 (MQ=255)
                                         tcttGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCt                                                                                                                                  <  6:60817/105‑1 (MQ=255)
                                              ccGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGACGGCTCGCCGCTGa                                                                                 <  5:253988/148‑1 (MQ=255)
                                               cggcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCg                                                                                                                >  2:290939/1‑116 (MQ=255)
                                               cggcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCg                                                                                                                <  1:290939/116‑1 (MQ=255)
                                                ggcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCggg                                                                                                       >  1:2804/1‑124 (MQ=255)
                                                ggcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCggg                                                                                                       <  2:2804/124‑1 (MQ=255)
                                                ggcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGACGGCTCGCCGCTGAccc                                                                              <  6:4741/149‑2 (MQ=255)
                                                       ggCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAg                                                                       <  5:21133/149‑1 (MQ=255)
                                                       ggCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAg                                                                       >  6:21133/1‑149 (MQ=255)
                                                             tCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGa                                                                  <  4:257069/148‑1 (MQ=255)
                                                                  gTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTc                                                                          >  4:117672/1‑135 (MQ=255)
                                                                  gTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTc                                                                          <  3:117672/135‑1 (MQ=255)
                                                                    cTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTcg                                                          >  4:4281/1‑149 (MQ=255)
                                                                       ccTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCg                                                                                                         <  3:256408/99‑1 (MQ=255)
                                                                       ccTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCg                                                                                                         >  4:256408/1‑99 (MQ=255)
                                                                        cTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTcc                                                                         <  4:299819/130‑1 (MQ=255)
                                                                        cTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTcc                                                                         >  3:299819/1‑130 (MQ=255)
                                                                                   aGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTc                                                           >  1:197624/1‑133 (MQ=255)
                                                                                   aGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTc                                                           <  2:197624/133‑1 (MQ=255)
                                                                                     gTGCAGCTGAGGTTGCAACCGGTGTTGAGGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGc                                                                            >  1:196251/1‑114 (MQ=255)
                                                                                     gTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGc                                                                            <  2:196251/114‑1 (MQ=255)
                                                                                      tGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGata                                        >  3:42021/1‑149 (MQ=255)
                                                                                        cAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATc                                      <  4:42021/149‑1 (MQ=255)
                                                                                             gAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGc                                 >  5:250996/1‑149 (MQ=255)
                                                                                                       ccGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActg                       >  6:240989/1‑149 (MQ=255)
                                                                                                        cGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActgc                      <  5:240989/149‑1 (MQ=255)
                                                                                                          gtgtTGACGTTGAGCACCACGGTGTTGAGCGCTGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActgctg                    <  6:250996/149‑1 (MQ=255)
                                                                                                            gtTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGc                  >  8:334748/1‑149 (MQ=255)
                                                                                                                    tGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGc                                 <  7:204791/126‑1 (MQ=255)
                                                                                                                    tGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGc                                 >  8:204791/1‑126 (MQ=255)
                                                                                                                      aGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTg        <  3:4281/149‑1 (MQ=255)
                                                                                                                        caccacGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGGGATCAGTTCGAGCAGGGTGTCCTGGATATc                                      <  1:32485/117‑1 (MQ=255)
                                                                                                                        caccacGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGAGCAGGGTGTCCTGGATATc                                      >  2:32485/1‑117 (MQ=255)
                                                                                                                           cacGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCAt     >  2:6331/1‑147 (MQ=255)
                                                                                                                            acGGTGTTGAGCGCGGTGCGCTCGAATTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAg  <  3:145018/149‑1 (MQ=255)
                                                                                                                                    gAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGcc                 >  2:261799/1‑126 (MQ=255)
                                                                                                                                           |                                                                                                                                      
CAGCATCTCCACCGAGGCTTCGGCGGTGGTGGTACTCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAG  >  NC_002947/3924379‑3924651

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