Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,824,139 G→T 7.5% A359S (GCC→TCC)  yggW → coproporphyrinogen/heterocyclic compound oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,824,1390GT 7.5% 161.5 / 3.7 53A359S (GCC→TCC) yggWcoproporphyrinogen/heterocyclic compound oxidase
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base T (2/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.89e-01

AGGGCCGGCGAGAAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATGGATCTGGTACTTGATCTGCTCGCGACG  >  NC_002947/5823992‑5824279
                                                                                                                                                   |                                                                                                                                            
aGGGCCGGCGAGAAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTgccgc                                                                                                                                             <  2:35491/149‑1 (MQ=255)
    ccGGCGAGAAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCTCAGCTACAGGAAGCACGCCGTGCCGCCGaa                                                                                                                                         <  1:23278/149‑1 (MQ=255)
    ccGGCGAGAAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGaa                                                                                                                                         >  2:207117/1‑149 (MQ=255)
     cGGCGAGAAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAAc                                                                                                                                        <  1:207117/149‑1 (MQ=255)
       gCGAGAAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACa                                                                                                                                       >  3:74807/1‑148 (MQ=255)
           gaAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAgg                                                                                                                                  >  7:25589/1‑149 (MQ=255)
            aaGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCTCCGAACAAAAggg                                                                                                                                 >  7:7963/1‑149 (MQ=255)
            aaGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCTCCGAACAAAAggg                                                                                                                                 >  7:11901/1‑149 (MQ=255)
             aGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCTCCGAACAAAAGGGc                                                                                                                                <  8:7963/149‑1 (MQ=255)
             aGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCTCCGAACAAAAGGGc                                                                                                                                <  8:11901/149‑1 (MQ=255)
                     ccAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTa                                                                                                                          >  2:217702/1‑147 (MQ=255)
                      cAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAg                                                                                                                       <  8:164299/149‑1 (MQ=255)
                             tGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAAc                                                                                                                <  1:159389/149‑1 (MQ=255)
                                    ccGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGaac                                                                                                          >  4:20316/1‑146 (MQ=255)
                                      gTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGAc                                                                                                       <  3:195859/149‑1 (MQ=255)
                                          gAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTgg                                                                                                    >  2:118544/1‑148 (MQ=255)
                                          gAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTgg                                                                                                    >  6:224235/1‑148 (MQ=255)
                                                 tgatgaACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACGGGTAGCCACg                                                                                            >  8:174856/1‑149 (MQ=255)
                                                   atgaACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGcc                                                                                          >  6:111941/1‑149 (MQ=255)
                                                        cGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGACGCGGATCTCCTGGTAGCCACGCCGAggg                                                                                     >  8:186455/1‑149 (MQ=255)
                                                        cGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGACCGACTGGTAGCCACGCCGAggg                                                                                     <  3:20316/149‑1 (MQ=255)
                                                           ccTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAgg                                                                                      <  8:77215/145‑1 (MQ=255)
                                                           ccTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAgg                                                                                      >  7:77215/1‑145 (MQ=255)
                                                               cgcCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAgtt                                                                              >  5:22000/1‑149 (MQ=255)
                                                                  cTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAgttgtt                                                                           <  3:233346/149‑1 (MQ=255)
                                                                  cTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAgttgtt                                                                           >  8:163369/1‑149 (MQ=255)
                                                                  cTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAgttgct                                                                           <  5:224235/149‑3 (MQ=255)
                                                                       cGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCTTACCGAGGGGCCAGTTGTTCCTCa                                                                      <  7:40391/149‑1 (MQ=255)
                                                                             cGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGAcc                                                                <  1:217702/149‑1 (MQ=255)
                                                                                ggAAGCCGAACTCGTCACCCAGCGCACCGGCTTGCCACGGGCACAGCTACAGGAAGCACGCCGGGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCGACGCCGAGGGGCCAGTTGTTCCTCAATGACctg                                                              <  8:183433/148‑1 (MQ=255)
                                                                                 gAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACctgct                                                            <  1:118544/149‑1 (MQ=255)
                                                                                 gAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACctgct                                                            <  6:22000/149‑1 (MQ=255)
                                                                                 gAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACctgc                                                             <  8:25589/148‑1 (MQ=255)
                                                                                   aGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACctgctg                                                           >  8:3707/1‑148 (MQ=255)
                                                                                    gCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACctgctg                                                           >  5:97621/1‑147 (MQ=255)
                                                                                       gAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTa                                                      >  8:217916/1‑149 (MQ=255)
                                                                                        aaCTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGAt                                                                                                           >  2:81357/1‑95 (MQ=255)
                                                                                        aaCTGTTCACCCACCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGAt                                                                                                           <  1:81357/95‑1 (MQ=255)
                                                                                            gTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCt                                                 >  8:14039/1‑149 (MQ=255)
                                                                                                  ccAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGAc                                                                                                       <  4:219108/89‑1 (MQ=255)
                                                                                                  ccAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGAc                                                                                                       >  3:219108/1‑89 (MQ=255)
                                                                                                      cgcACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAgg                                       <  4:74807/149‑1 (MQ=255)
                                                                                                          ccGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGa                                   <  7:3707/149‑1 (MQ=255)
                                                                                                             gCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCAGAACAAAAGGGCCATTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGAGCTGAGGATGAcc                                 >  5:78351/1‑148 (MQ=255)
                                                                                                                   cACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATGGAt                          >  3:78136/1‑149 (MQ=255)
                                                                                                                        gCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATTGATCTgg                      <  7:174856/148‑1 (MQ=255)
                                                                                                                              cTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAgg                                                                                      >  7:77664/1‑78 (MQ=255)
                                                                                                                              cTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAgg                                                                                      <  8:77664/78‑1 (MQ=255)
                                                                                                                               tACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCATTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCAAGTTGTTCCTCAATGACCTGCTGAAGTATTTCTAGACATGAGGATGACCCATGCATCTGGTACTTGa               >  4:242976/1‑148 (MQ=255)
                                                                                                                                 cAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGTGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATGGATCTGGTACTTGATc             <  3:242976/148‑1 (MQ=255)
                                                                                                                                    gAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATGGATCTGGTa                    >  4:146548/1‑138 (MQ=255)
                                                                                                                                    gAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATGGATCTGGTa                    <  3:146548/138‑1 (MQ=255)
                                                                                                                                    gAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCGTGACCTGAGGATGACCCATGGATCTGGTACGTGATCTGCt         >  1:210561/1‑149 (MQ=255)
                                                                                                                                           gCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATGGATCTGGTACTTGATCTGCTCGCGACg  <  7:217916/149‑1 (MQ=255)
                                                                                                                                                   gccgAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGcc                                                                                  <  2:182958/61‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
AGGGCCGGCGAGAAGCTGCTGCCAGTCGATGAACTGCCGTTCGAGTTCCTGATGAACGCCCTGCGCCTGACCGATGGCGTGGAAGCCGAACTGTTCACCCAGCGCACCGGCTTGCCACTGGCACAGCTACAGGAAGCACGCCGTGCCGCCGAACAAAAGGGCCTTTTACAGGTCGAACCGGATCGACTGGTAGCCACGCCGAGGGGCCAGTTGTTCCTCAATGACCTGCTGCAGTATTTCTTGACCTGAGGATGACCCATGGATCTGGTACTTGATCTGCTCGCGACG  >  NC_002947/5823992‑5824279

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