Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,902,220 A→G 15.9% D924G (GAC→GGC)  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,902,2200AG15.9% 99.3 / 11.6 44D924G (GAC→GGC) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base G (3/4);  total (21/23)
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.71e-01

CCGACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACC  >  NC_002947/1902088‑1902364
                                                                                                                                    |                                                                                                                                                 
ccGACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGGCCAGGCGCTGGTCGAg                                                                                                                                   >  5:32938/1‑149 (MQ=255)
 cGACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGGCCAGGCGCTGGTCGAGc                                                                                                                                  <  5:320917/149‑1 (MQ=255)
   aCAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGGCCAGGCGCTGGTCGAGCtt                                                                                                                                <  6:32938/149‑2 (MQ=255)
        ggcgcggcCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTTCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCTCAAGCTCGACCAGCGGCTGGCCGATCTGGGAgg                                                                                                                           <  8:284324/149‑1 (MQ=255)
                 ttCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGGCCAGGCGCTGGTCGAGCttgcg                                                                                                                             <  4:7683/138‑5 (MQ=255)
                 ttCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGGCCAGGCGCTGGTCGAGCttgcg                                                                                                                             >  3:7683/1‑134 (MQ=255)
                      gCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCctgctg                                                                                                              <  4:303734/148‑1 (MQ=255)
                      gCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTgg                                                                                                             >  1:168819/1‑149 (MQ=255)
                      gCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTgg                                                                                                             >  3:138499/1‑149 (MQ=255)
                           cTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCgg                                                                                                        <  2:168819/149‑1 (MQ=255)
                                cAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGa                                                                                                   <  6:194915/149‑1 (MQ=255)
                                    gCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGc                                                                                               <  1:96026/149‑1 (MQ=255)
                                          aggGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGccc                                                                                         <  2:225456/149‑1 (MQ=255)
                                          aggGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGccc                                                                                         <  4:138499/149‑1 (MQ=255)
                                                cTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTc                                                                                   <  8:55935/149‑1 (MQ=255)
                                                cTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTc                                                                                   >  8:102305/1‑149 (MQ=255)
                                                              cTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTgc                                                                     >  6:320373/1‑149 (MQ=255)
                                                                  ccTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTgctggc                                                                 >  3:2701/1‑149 (MQ=255)
                                                                   cTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTgctggct                                                                >  4:112591/1‑149 (MQ=255)
                                                                   cTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTgctggct                                                                <  3:112591/149‑1 (MQ=255)
                                                                    tATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCAAGCCCGACTTCGACGAAGCCTTTgctggctg                                                               >  5:92775/1‑149 (MQ=255)
                                                                      tGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTgctggctg                                                               <  1:8789/147‑1 (MQ=255)
                                                                      tGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCtggctggc                                                             >  7:223065/1‑149 (MQ=255)
                                                                        ccGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTc                                                           >  4:18584/1‑149 (MQ=255)
                                                                         cGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTc                                                           >  7:197317/1‑148 (MQ=255)
                                                                           tACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑aa                                                       <  5:320373/149‑1 (MQ=255)
                                                                              tGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCg                                                    >  7:177660/1‑149 (MQ=255)
                                                                                   cTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTAtt                                               <  4:2701/149‑1 (MQ=255)
                                                                                      ggCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGcc                                            >  2:11630/1‑149 (MQ=255)
                                                                                               aGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTggg                                   <  3:18584/149‑1 (MQ=255)
                                                                                                 ttACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGc                                  >  1:45167/1‑148 (MQ=255)
                                                                                                     gACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGc                             <  7:102305/149‑1 (MQ=255)
                                                                                                       ccAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGCGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCcg                           >  4:221009/1‑149 (MQ=255)
                                                                                                            tCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCgcggc                       <  8:177660/148‑1 (MQ=255)
                                                                                                                cGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGgcgc                     <  8:197317/146‑1 (MQ=255)
                                                                                                                 ggTTTTGGCCGCAAGCTCGACCAGCGCCTGGCGGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCAtc                 <  8:223065/149‑1 (MQ=255)
                                                                                                                 ggTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCAAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCAtc                 <  6:92775/149‑1 (MQ=255)
                                                                                                                    tttGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTcc              >  2:144211/1‑149 (MQ=255)
                                                                                                                         ctcccaGCTCGGCCAGGCGCTGGTCGAGCTTTGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCGCAAGCGTTATTGCCTGCACTGGGCACTGCCgcggc                       <  1:205708/131‑1 (MQ=255)
                                                                                                                         ctcccaGCTCGGCCAGGCGCTGGTCGAGCTTGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCgcggc                       >  2:205708/6‑135 (MQ=255)
                                                                                                                         ccGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGc           >  3:105606/1‑147 (MQ=255)
                                                                                                                                cTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAAcc  >  2:27341/1‑149 (MQ=255)
                                                                                                                                 tCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCAtctc               <  6:204121/135‑1 (MQ=255)
                                                                                                                                 tCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCAtctc               >  5:204121/1‑135 (MQ=255)
                                                                                                                                    |                                                                                                                                                 
CCGACAGCGGCGCGGCCTTCTGGCAGGCTTTGCAGGGCGAGGAGGGCGCTCACTGCGCCGGGCTGCCCTATGCCGTACTGGCCCTGGGCGACTCCAGTTACGACCAGTTCTGCGGTTTTGGCCGCAAGCTCGACCAGCGCCTGGCCGAGCTGGGAGGGCATCGCCTGCTGGCGCGGGTGGACTGCGAGCCCGACTTCGACGAAGCCTTTGCTGGCTGGCTCG‑AAGCGTTATTGCCTGCACTGGGCACTGCCGCGGCGCCATCTCCGGCGAGTGAACC  >  NC_002947/1902088‑1902364

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