Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,531,275 G→A L98L (CTG→TTG)  narW ← nitrate reductase molybdenum cofactor assembly chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,531,2750GA76.3% 74.4 / 19.4 38L98L (CTG→TTG) narWnitrate reductase molybdenum cofactor assembly chaperone
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (18/11);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAAC  >  NZ_CP009273/1531191‑1531360
                                                                                    |                                                                                     
aTCCTTCTTTCGCCTGATGGTCAAGCAGTACACTTAAATACACCAGATACAGCGGTAGGTAATCGGGCAGTTTGCGAGAATTCAAGTGCa                                                                                  >  2:221849/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  1:233071/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                  >  1:40674/1‑90 (MQ=255)
   cttcttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGc                                                                               <  2:1033231/90‑1 (MQ=255)
      cttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGa                                                                            <  1:127980/90‑1 (MQ=255)
       tttCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           <  2:574354/90‑1 (MQ=255)
          cGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCtt                                                                        >  1:767088/1‑90 (MQ=255)
           gCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCttt                                                                       >  2:865789/1‑90 (MQ=255)
            ccTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCtttt                                                                      >  2:136739/1‑90 (MQ=255)
               gATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCAt                                                                   >  1:627023/1‑90 (MQ=255)
                          aGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCcagc                                                        <  2:429327/90‑1 (MQ=255)
                          aGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagc                                                        <  1:865789/90‑1 (MQ=255)
                          aGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagc                                                        >  2:1050049/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  1:933217/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       <  1:1050049/90‑1 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                          >  1:395006/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  2:91025/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  1:454715/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                        >  2:173852/1‑90 (MQ=255)
                                            aGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTgg                                      >  1:966375/1‑90 (MQ=255)
                                                      gTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           <  2:1092983/43‑1 (MQ=255)
                                                      gTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                           >  1:1092983/1‑43 (MQ=255)
                                                           tAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGAc                       >  2:314347/1‑90 (MQ=255)
                                                            aaTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACt                      <  1:711409/90‑1 (MQ=255)
                                                             aTCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     >  2:746191/1‑90 (MQ=255)
                                                             aTCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     <  1:55579/90‑1 (MQ=255)
                                                             aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                     >  1:296150/1‑90 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       >  1:514879/1‑53 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCcagca                                                       <  2:514879/53‑1 (MQ=255)
                                                                      ttCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAAc            <  2:966375/90‑1 (MQ=255)
                                                                      ttCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAAc            <  2:556057/90‑1 (MQ=255)
                                                                      tccccGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAAc            <  2:40674/88‑1 (MQ=255)
                                                                           gACAATCCAACTGCAGGCCGCCCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGtt       <  2:1017412/90‑1 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa   >  2:650988/1‑90 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa   >  2:818641/1‑90 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa   >  1:114022/1‑90 (MQ=255)
                                                                                tCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAAc  <  1:1086952/90‑1 (MQ=255)
                                                                                tCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAAc  <  2:767088/90‑1 (MQ=255)
                                                                                    |                                                                                     
ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAAC  >  NZ_CP009273/1531191‑1531360

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