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% 65.5 / NA 19D411G (GAC→GGC) cysJNADPH‑dependent assimilatory sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/11);  total (8/11)

CCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGG  >  NZ_CP009273/2883895‑2884157
                                                                                                                                   |                                                                                                                                   
ccGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACg                                                                                                                              <  1:219794/139‑1 (MQ=255)
 cGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACg                                                                                                                              >  2:218354/1‑138 (MQ=255)
             aTTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACg                                                                                                                 >  2:557893/1‑139 (MQ=255)
             aTTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTCCCACGCGGTCAGCGAGGAAGCTGGAGGCACCCCCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACg                                                                                                                 <  2:15958/139‑1 (MQ=255)
               tggctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGt                                                                                                               >  1:274774/1‑139 (MQ=255)
               tggctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGt                                                                                                               >  1:459946/1‑139 (MQ=255)
                ggctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGt                                                                                                               <  2:226908/138‑1 (MQ=255)
                ggctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGt                                                                                                               <  1:218354/138‑1 (MQ=255)
                                ttATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtct                                                                                              <  2:104641/139‑1 (MQ=255)
                                          gTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGc                                                                                    >  2:290537/1‑139 (MQ=255)
                                                        cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGTTGGAGGCACCCCCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                                                      <  1:340509/139‑1 (MQ=255)
                                                        cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCCCCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                                                      <  1:548951/139‑1 (MQ=255)
                                                        cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                                                      <  1:529725/139‑1 (MQ=255)
                                                          gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATg                                                                    <  1:340734/139‑1 (MQ=255)
                                                                                  gTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGcac                                                                                                                           >  2:266399/1‑60 (MQ=255)
                                                                                  gTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGcac                                                                                                                           <  1:266399/60‑1 (MQ=255)
                                                                                         aGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcag                                       >  2:7862/1‑137 (MQ=255)
                                                                                                                              ccACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCggg  <  2:287860/137‑1 (MQ=255)
                                                                                                                              ccACGCCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCggg  >  2:8761/1‑137 (MQ=255)
                                                                                                                                   |                                                                                                                                   
CCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGG  >  NZ_CP009273/2883895‑2884157

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