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,2750GA89.2% 112.1 / 2.7 37L98L (CTG→TTG) narWnitrate reductase molybdenum cofactor assembly chaperone
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (18/15);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAA  >  NZ_CP009273/1531191‑1531359
                                                                                    |                                                                                    
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                 >  1:375325/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                 >  2:550512/1‑90 (MQ=255)
aTCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCa                                                                                 >  2:82726/1‑90 (MQ=255)
 tCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAg                                                                                >  2:628680/1‑90 (MQ=255)
     tcttTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCg                                                                            <  1:496664/90‑1 (MQ=255)
       tttCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                          <  2:569986/90‑1 (MQ=255)
       tttCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGAc                                                                          <  1:667067/90‑1 (MQ=255)
         tCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCt                                                                        <  1:545686/90‑1 (MQ=255)
         tCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCt                                                                        <  1:628680/90‑1 (MQ=255)
                           gTACGCTTAAATTCTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGTCAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      <  2:519426/90‑1 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      >  1:406247/1‑90 (MQ=255)
                           gTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      >  2:808856/1‑90 (MQ=255)
                             aCGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAgg                                                    >  2:83686/1‑90 (MQ=255)
                             aCGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAgg                                                    >  1:613158/1‑90 (MQ=255)
                              cGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCcagcggtc                                                   <  1:77551/90‑5 (MQ=255)
                              cGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGt                                                   >  2:402111/1‑90 (MQ=255)
                                 ttAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGGcca                                                >  2:602574/1‑90 (MQ=255)
                                   aaaTACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTccacc                                              <  1:82726/90‑1 (MQ=255)
                                     aTACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCAt                                            <  2:150309/90‑1 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  2:423530/1‑90 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  1:491316/1‑90 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  1:171136/1‑90 (MQ=255)
                                        cTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGc                                         >  2:747458/1‑90 (MQ=255)
                                         tCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGcc                                        <  1:760893/90‑1 (MQ=255)
                                          ccAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCt                                       >  2:813749/1‑90 (MQ=255)
                                                 cAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCcgcg                                >  1:409546/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                               >  2:714409/1‑90 (MQ=255)
                                                  aGCGGTAGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGa                               >  2:85110/1‑90 (MQ=255)
                                                        aGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATa                                                                 >  2:252118/1‑50 (MQ=255)
                                                        aGATAATCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATa                                                                 <  1:252118/50‑1 (MQ=255)
                                                            aaTCGGGCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACt                     <  1:402111/90‑1 (MQ=255)
                                                                gggCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGc                 >  1:785545/1‑90 (MQ=255)
                                                                 ggCAGTTCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCa                <  1:101999/90‑1 (MQ=255)
                                                                      ttCCCGACAATCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAAc           <  1:490045/90‑1 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      >  2:158727/1‑38 (MQ=39)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCcagca                                                      <  1:158727/38‑1 (MQ=255)
                                                                               aTCCAACTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGaa  >  1:722918/1‑90 (MQ=255)
                                                                                    |                                                                                    
ATCCTTCTTTCGCCTGATCGTCAGGCAGTACGCTTAAATACTCCAGATACAGCGGTAGATAATCGGGCAGTTCCCGACAATCCAGCTGCAGGCCGACCTTTTCATACTCCGCCAGCAGGTCCACCATTGCCTGGCCGCGATCGCGGGACTCGGCATGAACATGTTCGAA  >  NZ_CP009273/1531191‑1531359

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