Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 983,587 T→A 7.0% intergenic (+45/‑44) iscX → / → ndk cysteine desulfurase activity regulator/nucleoside diphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947983,5870TA 7.0% 173.9 / 3.0 57intergenic (+45/‑44)iscX/ndkcysteine desulfurase activity regulator/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base A (2/2);  total (31/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.36e-01

GTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCC  >  NC_002947/983447‑983724
                                                                                                                                            |                                                                                                                                         
gTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGtt                                                                                                                                   <  1:155382/149‑1 (MQ=255)
           gCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCggg                                                                                                                                     <  3:191279/136‑1 (MQ=255)
           gCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCggg                                                                                                                                     >  4:191279/1‑136 (MQ=255)
            ccAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGt                                                                                                                        >  8:211896/1‑148 (MQ=255)
             cAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTaa                                                                                                                      >  6:6006/1‑149 (MQ=255)
                agTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAacac                                                                                                                   <  2:113303/149‑1 (MQ=255)
                  ttCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCggg                                                                                                                                     <  6:257325/129‑1 (MQ=255)
                  ttCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCggg                                                                                                                                     >  5:257325/1‑129 (MQ=255)
                   tCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCa                                                                                                                >  5:34399/1‑149 (MQ=255)
                       cgacgaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATAtt                                                                                                            <  3:116277/149‑1 (MQ=255)
                       cgacgaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATAtt                                                                                                            >  7:266400/1‑149 (MQ=255)
                            aTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                       <  5:6006/149‑1 (MQ=255)
                            aTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                       <  4:54036/149‑1 (MQ=255)
                             tCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAAt                                                                                                                               >  8:258102/1‑124 (MQ=255)
                                   aCGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTc                                                                                                 <  6:34399/148‑1 (MQ=255)
                                   aCGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCAGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTc                                                                                                 >  3:252960/1‑148 (MQ=255)
                                    cgcgGCGGTGAGAACCTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGCAGCCGACTAAGCGTCACTTGCAGGTTAGCCAATCCCCCGGACCCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACCCTCATATTCTGGAGTTTTCCa                                                                                               <  4:252960/149‑1 (MQ=255)
                                     gcggcgGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTgg                                                                                                        <  4:32965/139‑1 (MQ=255)
                                     gcggcgGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTgg                                                                                                        >  3:32965/1‑139 (MQ=255)
                                     gcggcgGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCAt                                                                                              <  3:229084/149‑1 (MQ=255)
                                       ggcggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATgg                                                                                            >  4:199349/1‑149 (MQ=255)
                                          ggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCa                                                                                               <  4:145845/143‑1 (MQ=255)
                                          ggTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCa                                                                                               >  3:145845/1‑143 (MQ=255)
                                              agaAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTTCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCa                                                                                     <  8:266155/149‑1 (MQ=255)
                                                aaaGTGCTCGAGGAAATCCAGGACGCCTGGATCGACGAAGCCGAATAAGCGTCACTTGCAGGTTAGGCAATCCCGCGGAACCCGATTATAATTCGAGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGCTCAAc                                                                                   >  2:200347/1‑149 (MQ=255)
                                                   gTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTc                                                                                                 >  6:235287/1‑132 (MQ=255)
                                                   gTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTc                                                                                                 <  5:235287/132‑1 (MQ=255)
                                                    tGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTAc                                                                               <  2:297213/149‑1 (MQ=255)
                                                    tGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTAc                                                                               >  5:127052/1‑149 (MQ=255)
                                                        cGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                       <  5:123343/121‑1 (MQ=255)
                                                        cGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                       >  6:123343/1‑121 (MQ=255)
                                                              aaTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcat                                                                      >  7:124193/1‑148 (MQ=255)
                                                              aaTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcat                                                                      >  7:119952/1‑148 (MQ=255)
                                                                         ccTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcatc                                                                     >  3:23792/1‑138 (MQ=255)
                                                                         ccTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcatc                                                                     <  4:23792/138‑1 (MQ=255)
                                                                                 gacgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGtt                                                   <  8:266400/148‑1 (MQ=255)
                                                                                  acgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                 <  7:211896/149‑1 (MQ=255)
                                                                                  acgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                 >  7:73294/1‑149 (MQ=255)
                                                                                  acgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGc                                                 <  4:73151/149‑1 (MQ=255)
                                                                                      aGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCGATCATCAAGCCTGACGCTGTTGCCaa                                              >  5:115647/1‑148 (MQ=255)
                                                                                          gACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACg                                         >  3:58000/1‑149 (MQ=255)
                                                                                               aGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGcc                                                             <  4:151717/124‑1 (MQ=255)
                                                                                               aGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGcc                                                             >  3:151717/1‑124 (MQ=255)
                                                                                                gCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCg                                   >  8:247804/1‑149 (MQ=255)
                                                                                                gCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCg                                   >  8:295094/1‑149 (MQ=255)
                                                                                                gCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCg                                   >  7:150069/1‑149 (MQ=255)
                                                                                                       ttGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTATTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGa                             >  1:68478/1‑148 (MQ=255)
                                                                                                         gCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATca                          >  6:181897/1‑149 (MQ=255)
                                                                                                         gCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCACCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATca                          >  4:110915/1‑149 (MQ=255)
                                                                                                                         cccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaag           >  8:157259/1‑148 (MQ=255)
                                                                                                                          ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCgaagaa         <  7:295094/149‑1 (MQ=255)
                                                                                                                              aaCCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCg     <  8:73294/149‑1 (MQ=255)
                                                                                                                               aCCCGCGTATAATTCGCGTGTTTAATTTTTCTTAACACTCATTTTCTGGAGTTATCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCgg    <  6:115647/149‑1 (MQ=255)
                                                                                                                                cccGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGc   <  5:181897/149‑1 (MQ=255)
                                                                                                                                 ccGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCAGCACTCGCTTCGAAGAAGCCGGcc  >  7:120167/1‑149 (MQ=255)
                                                                                                                                   gcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATc                           >  8:115175/1‑122 (MQ=255)
                                                                                                                                   gcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATc                           <  7:115175/122‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                         
GTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCGGCC  >  NC_002947/983447‑983724

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