Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,444,344 C→G 8.2% D178H (GAC→CAC)  nicF ← maleamate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,444,3440CG 8.2% 148.5 / 5.1 49D178H (GAC→CAC) nicFmaleamate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base G (2/2);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

CGGTCCTTGATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGTGCATGCAGCAGCGGGGCCAGGCTGGTGGCG  >  NC_002947/4444196‑4444492
                                                                                                                                                    |                                                                                                                                                    
cGGTCCTTGATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCa                                                                                                                                                     >  2:650520/1‑150 (MQ=255)
      ttGATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGTCGGTGACCgcca                                                                                                                                               >  2:682145/1‑147 (MQ=255)
      ttGATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCNaca                                                                                                                                               <  2:704408/150‑1 (MQ=255)
      ttGATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGaca                                                                                                                                               >  2:841309/1‑150 (MQ=255)
      ttGATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGaca                                                                                                                                               <  2:727479/150‑1 (MQ=255)
        gATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGTCGGTGACCGccaca                                                                                                                                             <  1:682145/150‑1 (MQ=255)
             ttATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTcgc                                                                                                                                        >  2:296358/1‑150 (MQ=255)
             ttATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGATGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTcgc                                                                                                                                        <  1:120765/150‑1 (MQ=255)
              tATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGGTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTCCTGTCGATGTCGCACCGAGGATCCGCGGGCGGGCCGCCTTGGCGGTCACCGACACACTcgcg                                                                                                                                       >  2:242471/1‑150 (MQ=255)
                ttGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGca                                                                                                                                     <  1:903242/150‑1 (MQ=255)
                           cTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGtt                                                                                                                          <  1:296358/150‑1 (MQ=255)
                               ccGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGaa                                                                                                                      >  2:1110927/1‑150 (MQ=255)
                                 ggCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAAc                                                                                                                    <  2:993132/150‑1 (MQ=255)
                                          ccGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGtgctgc                                                                                                           <  1:242471/150‑1 (MQ=255)
                                            ggCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTGGAACAGGTTGGCTTCGTGCGGGTCGGTGTGGCGGTCACCGGCACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGAAt                                                                                                         >  1:209726/1‑150 (MQ=255)
                                            ggCCAGGTGCCTGAGTTGCTGCATGGCATCCCGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCAt                                                                                                         <  2:190243/150‑1 (MQ=255)
                                                       tGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCg                                                                                              >  1:211414/1‑150 (MQ=255)
                                                       tGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATACACCg                                                                                              >  1:864232/1‑150 (MQ=255)
                                                               tgcATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGgcg                                                                                      >  1:681065/1‑150 (MQ=255)
                                                               tgcATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATACACCGCACTGgcg                                                                                      <  2:864232/150‑1 (MQ=255)
                                                                gcATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGgcgc                                                                                     <  2:211414/150‑1 (MQ=255)
                                                                 cATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGcgcg                                                                                    >  1:351830/1‑150 (MQ=255)
                                                                   tGGCATCCTGCCGGTTGCCCACATCCCCGTACTCGCTGCCGCTGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGa                                                                                  <  2:181991/150‑1 (MQ=255)
                                                                   tGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGa                                                                                  >  1:520440/1‑150 (MQ=255)
                                                                   tGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGa                                                                                  <  1:414964/150‑1 (MQ=255)
                                                                     gCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATg                                                                                >  2:74924/1‑150 (MQ=255)
                                                                       aTCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGCGACGCGTTCTTCATGGCATCCACCGCACTGGCGCGGATGCa                                                                              >  2:1112451/1‑150 (MQ=255)
                                                                        tCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAg                                                                             <  1:841309/150‑1 (MQ=255)
                                                                        tCCTGCCGCGTGATGACAGCACCTGACGTTCTTTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAg                                                                             <  2:238091/150‑1 (MQ=255)
                                                                         ccTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGc                                                                            <  1:514380/150‑1 (MQ=255)
                                                                              cGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGTTGGGTAGCATCCACCGCACTGGCGCGGATGCAGCCGCTg                                                                       <  1:1112451/150‑1 (MQ=255)
                                                                                 gTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTc                                                                    <  2:292982/150‑1 (MQ=255)
                                                                                 gTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTc                                                                    <  1:74924/150‑1 (MQ=255)
                                                                                 gTGACGACATCACCGTACTTGCTGACGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTAGTGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTc                                                                    <  2:564312/150‑1 (MQ=255)
                                                                                    acgacATCACCGTACTTGCTGTCGCTGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTagag                                                                 <  1:1007209/150‑1 (MQ=255)
                                                                                                           gatgtcgaACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATc                                          <  2:399450/150‑1 (MQ=255)
                                                                                                           aatgtcgaACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATc                                          <  2:407079/149‑1 (MQ=255)
                                                                                                                cgaACAGGTTGGCTTCGGGCGGGTCGCTGTGTCGTTCAGCGACACACTCGCGCACCACGGTGGTTAGTTTTGCTTTCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCacac                                     <  1:748664/150‑1 (MQ=255)
                                                                                                                       gTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGtgc                              <  2:520440/150‑1 (MQ=255)
                                                                                                                        ttGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGtgca                             >  1:1074214/1‑150 (MQ=255)
                                                                                                                           gCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGtgcatgc                          <  1:650520/150‑1 (MQ=255)
                                                                                                                             ttCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGTGCATgcag                        >  2:842799/1‑150 (MQ=255)
                                                                                                                                 tGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGTGCATGCagcagc                    <  2:351830/150‑1 (MQ=255)
                                                                                                                                       tCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGTGCATGCAGCAGCGGGGcc              >  1:963018/1‑150 (MQ=255)
                                                                                                                                          gtgtgTCGGTGACCGCCACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTg                                                                       >  1:892387/2‑90 (MQ=255)
                                                                                                                                          gtgtgTCGGTGACCGCCACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTg                                                                       <  2:892387/89‑1 (MQ=255)
                                                                                                                                           tgtgGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTtgtg                               >  2:720045/1‑129 (MQ=255)
                                                                                                                                           tgtgGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTtgtg                               <  1:720045/129‑1 (MQ=255)
                                                                                                                                                   tCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGTGCATGCAGCAGCGGGGCCAGGCTGGTGGCg  >  2:223400/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGTCCTTGATTGTTATTGGGTAGGCGCTGCCCGGCGTTCAGCCGGCCAGGTGCCTGAGTTGCTGCATGGCATCCTGCCGGGTGACGACATCACCGTACTTGCTGTCGATGTCGAACAGGTTGGCTTCGTGCGGGTCGCTGTGGCGGTCACCGACACACTCGCGCACCACGATGGTTCGAAACCCGTGCTGCATGGCATCCACCGCACTGGCGCGGATGCAGCCGCTGGTCGAGCAGCCAGCCAGCACCACGGTATCCACACCTTGTGCATGCAGCAGCGGGGCCAGGCTGGTGGCG  >  NC_002947/4444196‑4444492

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