Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,446,491 G→C 10.3% Y358* (TAC→TAG)  nicC ← 6‑hydroxynicotinate 3‑monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,446,4910GC10.3% 113.1 / 5.0 39Y358* (TAC→TAG) nicC6‑hydroxynicotinate 3‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base C (2/2);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTC  >  NC_002947/4446350‑4446631
                                                                                                                                             |                                                                                                                                            
cgccATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTTCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGc                                                                                                                                        >  3:51548/1‑148 (MQ=255)
   cATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAACCCCAGGCTGGATCTTCCTGGCTGTAGAGCCAgg                                                                                                                                    >  7:116899/1‑149 (MQ=255)
    aTAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTTCCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTCTACAGCCAgga                                                                                                                                   <  7:158015/149‑2 (MQ=255)
     tAGCGCAGGTAATGCTTGCGGATGCCATTGGCGCTGACGTTACCCCCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTCTACAGCCAggaa                                                                                                                                  <  4:72928/149‑3 (MQ=255)
        cgcAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTg                                                                                                                               <  7:170261/149‑1 (MQ=255)
              tAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTg                                                                                                                         <  7:45670/149‑1 (MQ=255)
                aTGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGa                                                                                                                       <  5:166356/149‑1 (MQ=255)
                   cTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATa                                                                                                                     <  6:269560/148‑1 (MQ=255)
                      gCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGa                                                                                                                 <  7:13728/149‑1 (MQ=255)
                         gATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACGGATTg                                                                                                              >  7:86922/1‑149 (MQ=255)
                              cATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCt                                                                                                         >  3:235485/1‑149 (MQ=255)
                                      cTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATAc                                                                                                                    <  8:129156/130‑1 (MQ=255)
                                      cTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATAc                                                                                                                    >  7:129156/1‑130 (MQ=255)
                                       tGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGccc                                                                                                >  6:120198/1‑149 (MQ=255)
                                         aCGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCg                                                                                               <  4:235485/148‑1 (MQ=255)
                                                   caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGcc                                                                                                 <  3:111963/136‑1 (MQ=255)
                                                   caGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGcc                                                                                                 >  4:111963/1‑136 (MQ=255)
                                                                    aTGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAAcgc                                                                    <  1:110957/148‑1 (MQ=255)
                                                                          ccTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTCGCGGTTGGCTTCATACAACGCGAACGcc                                                             <  4:51548/149‑1 (MQ=255)
                                                                             ccccGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATg                                                                                                   <  5:122182/108‑1 (MQ=255)
                                                                             ccccGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATg                                                                                                   >  6:122182/1‑108 (MQ=255)
                                                                                   tttCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCg                                                        <  8:114149/145‑1 (MQ=255)
                                                                                   tttCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCg                                                        >  7:114149/1‑145 (MQ=255)
                                                                                      cAAGCTGCTGCCCAAAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCg                                                 >  5:224559/1‑149 (MQ=255)
                                                                                             ctgcCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                          >  3:4638/1‑149 (MQ=255)
                                                                                             ctgcCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                          >  7:139050/1‑149 (MQ=255)
                                                                                             ctgcCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTCTTGGCGTAGGAGACCGCTTGCACCTGCGATGCACGCTCATTGCGGTAGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGcc                                          >  1:217730/1‑149 (MQ=255)
                                                                                                cccATAGAGGTCGTAGCCGTAAACGCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                        <  8:116899/148‑1 (MQ=255)
                                                                                                  cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                        >  1:115778/1‑146 (MQ=255)
                                                                                                  cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCg                                        <  2:115778/146‑1 (MQ=255)
                                                                                                  cATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGtc                                     >  2:58595/1‑149 (MQ=255)
                                                                                                                        ccAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGc                                           <  5:109554/121‑1 (MQ=255)
                                                                                                                        ccAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGc                                           >  6:109554/1‑121 (MQ=255)
                                                                                                                         cAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCg                       >  5:16508/1‑140 (MQ=255)
                                                                                                                         cAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCg                       <  6:16508/140‑1 (MQ=255)
                                                                                                                              tGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCcg         >  3:117650/1‑149 (MQ=255)
                                                                                                                                ccaacaccTGGCTCTACAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCcgcg       >  3:186532/7‑149 (MQ=255)
                                                                                                                                     accTGGCTCTAAAGCAAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTc  >  7:254101/2‑149 (MQ=255)
                                                                                                                                      ccTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCt   <  4:4638/147‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CGCCATAGCGCAGGTAATGCTGGCGGATGCCATTGGCGCTGACGTTACCACCAGCGACGAAGGTGCTCATGCCGCCTCCCCGCTTTCAAGCTGCTGCCCATAGAGGTCGTAGCCGTAAACCCAGGCTGGATCTTCCTGGCTGTAGAGCCAGGTGTTGGCGTTGGATACCGATTGCACCTGCGATGCCCGCTCTTTGCGGTTGGCTTCATACAACGCGAACGCCGTGCGGTGGTCGGACAGCCCGGTCTCCTGCAGGCAGCGGGTCAGCATGGCCGCGTCCTC  >  NC_002947/4446350‑4446631

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