Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,849 A→T 10.9% K120* (AAG→TAG) ‡ PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,8490AT10.9% 127.7 / 7.7 46K120* (AAG→TAG) ‡PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base T (2/3);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

CAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTG  >  NC_002947/1899705‑1899998
                                                                                                                                                |                                                                                                                                                      
cagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCAACCGCCGAGCGTCTGCGCGCCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGGGTCCCAGTCCCTGTCCCACAAGTTg                                                                                                                                                   >  2:255676/1‑150 (MQ=255)
   ccgccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGcc                                                                                                                                                >  2:648477/1‑150 (MQ=255)
    cgccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACTTGTTGGCCa                                                                                                                                               <  2:466322/150‑1 (MQ=255)
         gTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGc                                                                                                                                          <  1:648477/150‑1 (MQ=255)
           gcgcACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACTTGTTGGCCaggtactg                                                                                                                                        <  1:116610/150‑8 (MQ=255)
            cgcACCAGCATGGACCAGGCCATAGCCGCGACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCttt                                                                                                                                       >  1:150946/1‑150 (MQ=255)
                   gCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTa                                                                                                                                <  1:94821/150‑1 (MQ=255)
                    cATGGACCAGGCCATCGCTGCCACAGCCGAGCGTCATCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTAc                                                                                                                               >  1:88876/1‑150 (MQ=255)
                     aTGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTAcg                                                                                                                              >  2:1090408/1‑150 (MQ=255)
                         aCCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATAATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAg                                                                                                                                                      <  2:283018/122‑1 (MQ=255)
                         aCCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATAATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAg                                                                                                                                                      >  1:283018/1‑122 (MQ=255)
                           cAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCcac                                                                                                                        <  1:788044/150‑1 (MQ=255)
                              gCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATc                                                                                                                     <  1:223198/150‑1 (MQ=255)
                                 aTCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAg                                                                                                                  >  2:855402/1‑150 (MQ=255)
                                 aTCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAg                                                                                                                  >  1:733064/1‑150 (MQ=255)
                                  tCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGt                                                                                                                 <  1:823106/150‑1 (MQ=255)
                                             gccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCgcgc                                                                                                      <  1:450859/150‑1 (MQ=255)
                                             gccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCgcgc                                                                                                      >  1:766502/1‑150 (MQ=255)
                                                 aGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCtgt                                                                                                  <  2:94918/150‑1 (MQ=255)
                                                   cGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGTCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGc                                                                                                <  2:968125/150‑1 (MQ=255)
                                                       tGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATgg                                                                                            >  2:473508/1‑150 (MQ=255)
                                                            ggCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgcca‑gc                                                                                      <  1:855402/150‑1 (MQ=255)
                                                               atcatcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgcca‑gcgcc                                                                                   >  2:929339/1‑150 (MQ=255)
                                                                  atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgcca‑gcgccagc                                                                                >  1:186661/1‑150 (MQ=255)
                                                                     gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                                >  2:329065/1‑116 (MQ=255)
                                                                     gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                                <  1:329065/116‑1 (MQ=255)
                                                                                     aTGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGc                                                             <  1:1090408/150‑1 (MQ=255)
                                                                                                cTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGAt                                                  <  1:375343/150‑1 (MQ=255)
                                                                                                  gTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGCCCATCTCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGCTgg                                                <  1:900115/150‑1 (MQ=255)
                                                                                                     cGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATgg                                                >  2:53919/1‑147 (MQ=255)
                                                                                                     cGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATgg                                                <  1:53919/147‑1 (MQ=255)
                                                                                                        gtCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGgcgc                                                      <  1:628194/138‑1 (MQ=255)
                                                                                                        gtCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGgcgc                                                      >  2:628194/1‑138 (MQ=255)
                                                                                                         tCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCt                                         <  2:766502/150‑1 (MQ=255)
                                                                                                          ccGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTg                                        <  1:929339/150‑1 (MQ=255)
                                                                                                            gggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGACCAAATCGCGACTGTGCATGGCCA‑GCGGCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGc                                      >  2:124272/1‑150 (MQ=255)
                                                                                                              gCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAg                                    >  1:937037/1‑150 (MQ=255)
                                                                                                                  aTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTAc                                <  2:1077993/150‑1 (MQ=255)
                                                                                                                        cTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGAc                          >  2:685672/1‑150 (MQ=255)
                                                                                                                         tGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACt                         >  1:716473/1‑150 (MQ=255)
                                                                                                                             ggCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGa                     <  1:685672/150‑1 (MQ=255)
                                                                                                                                 cAGTACCTGGCCAACAAGTTGGCCAATGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCAc                 >  2:372842/1‑150 (MQ=255)
                                                                                                                                  aGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACg                <  2:937037/150‑1 (MQ=255)
                                                                                                                                     ccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATg                                                                                             >  2:110119/4‑71 (MQ=38)
                                                                                                                                     ccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATg                                                                                             <  1:110119/68‑1 (MQ=38)
                                                                                                                                      cttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACCCGATGTCAGCCGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGa            >  2:870262/3‑150 (MQ=255)
                                                                                                                                            caacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAtt      >  1:488736/1‑150 (MQ=255)
                                                                                                                                            caacaaGTTGGCCAAGGGCTTTATCCGTACGAGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAtt      >  1:113159/1‑150 (MQ=255)
                                                                                                                                                aaGTTGGACATGGGATTAAACCAAACTCTCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGGCGGCTAC‑GTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTg  <  1:1092978/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCA‑GCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTG  >  NC_002947/1899705‑1899998

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