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% 76.0 / NA 21L98L (CTG→TTG) narWnitrate reductase molybdenum cofactor assembly chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/6);  total (15/6)

ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAA  >  NZ_CP009273/1531191‑1531359
                                                                                    |                                                                                    
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                 >  2:446097/1‑90 (MQ=255)
     tcttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCg                                                                            >  2:67779/1‑90 (MQ=255)
     tcttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCg                                                                            <  1:422608/90‑1 (MQ=255)
     tcttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCg                                                                            <  1:46301/90‑1 (MQ=255)
        ttCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAcc                                                                         <  2:79431/90‑1 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      >  2:457005/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      <  2:229057/90‑1 (MQ=255)
                             aCGCTTAAATACTCCAGATACAGCGGTAGAGAAACGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAgg                                                    >  2:8338/1‑90 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  1:286854/1‑90 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  2:345858/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                       >  2:461774/1‑90 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                       >  2:438137/1‑90 (MQ=255)
                                                 cAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCcgcg                                >  1:440289/1‑90 (MQ=255)
                                                      gTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgc                           <  2:225755/90‑1 (MQ=255)
                                                      gTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgc                           >  2:332394/1‑90 (MQ=255)
                                                      gTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATcgc                           >  2:179996/1‑90 (MQ=255)
                                                             aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                    >  1:297684/1‑90 (MQ=255)
                                                             aTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTc                    <  1:229492/90‑1 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                 >  1:93588/1‑90 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa  >  1:115594/1‑90 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa  >  1:349908/1‑90 (MQ=255)
                                                                                    |                                                                                    
ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAA  >  NZ_CP009273/1531191‑1531359

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