Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,284,519 T→C S181P (TCC→CCC)  narP → nitrate/nitrite response regulator protein NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,284,5190TC98.1% 169.4 / ‑4.1 53S181P (TCC→CCC) narPnitrate/nitrite response regulator protein NarP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/21);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

GCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTT  >  NZ_CP009273/2284433‑2284601
                                                                                      |                                                                                  
gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCg                                                                                 >  1:1002831/1‑90 (MQ=255)
gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCg                                                                                 >  1:789703/1‑90 (MQ=255)
gCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCg                                                                                 >  1:20911/1‑90 (MQ=255)
  gTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAg                                                                               >  2:1066587/1‑90 (MQ=255)
   tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGc                                                                              >  2:882259/1‑90 (MQ=255)
   tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGc                                                                              <  1:1066587/90‑1 (MQ=255)
   tGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGc                                                                              <  1:352426/90‑1 (MQ=255)
     cTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                            <  2:1401584/90‑1 (MQ=255)
        aCGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACa                                                                         <  1:1044862/90‑1 (MQ=255)
        aCGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGACCAGACa                                                                         <  2:1041675/90‑1 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  2:1012292/1‑90 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  1:350510/1‑90 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    <  2:714392/90‑1 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  2:734565/1‑90 (MQ=255)
             gcgcGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTaaa                                                                    >  1:119436/1‑90 (MQ=255)
                    cTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACat                                                             <  1:5332/90‑1 (MQ=255)
                       gATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATAtt                                                          <  2:548123/90‑1 (MQ=255)
                            tCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCaa                                                     >  2:1545208/1‑90 (MQ=255)
                            tCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCaa                                                     >  1:1447280/1‑90 (MQ=255)
                              tGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATc                                                   >  2:1465242/1‑90 (MQ=255)
                              tGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATc                                                   >  1:986654/1‑90 (MQ=255)
                                cACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATctg                                                 <  1:1415319/90‑1 (MQ=255)
                                  cGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGCTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATctgct                                               <  1:1545208/90‑1 (MQ=255)
                                    aGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACat                                                             >  2:375593/1‑74 (MQ=255)
                                    aGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACat                                                             <  1:375593/74‑1 (MQ=255)
                                    aGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCtgctgc                                             >  1:692588/1‑90 (MQ=255)
                                      cTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGt                                           <  2:892853/90‑1 (MQ=255)
                                        ggCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                         <  2:351952/90‑1 (MQ=255)
                                         gCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaaa                                        >  2:1339626/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                            >  2:1030542/1‑49 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAg                                                                            <  1:1030542/49‑1 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                         <  2:1293386/84‑1 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTaa                                         >  1:1293386/1‑84 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  2:1160533/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  2:1200916/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  1:49455/1‑90 (MQ=255)
                                              ggggCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCaa                                   >  1:1197445/1‑90 (MQ=255)
                                                  cTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTc                               >  2:1539381/1‑90 (MQ=255)
                                                   tGTCAAATAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTcc                              >  2:594295/1‑90 (MQ=255)
                                                       aaaTAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTc                          <  1:707964/90‑1 (MQ=255)
                                                        aaTAAACAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCa                         >  2:1436788/1‑90 (MQ=255)
                                                              cAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTg                   >  2:741816/1‑90 (MQ=255)
                                                              cAGATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTg                   >  1:113278/1‑90 (MQ=255)
                                                                gATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                 >  1:177299/1‑90 (MQ=255)
                                                                gATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                 >  2:1469225/1‑90 (MQ=255)
                                                                gATTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggc                 >  2:679312/1‑90 (MQ=255)
                                                                 aTTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcg                <  2:78517/90‑1 (MQ=255)
                                                                 aTTGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcg                <  1:594295/90‑1 (MQ=255)
                                                                   tGCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTggcggc              <  2:177299/90‑1 (MQ=255)
                                                                    gCCTCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGcc             <  1:259197/90‑1 (MQ=255)
                                                                       tCGGTGTTGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGccacc          <  2:1581328/90‑1 (MQ=255)
                                                                             ttGAATATTCCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTg    <  1:208914/90‑1 (MQ=255)
                                                                               gAATATTCCCGAGCAGACAATAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGtt  >  1:462762/1‑90 (MQ=255)
                                                                                      |                                                                                  
GCGTGCTGACGGAGCGCGAGCTGGATGTTCTGCACGAGCTGGCACAGGGGCTGTCAAATAAACAGATTGCCTCGGTGTTGAATATTTCCGAGCAGACAGTAAAAGTACATATTCGCAATCTGCTGCGTAAACTCAATGTCCGCTCACGCGTGGCGGCCACCATTCTGTT  >  NZ_CP009273/2284433‑2284601

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