Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,282,627 C→T Q922* (CAG→TAG)  narG → nitrate reductase A subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,282,6270CT89.3% 69.0 / ‑4.2 30Q922* (CAG→TAG) narGnitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base T (6/19);  minor base A (1/2);  total (9/21)
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

CATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTC  >  NC_000913/1282404‑1282864
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
cATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGAGTGCCCGGGACACCCagaaggggaagggtctgctggggcaac                                                                                                                                                                                                                       >  1:502025/1‑227 (MQ=255)
                                     gAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCa                                                                                                                                                                                  <  1:113231/251‑1 (MQ=255)
                                            cgcAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGGATGCCTGGGACACCAAGAACGAGAAGGATCTGCTGCGTAAGCTCAACTACCACAAAACGGAAGGTCCGGCGGAAGGGCAGCCggtg                                                                                                                                                                           >  1:310758/1‑248 (MQ=255)
                                                  cGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATTCCAGGTAAAACCGCGACACATTTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAac                                                                                                                                                                     <  2:419866/251‑1 (MQ=255)
                                                  cGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAac                                                                                                                                                                     <  1:334330/251‑1 (MQ=255)
                                                      ggATGAGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAACCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACTCCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAacac                                                                                                                                                                   <  2:328187/249‑1 (MQ=255)
                                                                       aaaaaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGTCACACATTATGGTCGTAGAGCGCGAGTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCTAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCagag                                                                                                                                               <  1:99649/251‑1 (MQ=255)
                                                                         aaaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAgat                                                                                                                                             <  1:178114/251‑1 (MQ=255)
                                                                            ggCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAgatgat                                                                                                                                          <  1:295603/251‑1 (MQ=255)
                                                                            ggCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAgatgat                                                                                                                                          <  2:131922/251‑1 (MQ=255)
                                                                               gAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTGGCCGGAACACCTAGAGCGAGATGGGGCTGCTGGGTAAGCTCAACTACACCAAAACGGAAGGGCCGGCGAAAGGCCCGCCCATGCTGAACACCCCCATTGATGCGGCCGAgaagatacc                                                                                                                                       >  1:320348/1‑248 (MQ=255)
                                                                                     gACCTGATCCCAGGTAAAACCTCGCCAAACATTATGTGCGCAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACt                                                                                                                                 <  2:564398/251‑1 (MQ=255)
                                                                                           aTCCAAGGTAAAACCGCGCCACACATTATGTTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCAcc                                                                                                                           <  2:149627/251‑1 (MQ=255)
                                                                                                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGt                                                                                                  <  1:521242/251‑1 (MQ=255)
                                                                                                                            cacgccgcgctgccgatctcGACGTACGAACGTTTTACCTCTATCGGCCCGCTGATGGAGAAAAACGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTg                                                                                         <  2:265737/232‑1 (MQ=255)
                                                                                                                                         aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAAGGATTGCCCGGGACACCTAGAGCGAGATGGATCTGCGGCGGAAGCTCAAATACACCACAGAGGAAAGTCCGGCGGAAGGGCAGCCGGGGCTGAACACCGCAATTGAAGTGGGCAGAAGGATCCTGACCCCGGACCCGGACACCAACGGTCAAGGAGGCCGGAAAGCCTGGGGTGCCCCGAGc                                                                            >  1:301007/1‑251 (MQ=255)
                                                                                                                                             ggCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAAt                                                                       <  1:350737/251‑1 (MQ=255)
                                                                                                                                                           gCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGGACACCTAGAGCGAGATGGGTCTGCTGCGTAAGCTCAACTACACCAAAGAGGAAGGGCCGGCGCAAGGCCAGCCGATGCTGAACACCGCAATTTATGCGGCAGAGATGATCCTGACACTGGCACCGGGAAACAACGGGCAGGGAGCCGTGAAAGCCTGGGCCGCCCTTAGCGAATTTACCGGTCggggc                                                         >  1:298315/1‑247 (MQ=255)
                                                                                                                                                               tACCGCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGACCGGCGAAAGGCCAGCCTATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTg                                                                                                                                      <  2:510670/172‑1 (MQ=255)
                                                                                                                                                                     tATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGGTTTCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAACGGAAGGGCCGGCGAAAGGGCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGGAAGCCTGGGCTGCCCTGAGCGGATTTACCGGTCGGGACCATACGCATc                                               >  1:347423/1‑251 (MQ=255)
                                                                                                                                                                      aTCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCAt                                                <  1:293756/249‑1 (MQ=255)
                                                                                                                                                                         ggCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGgc                                           <  1:8358/251‑1 (MQ=255)
                                                                                                                                                                           cccGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTTCCTGGGACACCCAGAGAGAGGTGGATCTGCTGCGTAAGGTCAACCAAAACAAAGCGGAAGGTCCGGCGGAAGGCCAGCCGGTGCTGAACACCGCAATTTATGGGGCAGAGATGAACCTTACACTGGCCCCGGAAAACAACGGGCAGGGAGCCGGGAAAGCCCGGGCGTCCCTGAGCGAGATTTAcgggcgggacccaacgcctctggggc                                          >  1:524823/1‑226 (MQ=255)
                                                                                                                                                                                                  aTGGCGGTAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATc                  <  1:213111/251‑1 (MQ=255)
                                                                                                                                                                                                     gCGGTAAAGGTATGGCCTGGAACACCTGGAGCGAGAGGGATCTGCTGCGTAAGCACAACTACCCCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGc               <  2:447792/251‑1 (MQ=255)
                                                                                                                                                                                                         tAAAGGGATTGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTAGTGACCATACGCATCTGGCTCGGAATAAAGAAGAAGAGAAGATCCGCTTcc           >  2:443023/1‑251 (MQ=255)
                                                                                                                                                                                                               ggATTGCCTGGAACACCAAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat      <  1:487373/251‑1 (MQ=255)
                                                                                                                                                                                                                   tGCCTGGACCACCAAGAGCGTGATGGATCTGCTGCGTAAGCTCAACCACCCCAACGCGGAAGGTCCGGCGAAATGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGGCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGatat  <  2:298315/250‑2 (MQ=255)
                                                                                                                                                                                                                   tGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt   <  1:528225/250‑1 (MQ=255)
                                                                                                                                                                                                                   tGCCTGGAACACCTAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTc  >  1:233554/1‑251 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
CATCGTCACGCTGCCTATCCAGCATGACTCTGCCGCTGAACTGGCGCAGCCGCTGGATGTGAAAGACTGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGAAAGGCCAGCCGATGCTGAACACCGCAATTGATGCGGCAGAGATGATCCTGACACTGGCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTC  >  NC_000913/1282404‑1282864

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