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,2750GA100.0% 154.0 / NA 40L98L (CTG→TTG) narWnitrate reductase molybdenum cofactor assembly chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/13);  total (27/13)

CAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAA  >  NZ_CP009273/1531189‑1531359
                                                                                      |                                                                                    
cAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTg                                                                                   <  1:779136/90‑1 (MQ=255)
  aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                 >  2:1147399/1‑90 (MQ=255)
  aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                 >  1:295582/1‑90 (MQ=255)
         tttCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                          <  1:6183/90‑1 (MQ=255)
            cGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCtt                                                                       >  2:133366/1‑90 (MQ=255)
                 gATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCAt                                                                  >  2:727328/1‑90 (MQ=255)
                     gTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTc                                                              <  1:231328/90‑1 (MQ=255)
                            aGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagc                                                       >  1:81772/1‑90 (MQ=255)
                             gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      >  2:701563/1‑90 (MQ=255)
                                 gCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTc                                                  <  2:296170/90‑1 (MQ=255)
                                 gCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTc                                                  <  1:936348/90‑1 (MQ=255)
                                   ttAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTcca                                                >  1:29728/1‑90 (MQ=255)
                                         aCTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTg                                          >  2:574381/1‑90 (MQ=255)
                                          cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  2:517311/1‑90 (MQ=255)
                                          cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  2:1000052/1‑90 (MQ=255)
                                            ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                       >  2:810935/1‑90 (MQ=255)
                                            ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                       >  2:1189246/1‑90 (MQ=255)
                                            ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                       >  2:1112116/1‑90 (MQ=255)
                                                  aCAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCcgc                                 >  2:233084/1‑90 (MQ=255)
                                                   cAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCcgcg                                >  1:1205864/1‑90 (MQ=255)
                                                    aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAcc                                                                         >  1:893332/1‑48 (MQ=255)
                                                    aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAcc                                                                         <  2:893332/48‑1 (MQ=255)
                                                    aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                               >  1:701611/1‑90 (MQ=255)
                                                    aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                               >  1:402400/1‑90 (MQ=255)
                                                             tAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGAc                      >  1:1035098/1‑90 (MQ=255)
                                                               aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                    <  1:1139155/90‑1 (MQ=255)
                                                               aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                    >  2:315831/1‑90 (MQ=255)
                                                               aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                    <  2:1139174/90‑1 (MQ=255)
                                                               aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                    >  1:467953/1‑90 (MQ=255)
                                                                  gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagc                                                       <  1:966106/52‑1 (MQ=255)
                                                                  gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagc                                                       >  2:966106/1‑52 (MQ=255)
                                                                  gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                 >  2:134430/1‑90 (MQ=255)
                                                                  gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                 >  1:920626/1‑90 (MQ=255)
                                                                  gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                 >  1:222644/1‑90 (MQ=255)
                                                                   ggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCa                <  2:975468/90‑1 (MQ=255)
                                                                        ttCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAAc           <  2:198789/90‑1 (MQ=255)
                                                                        ttCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAAc           <  2:16914/90‑1 (MQ=255)
                                                                         tCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACa          <  2:295582/90‑1 (MQ=255)
                                                                                 aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa  >  2:707206/1‑90 (MQ=255)
                                                                                 aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa  >  2:95101/1‑90 (MQ=255)
                                                                                      |                                                                                    
CAATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAA  >  NZ_CP009273/1531189‑1531359

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