Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,727,292 G→A 8.7% A212T (GCA→ACA)  nemA → chromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,727,2920GA 8.7% 138.1 / 3.2 46A212T (GCA→ACA) nemAchromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base A (2/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCAGGCCATTGCTAACGCGCGTGAAGCCGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAA  >  NC_000913/1727152‑1727435
                                                                                                                                            |                                                                                                                                               
gTCAGGCCATTGCTAACGCGCGTGAAGCCGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTg                                                                                                                                         >  1:4562/1‑149 (MQ=255)
gTCAGGCCATTGCTAACGCGCGTGAAGCCGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTg                                                                                                                                         <  2:30457/149‑1 (MQ=255)
         ttGCTAATGCGCGTGAAGCCGGTTTTGATCTGGTAGAGCTCCTCTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGaa                                                                                                                                <  1:208976/149‑1 (MQ=255)
           gCTAACGCGCGTGAAGCCGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGt                                                                                                                              >  1:134557/1‑149 (MQ=255)
                cgcgcgTGAAGCCGGTTTTGATCTGGTAGAGCTCCACTCTGGTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCACACGTTTGGTACTGGAAGTGGTc                                                                                                                          >  3:53220/1‑148 (MQ=255)
                          gCCGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTgg                                                                                                                            >  2:91109/1‑136 (MQ=255)
                          gCCGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTgg                                                                                                                            <  1:91109/136‑1 (MQ=255)
                            cGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTg                                                                                                             <  2:278282/149‑1 (MQ=255)
                                       tGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATgggg                                                                                                   <  3:318221/148‑1 (MQ=255)
                                        ggTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGt                                                                                                  >  4:200310/1‑148 (MQ=255)
                                             agCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGAt                                                                                            <  2:134557/149‑1 (MQ=255)
                                                 ccACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGc                                                                                         >  3:132195/1‑148 (MQ=255)
                                                  cACTCTGGTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCACACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCa                                                                                        <  4:53220/148‑1 (MQ=255)
                                                  cACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCACACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCAt                                                                                       >  2:213883/1‑149 (MQ=255)
                                                   aCTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTTGATGCCGGGATTGAAGAATGGGGTGCCGATCGCAt                                                                                       >  3:347915/1‑148 (MQ=255)
                                                          tCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTc                                                                               >  2:135859/1‑149 (MQ=255)
                                                            aCGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTcgc                                                                             >  1:88817/1‑149 (MQ=255)
                                                             cGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTcgc                                                                             <  3:147225/148‑1 (MQ=255)
                                                                   tttGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGcc                                                                                               <  2:188018/124‑1 (MQ=255)
                                                                   tttGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGcc                                                                                               >  1:188018/1‑124 (MQ=255)
                                                                    ttGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCAc                                                                      >  4:240588/1‑148 (MQ=255)
                                                                           aTCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGTTGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCgg                                                               >  4:230749/1‑148 (MQ=255)
                                                                           aTCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCgg                                                               <  3:240588/148‑1 (MQ=255)
                                                                           aTCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGt                                                              >  2:314992/1‑149 (MQ=255)
                                                                           aTCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGt                                                              <  1:167563/149‑1 (MQ=255)
                                                                            tCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTa                                                             <  1:135859/149‑1 (MQ=255)
                                                                            tCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTa                                                             >  2:88364/1‑149 (MQ=255)
                                                                            tCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTa                                                             >  2:311903/1‑149 (MQ=255)
                                                                               gTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACt                                                           >  4:289816/1‑148 (MQ=255)
                                                                                ttCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACtt                                                          >  4:312994/1‑148 (MQ=255)
                                                                                ttCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACtt                                                          >  4:105345/1‑148 (MQ=255)
                                                                                     ttCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCggg                                                                                                                 >  3:208193/1‑88 (MQ=255)
                                                                                     ttCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCggg                                                                                                                 <  4:208193/88‑1 (MQ=255)
                                                                                                 ataCCATCGTACCGATCAGTACGGCGGCAGCGTGCAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAAc                                         <  3:330564/146‑1 (MQ=255)
                                                                                                        cGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGa                                  >  4:350128/1‑148 (MQ=255)
                                                                                                         gTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTTGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGaa                                 <  4:347915/148‑1 (MQ=255)
                                                                                                         gTACCGATCAGGACGGCGGCAGCGTGGAAAATCGCACACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAAt                                <  1:213883/149‑1 (MQ=255)
                                                                                                           aCCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATg                               >  4:173750/1‑148 (MQ=255)
                                                                                                              gATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTg                                                                                                                                         >  2:373663/1‑39 (MQ=255)
                                                                                                              gATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTg                                                                                                                                         <  1:373663/39‑1 (MQ=255)
                                                                                                               aTCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATgaagaa                          <  1:88364/149‑1 (MQ=255)
                                                                                                                tCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATgaagaa                          <  3:350128/148‑1 (MQ=255)
                                                                                                                  aGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGc                        <  4:16431/148‑1 (MQ=255)
                                                                                                                                  ggAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTg        <  3:98411/148‑1 (MQ=255)
                                                                                                                                       aTCGCGCACGTTTGTTGCTGGAAGTGGTCGATGCGGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGATGAAGCCGATGCACTATATCTGATTGa   >  4:250521/1‑148 (MQ=255)
                                                                                                                                        tCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGaa  >  3:186591/1‑148 (MQ=255)
                                                                                                                                          gcgcACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGaa  <  1:252474/146‑1 (MQ=255)
                                                                                                                                          gcgcACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGaa  >  2:252474/1‑146 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GTCAGGCCATTGCTAACGCGCGTGAAGCCGGTTTTGATCTGGTAGAGCTCCACTCTGCTCACGGTTATTTGCTGCATCAGTTCCTTTCTCCTTCTTCAAACCATCGTACCGATCAGTACGGCGGCAGCGTGGAAAATCGCGCACGTTTGGTACTGGAAGTGGTCGATGCCGGGATTGAAGAATGGGGTGCCGATCGCATTGGCATTCGCGTTTCACCAATCGGTACTTTCCAGAACACAGATAACGGCCCGAATGAAGAAGCCGATGCACTGTATCTGATTGAA  >  NC_000913/1727152‑1727435

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