Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,727,494 (G)6→7 coding (836/1098 nt) nemA → N‑ethylmaleimide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,727,4941.G85.7% 76.6 / 6.9 28coding (836/1098 nt)nemAN‑ethylmaleimide reductase
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (7/17);  total (9/19)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGG‑TG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGAT  >  NC_000913/1727250‑1727741
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
aaCCATCGTACCGATCAGTACGGCGGCAGCGTGGCAAATCGCGCACGTTTGGTCCTGGCAGCGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCCCCAATCGGTACTTTCCAGAACACAAATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAAACAGATTGGG‑CG‑G‑GGGGGTG‑AAc                                                                                                                                                                                                                                                    >  2:202140/1‑251 (MQ=255)
aaCCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGGCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATACACTGTATCTGATTGAACCACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AAc                                                                                                                                                                                                                                                    >  1:226002/1‑251 (MQ=255)
              tCAGTACGGCGGCAGCGTGGAAAATCGCGCCCGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGCTCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGG‑GGGAACCGTATACTGATGc                                                                                                                                                                                                                                       >  2:280774/1‑250 (MQ=255)
              tCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGc                                                                                                                                                                                                                                       <  1:280774/250‑1 (MQ=255)
                 gTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTc                                                                                                                                                                                                                                   <  1:35456/251‑1 (MQ=255)
                                ggAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTAcg                                                                                                                                                                                                                     <  1:83286/250‑1 (MQ=255)
                                     aTCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATGGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTAGCGGCATATGTCAGAACCAGATTGGG‑CG‑GTGGGG‑TG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCg                                                                                                                                                                                                               <  1:66633/251‑1 (MQ=255)
                                                  ggTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTGCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTAGCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTTCATATGTCAGAACCAGATGGGG‑CG‑G‑GGGGGTT‑AATTGTATACTGATGGGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCg                                                                                                                                                                                                  <  1:409274/251‑1 (MQ=255)
                                                     aCTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATGGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAtt                                                                                                                                                                                               <  1:368535/251‑1 (MQ=255)
                                                     aCTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGAtt                                                                                                                                                                                               <  1:100724/251‑1 (MQ=255)
                                                                      ccGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATa                                                                                                                                                                              <  1:214368/251‑1 (MQ=255)
                                                                         ggATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTAGCTGCATATGTCAGAACCAGATTGGT‑CG‑G‑GGGGGTG‑ATCCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCTCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATacac                                                                                                                                                                           <  1:177871/251‑1 (MQ=255)
                                                                              gaagaaTGGGGAGCCGATCGCATTGGACTTCCCGTCTCAGCAATCGGTACTTTCCAGAAACCAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAAACACTGGGCAAACGCGGCATTGCTTATCTGCATATGTCAGAAACCGATTGGGAGG‑G‑GGGG‑TG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGAGCCCGTTTCCACGGACCGATTATCGGCGCAGGGGCATACACAGt                                                                                                                                                                        >  2:278087/1‑249 (MQ=255)
                                                                                              aTCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGCTTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACACTAGAAAAAGCTGAAACGCt                                                                                                                                                      >  2:334301/1‑251 (MQ=255)
                                                                                                      ggCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGc                                                                                                                                               <  2:371347/250‑1 (MQ=255)
                                                                                                       gCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTAGCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCaa                                                                                                                                             <  1:431638/251‑1 (MQ=255)
                                                                                                                                 ttCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCTATGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACCCAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAGTTGATGCGGTGGCATTTg                                                                                                                   >  2:389619/1‑251 (MQ=255)
                                                                                                                                   ccAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAGCTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCGGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTTGCATTTgg                                                                                                                  <  1:215216/250‑1 (MQ=255)
                                                                                                                                   ccAGAACACAGATAACGGCCCGAATGAAGAAGCCGATCCTCTTTATCTCCTTGAACAACTGGGTAAACGCCGCATTGCTTATCTGCATATGTCCGAACCCGATTGGG‑CG‑G‑GGGGGCC‑CACCGTATACTGCTGCGTTCCGCGACAACCTACGCGCCCCTCTCCCCGGTCCGCCCCCCGGCGCAGGTGCATACACAGTAGCATCAGCTGAAACGCTGCTCTGCAAAGGGCTAATTGACGCGCTGGCacctgct                                                                                                                 >  1:85796/1‑245 (MQ=255)
                                                                                                                                                   ggCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCTGTGGCATATGGTCGTGACTGGATTGCGa                                                                                                 <  1:105968/251‑1 (MQ=255)
                                                                                                                                                   ggCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGATTGCg                                                                                                  >  1:159203/1‑250 (MQ=255)
                                                                                                                                                                         cACTTTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGT‑CGTG‑GGGG‑TT‑AACCGTATACCGATTCGTTCCGCGAAAAAGGACGCGCCCGTTTCCACGGTCCTATTATCGGCGCAGGTGCATACACATTAGAAAAAGCTGAAACGCTGATCGGCAAAGTGTTAATTGATGCGTTGGCATTTGTTCGTGACTGTATTGCGAACCCGGATCTGGTCGCCCGCtt                                                                           <  1:374962/251‑1 (MQ=255)
                                                                                                                                                                                         acaacTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCAGTTTCCGCGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGCCGTGACTCTATGGCGAACCCGGATCTGTTCGCCCGGTTGCAGCGCAAAGCTGAg                                                           <  1:134455/251‑1 (MQ=255)
                                                                                                                                                                                            acTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTGCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCGGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGAGTGCTAACCCGGATCTGGTCGCCCGCGTGCATCGCAAAGCTGAGCtt                                                        <  2:159203/251‑1 (MQ=255)
                                                                                                                                                                                                        cgggaTTGCTTAACTGCAAATGGCAGAAACAAGTTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGt                                            <  1:202642/247‑1 (MQ=255)
                                                                                                                                                                                                                   aTCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGttt                                 >  1:130128/1‑251 (MQ=255)
                                                                                                                                                                                                                         atatGTCAGAACCAGATGGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTATAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGTTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCGTTCAGCGCAAAGCTGAGCTTAACCCACAGCGGGCCGAAAGTTTCTACgg                           <  2:61681/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gATTGGG‑CG‑G‑GGGGGTG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAg             <  2:82952/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    gggg‑TG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTGTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGAt  <  1:244302/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
AACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAACAACTGGGTAAACGCGGCATTGCTTATCTGCATATGTCAGAACCAGATTGGG‑CG‑G‑GGGG‑TG‑AACCGTATACTGATGCGTTCCGCGAAAAAGTACGCGCCCGTTTCCACGGTCCGATTATCGGCGCAGGTGCATACACAGTAGAAAAAGCTGAAACGCTGATCGGCAAAGGGTTAATTGATGCGGTGGCATTTGGTCGTGACTGGATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGAT  >  NC_000913/1727250‑1727741

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