Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 132,471 T→C 9.4% intergenic (+45/‑107) cc → / → dsbA cytochrome c4/thiol:disulfide interchange protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947132,4710TC 9.4% 142.0 / 5.3 53intergenic (+45/‑107)cc/dsbAcytochrome c4/thiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base C (3/2);  total (27/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 = 6.00e-01

ACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGCCGCCAGC  >  NC_002947/132328‑132619
                                                                                                                                               |                                                                                                                                                    
aCCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCt                                                                                                                                                 >  7:23212/1‑149 (MQ=255)
   aaCGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTCGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCtttt                                                                                                                                              <  6:345486/149‑1 (MQ=255)
       aTGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttcatca                                                                                                                                          <  5:19718/149‑6 (MQ=255)
         ggCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgt                                                                                                                                        >  8:145716/1‑149 (MQ=255)
         ggCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgt                                                                                                                                        >  5:203887/1‑149 (MQ=255)
          gCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc                                                                                                                                       <  4:352391/149‑1 (MQ=255)
                 cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGGTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGACCTTTTTTCCGTCCGCAGcc                                                                                                                                >  4:341044/1‑149 (MQ=255)
                 cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGGTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCGTTTTTTCCGTCCGCAGcc                                                                                                                                >  4:42645/1‑149 (MQ=255)
                      ccATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTac                                                                                                                           <  6:6708/149‑1 (MQ=255)
                       cATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTaca                                                                                                                          >  7:214508/1‑149 (MQ=255)
                              aCCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTACCGTCCGCAGCCGTTACACTACAg                                                                                                                    >  7:319466/1‑148 (MQ=255)
                               ccATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGaa                                                                                                                  <  6:284239/149‑1 (MQ=255)
                                        gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc                                                                                                                                       >  2:338962/1‑119 (MQ=255)
                                        gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc                                                                                                                                       <  1:338962/119‑1 (MQ=255)
                                        gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGc                                                                                                           <  8:264408/147‑1 (MQ=255)
                                        gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGc                                                                                                           <  6:47684/147‑1 (MQ=255)
                                        gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttc                                                                                                                                              <  7:325336/112‑2 (MQ=255)
                                        gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttc                                                                                                                                              >  8:325336/1‑111 (MQ=255)
                                                aGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCg                                                                                                   <  6:19104/147‑1 (MQ=255)
                                                  cAACCATGACATCGAAGCGTAGGCCAGCTACATCCAGGGCCTGCAATAACTTTCGGTTAAAGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTCCCCTTTTTTCCATACGCAGCCGTTAAACTACAGAACACCAGCCCTTCCAGGcat                                                                                               >  3:5076/1‑147 (MQ=255)
                                                       aTGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtct                                                                                           >  1:99721/1‑148 (MQ=255)
                                                        tGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGa                                                                                                             >  2:201197/1‑129 (MQ=255)
                                                        tGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCCTCCGCAGCCGTTACACTACAGAACTCGa                                                                                                             <  1:201197/129‑1 (MQ=255)
                                                          aCATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGt                                                                                       >  1:94804/1‑149 (MQ=255)
                                                            aTCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGc                                                                                     <  3:243940/149‑1 (MQ=255)
                                                             tCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTg                                                                                      >  3:333093/1‑147 (MQ=255)
                                                               gAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAAGTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCAGCAGCCGTTACACTACAGAACTCGCGCCATTCCCGCCCTTTCTCAGTGCAgg                                                                                  >  3:264066/1‑149 (MQ=255)
                                                               gAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAgg                                                                                  >  6:164306/1‑149 (MQ=255)
                                                                 aGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTc                                                                                <  6:203887/149‑1 (MQ=255)
                                                                        ggaaGCTTCATCCAGGGGCTGCATTAACTTTCGTTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGcgtcgt                                                                         <  4:5076/146‑1 (MQ=255)
                                                                                cATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAcc                                                                  >  1:216409/1‑148 (MQ=255)
                                                                                  tCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCAt                                                               <  2:94804/149‑1 (MQ=255)
                                                                                             cATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg                                                     >  7:107245/1‑148 (MQ=255)
                                                                                             cATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGATACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg                                                     >  3:142440/1‑148 (MQ=255)
                                                                                                 aaCTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTACCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGt                                                 <  8:319466/148‑1 (MQ=255)
                                                                                                   cTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTaaa                                              <  8:214508/149‑1 (MQ=255)
                                                                                                    tttCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCGTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAg                                             >  5:265503/1‑149 (MQ=255)
                                                                                                       cGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCa                                           <  7:145716/148‑1 (MQ=255)
                                                                                                            aaTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGt                                      >  4:264569/1‑148 (MQ=255)
                                                                                                             aTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTaa                                    <  6:265503/149‑1 (MQ=255)
                                                                                                             aTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTa                                     >  1:296341/1‑148 (MQ=255)
                                                                                                                   tAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGAt                              >  3:274620/1‑149 (MQ=255)
                                                                                                                       tCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATtctc                          <  2:99721/149‑1 (MQ=255)
                                                                                                                             aaaaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgct                    >  4:143020/4‑149 (MQ=255)
                                                                                                                              aaaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctg                   >  4:44626/3‑149 (MQ=255)
                                                                                                                               gaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctgc                  >  7:40126/1‑149 (MQ=255)
                                                                                                                                   gggCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGcgctgcgctg              <  8:40126/149‑1 (MQ=255)
                                                                                                                                    ggcggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTgg             <  2:296341/149‑1 (MQ=255)
                                                                                                                                       ggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAAc                                  >  4:240061/1‑125 (MQ=255)
                                                                                                                                       ggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAAc                                  <  3:240061/125‑1 (MQ=255)
                                                                                                                                           gggCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTcgccgc      <  8:23212/149‑1 (MQ=255)
                                                                                                                                             tcTGCCCTTTTTTTCGTCCGTAGCCGTTAAACTACAGAACTGGATCCCTTCCCGGCCTGTCTCAGTGCAGGTCGAGTCGTGGCGACCCATTACTGATCAGGAGTAAAGCATGCGTTAACTGATTCTCAGGGCTGCGCTGGTCGCCGCCa    <  4:264066/148‑1 (MQ=255)
                                                                                                                                             gCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTATCAGTGCTGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTTATTCTGAGCGCTGCGCTGGTCGCCGCCa    <  4:142440/149‑1 (MQ=255)
                                                                                                                                               tGCCCTTTTTTCCGTACGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGCCGCCAGc  >  2:80370/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGGTCGCCGCCAGC  >  NC_002947/132328‑132619

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: ‑