Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,884,026 T→C D411G (GAC→GGC)  cysJ ← NADPH‑dependent assimilatory sulfite reductase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,884,0260TC100.0% 66.1 / NA 20D411G (GAC→GGC) cysJNADPH‑dependent assimilatory sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/12);  total (8/12)

CCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGG  >  NZ_CP009273/2883895‑2884157
                                                                                                                                   |                                                                                                                                   
ccGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACg                                                                                                                              <  2:513949/139‑1 (MQ=255)
 cGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACg                                                                                                                              >  2:36899/1‑138 (MQ=255)
   gggTTACTGGTTTGGCTGGCAGGCGGAAGTTAGCGTTATGTCCGAGAAATACGCGGTCTTCGCCCCCGTCTCCCACGCGTTCAGCGAGGAAGCTGGAGGCACCCCCGGCGCGGGCTGGGCCTGCCACGCCGTAACGcac                                                                                                                           <  2:273803/139‑1 (MQ=255)
        tCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGTGGACTTCGCCCTCTTCTTACACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGCAACGcaccacac                                                                                                                      >  1:273803/1‑139 (MQ=255)
            gATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCaa                                                                                                                   <  2:146322/138‑1 (MQ=255)
             aTTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACg                                                                                                                 <  1:541461/139‑1 (MQ=255)
                                          gTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGt                                                                                                               >  1:460156/1‑112 (MQ=255)
                                          gTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGt                                                                                                               <  2:460156/112‑1 (MQ=255)
                                                 aaaTACGCGGACTTCGCCCTCTTCTTCCACGCTGTCGGCGAGAAAGCTGGAGGCACCACCGGCAGGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgag                                                                             <  1:160638/139‑1 (MQ=255)
                                                        cgGACTTCGCCCTCGTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                                                      <  1:442504/139‑1 (MQ=255)
                                                            cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATgg                                                                   >  2:213366/1‑138 (MQ=255)
                                                                     cttcttCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcg                                                         >  1:402134/1‑139 (MQ=255)
                                                                     cttcttCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcg                                                         >  2:405582/1‑139 (MQ=255)
                                                                              cgcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAg                                                >  1:318793/1‑139 (MQ=255)
                                                                                      gggAGGAAGCGGGAGGCACACCCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgca                                        <  2:321518/137‑1 (MQ=255)
                                                                                                        cccccGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCtt                      <  1:249568/136‑1 (MQ=255)
                                                                                                                        ggCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGc      <  2:318793/139‑1 (MQ=255)
                                                                                                                         gCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgcgc                                         <  1:346193/103‑1 (MQ=255)
                                                                                                                         gCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgcgc                                         >  2:346193/1‑103 (MQ=255)
                                                                                                                              ccACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCggg  <  1:405582/137‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                   
CCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGG  >  NZ_CP009273/2883895‑2884157

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: