Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,903,567 G→C 17.9% intergenic (+50/‑221) PP_1703 → / → nirB assimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,903,5670GC17.9% 95.3 / 12.7 41intergenic (+50/‑221)PP_1703/nirBassimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base G (16/16);  new base C (3/4);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTA  >  NC_002947/1903420‑1903715
                                                                                                                                                   |                                                                                                                                                    
ggATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTgggcggg                                                                                                                                                    <  2:954619/149‑1 (MQ=255)
           gAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaatg                                                                                                                                         <  1:1259337/150‑6 (MQ=255)
               aTGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgc                                                                                                                                              >  1:875883/1‑141 (MQ=255)
               aTGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgc                                                                                                                                              <  2:875883/141‑1 (MQ=255)
                  gtcGGTTGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGa                                                                                                                                  >  2:330313/3‑150 (MQ=255)
                       atggatgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCggg                                                                                                                             <  1:213103/150‑1 (MQ=255)
                            tgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgcc                                                                                                                                             <  2:635808/129‑2 (MQ=255)
                            tgCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgcc                                                                                                                                             >  1:635808/1‑128 (MQ=255)
                             gcggcggcgGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTAcc                                                                                                                       <  2:395511/150‑1 (MQ=255)
                                     ggctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaa                                                                                                                                           >  1:1246186/1‑119 (MQ=255)
                                     ggctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCCCGCCCAgccaa                                                                                                                                           <  2:1246186/122‑4 (MQ=255)
                                      gctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg                                                                                                                              <  1:763584/134‑1 (MQ=255)
                                      gctgcggctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg                                                                                                                              >  2:763584/1‑134 (MQ=255)
                                            gctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCggg                                                                                                                             <  1:1087602/129‑1 (MQ=255)
                                            gctTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCggg                                                                                                                             >  2:1087602/1‑129 (MQ=255)
                                                     gAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaac                                                                                               <  1:1282901/150‑1 (MQ=255)
                                                           ttGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggc                                                                                         >  1:370860/1‑150 (MQ=255)
                                                               gCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggcaacc                                                                                     >  2:1096183/1‑150 (MQ=255)
                                                                cAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg                                                                                                                              >  2:344575/1‑108 (MQ=255)
                                                                cAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCgg                                                                                                                              <  1:344575/108‑1 (MQ=255)
                                                                 aaGGCCAGGCGTTATCGGCGCGATGTGTACTTAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAggcaacccaggcaaccca                                                                                   <  2:265979/150‑1 (MQ=255)
                                                                     ccAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTTTCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCa                                                                               >  1:1168669/1‑150 (MQ=255)
                                                                      cAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCaa                                                                              >  1:408252/1‑150 (MQ=255)
                                                                           gTTAGCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGa                                                                         >  2:548362/1‑150 (MQ=255)
                                                                             tATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGCCTGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGAtt                                                                       <  1:737290/150‑1 (MQ=255)
                                                                                  gcgcGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggcaa                                                                                       >  2:97850/1‑129 (MQ=255)
                                                                                  gcgcGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggcaacccaggcaa                                                                                       <  1:97850/129‑1 (MQ=255)
                                                                                         gtgtACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACTCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCt                                                           >  2:404814/1‑150 (MQ=255)
                                                                                           gtACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACTCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGc                                                         <  1:404814/150‑1 (MQ=255)
                                                                                            tACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCAc                                                                  <  2:1335210/140‑1 (MQ=255)
                                                                                            tACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCAc                                                                  >  1:1335210/1‑140 (MQ=255)
                                                                                                     gCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCgg                                               <  1:1283769/150‑1 (MQ=255)
                                                                                                       tCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGGGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag                                                                                                    >  1:904946/1‑95 (MQ=255)
                                                                                                       tCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGGGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag                                                                                                    <  2:904946/95‑1 (MQ=255)
                                                                                                        ccAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCcgc                                                                                                                    <  1:71645/78‑1 (MQ=255)
                                                                                                        ccAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCcgc                                                                                                                    >  2:71645/1‑78 (MQ=255)
                                                                                                               ggggTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGTAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACaaa                                     <  1:330313/150‑1 (MQ=255)
                                                                                                                 ggTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag                                                                                                    <  1:269664/85‑1 (MQ=255)
                                                                                                                 ggTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcag                                                                                                    >  2:269664/1‑85 (MQ=255)
                                                                                                                    gttgtGCACCTTGAATGGGCCTTGGCTGGCCATGCCCAACGCAAGTGTAGGAGCGGGTTTACGCAATTGCGAGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCa                                <  1:1096183/150‑1 (MQ=255)
                                                                                                                         gCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTCCCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGAt                           <  1:57931/150‑1 (MQ=255)
                                                                                                                                              tGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAGCCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGccc      >  2:1109207/1‑150 (MQ=255)
                                                                                                                                                 gcgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggc                                                                                                  <  1:544581/55‑1 (MQ=255)
                                                                                                                                                 gcgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACAcaggc                                                                                                  >  2:544581/1‑55 (MQ=255)
                                                                                                                                                  ggggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAGCCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTa  <  1:1109207/149‑1 (MQ=255)
                                                                                                                                                  cgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCt                                                           >  1:282232/1‑93 (MQ=255)
                                                                                                                                                  cgggcCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCt                                                           <  2:282232/93‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGATTGTCGCCGAGCATGGCCGGATGGATGCGGCGGCGGCTGCGGCTTATGTCGAAGGCTTGAGCAAGGCCAGGCGTTATCGGCGCGATGTGTACTGAAGTGCTCCAAGATGGGGTGTTGTGCACCTTGAATGTGCATTGGCTGGGCGGGCCCAGCGCAAGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTA  >  NC_002947/1903420‑1903715

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: