Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,148,546 G→A 8.5% A162T (GCC→ACC)  msuE → FMN reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,148,5460GA 8.5% 135.0 / 3.2 47A162T (GCC→ACC) msuEFMN reductase
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base A (2/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCG  >  NC_002947/3148408‑3148692
                                                                                                                                           |                                                                                                                                                  
ccGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCgg                                                                                                                                          <  2:299408/149‑1 (MQ=255)
ccGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCgg                                                                                                                                          >  4:41719/1‑149 (MQ=255)
  gCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCggtg                                                                                                                                        >  2:101027/1‑149 (MQ=255)
     tttGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTa                                                                                                                                     <  2:88717/149‑1 (MQ=255)
            gCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTagcgagcg                                                                                                                              <  2:94402/149‑1 (MQ=255)
              cATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAgcgagcga                                                                                                                             >  3:151553/1‑148 (MQ=255)
              cATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAgcgagcga                                                                                                                             >  3:77455/1‑148 (MQ=255)
              cATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGAc                                                                                                                            <  4:151553/149‑1 (MQ=255)
                    gCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCt                                                                                                                      >  1:382858/1‑149 (MQ=255)
                     cAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCtt                                                                                                                     >  2:145158/1‑149 (MQ=255)
                           gACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTtgatga                                                                                                                <  1:145158/148‑1 (MQ=255)
                           gACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTgatgat                                                                                                               <  4:77455/149‑1 (MQ=255)
                            aCCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatc                                                                                                              >  4:139886/1‑149 (MQ=255)
                             ccTGCTACTGGTCGGCAGCCCGGGGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcg                                                                                                             <  7:136710/149‑1 (MQ=255)
                             ccTGCTACTGGGCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATgatcg                                                                                                             <  7:142490/149‑1 (MQ=255)
                                    cTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCa                                                                                                      <  3:41719/149‑1 (MQ=255)
                                     tGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAg                                                                                                     >  8:6306/1‑149 (MQ=255)
                                       gTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGtt                                                                                                   <  6:184510/149‑1 (MQ=255)
                                       gTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGtt                                                                                                   >  5:184510/1‑149 (MQ=255)
                                                   gtgtATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTAtt                                                                                       <  7:6306/149‑1 (MQ=255)
                                                        tggcgcgTCCTACCCGGTTTTGTTCAAGCACCTGTTCGCCCTGGGCTATCACCCGTCGCTCAAGGCCGTCCCCGT‑GGTGCTTGCCGCCACCGGTGTTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCt                                                                                   <  1:101027/146‑1 (MQ=255)
                                                        tCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCtt                                                                                  >  3:145396/1‑149 (MQ=255)
                                                                              ttCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGc                                                             <  3:139886/148‑1 (MQ=255)
                                                                                cAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCg                                                           >  3:139166/1‑148 (MQ=255)
                                                                                 aaGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTc                                                                                      <  5:218492/120‑1 (MQ=255)
                                                                                 aaGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTc                                                                                      >  6:218492/1‑120 (MQ=255)
                                                                                    cACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGTGGGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACg                                                       >  1:96547/1‑149 (MQ=255)
                                                                                     aCCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGc                                                     <  4:49471/149‑1 (MQ=255)
                                                                                     aCCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGc                                                     <  2:64928/149‑1 (MQ=255)
                                                                                       cTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCt                                                   >  5:387333/1‑149 (MQ=255)
                                                                                              aCCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGcc                                            <  4:145396/149‑1 (MQ=255)
                                                                                               ccTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCa                                           <  6:387333/149‑1 (MQ=255)
                                                                                                   gTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTACCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAgtgt                                       >  4:348897/1‑149 (MQ=255)
                                                                                                   gTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTACCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAgtgt                                       >  4:354744/1‑149 (MQ=255)
                                                                                                    tCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTc                                      <  2:382858/149‑1 (MQ=255)
                                                                                                     cGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTc                                      >  5:244842/1‑148 (MQ=255)
                                                                                                      gATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTACCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGa                                    <  3:354744/149‑1 (MQ=255)
                                                                                                      gATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTACCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGa                                    <  3:348897/149‑1 (MQ=255)
                                                                                                                  cTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGc                                             >  5:390364/1‑128 (MQ=255)
                                                                                                                  cTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGc                                             <  6:390364/128‑1 (MQ=255)
                                                                                                                          cGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGcc                >  4:74506/1‑149 (MQ=255)
                                                                                                                           gTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCt               >  5:70814/1‑149 (MQ=255)
                                                                                                                             gCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTgg             >  4:184907/1‑149 (MQ=255)
                                                                                                                              cccGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTgg             <  6:244842/148‑1 (MQ=255)
                                                                                                                                  t‑ggtgCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGAc        <  6:70814/149‑1 (MQ=255)
                                                                                                                                     gtgCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGAcc       <  2:96547/148‑1 (MQ=255)
                                                                                                                                         ttGCCGCCACCGGAGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTgcg  >  8:285832/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CCGCATTTGCAGGCCATCGAGCAAGCCGACCTGCTACTGGTCGGCAGCCCGGTGTATCGCGCGTCCTACACCGGTTTGTTCAAGCACCTGTTCGACCTGGTCGATCACCAGTCGCTCAAGGGCGTGCCCGT‑GGTGCTTGCCGCCACCGGTGGTAGCGAGCGACATGCCTTGATGATCGATCACCAGTTGCGCCCGCTATTCGCTTTCTTCCAGGCCCACACCCTGCCGTACGGCCTGTACGCCAGTGTCGAGTCTTTCGACGATCAGCGCCTGGCTGACCCTGCG  >  NC_002947/3148408‑3148692

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