Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,853,466 T→C 10.8% L63P (CTG→CCG)  mazG → modified nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,853,4660TC10.8% 95.9 / 2.7 37L63P (CTG→CCG) mazGmodified nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base C (2/2);  total (17/20)
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.74e-01

CGCGACCCGCAGTACGGCTGCCCGTGGGACCTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCTGTACGCGCCGGTGGACACG  >  NC_002947/1853324‑1853614
                                                                                                                                              |                                                                                                                                                    
cgcgACCCGCAGTACGGCTGCCCGTGGGACCTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGt                                                                                                                                                 <  1:83071/148‑1 (MQ=255)
      ccGCAGTACGGCTGCCCGTGGGACCTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAggt                                                                                                                                          >  8:42722/1‑149 (MQ=255)
      ccGCAGTACGGCTGCCCGTGGGAACTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAggt                                                                                                                                          >  7:168552/1‑149 (MQ=255)
        gCAGTACGGCTGCCCGTGGGACCTAAAGCAGAACTACGCCAGCATCATCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAggtgg                                                                                                                                        >  6:158466/1‑149 (MQ=255)
          aGTACGGCTGCCCGTGGGACCTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTc                                                                                                                                      <  4:256736/149‑1 (MQ=255)
                     ccGTGGGACCTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTAcagcca                                                                                                                           <  7:249901/149‑1 (MQ=255)
                            aCCTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGcc                                                                                                                    >  3:1717/1‑149 (MQ=255)
                             ccTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGccc                                                                                                                   <  4:1717/149‑1 (MQ=255)
                                aaaGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCcgc                                                                                                                 >  1:200373/1‑148 (MQ=255)
                                   gCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCCGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCcgcg                                                                                                                >  2:263210/1‑146 (MQ=255)
                                   gCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCCGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCcgcg                                                                                                                >  2:267464/1‑146 (MQ=255)
                                       aaCTACGCCAGCATCGTCCTGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAgg                                                                                                         >  1:266693/1‑149 (MQ=255)
                                         cTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGc                                                                                                       <  1:68889/149‑1 (MQ=255)
                                                 gCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAg                                                                                               <  2:200373/149‑1 (MQ=255)
                                                  cATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGt                                                                                              <  7:42722/149‑1 (MQ=255)
                                                    tCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCCGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTc                                                                                            <  1:267464/149‑1 (MQ=255)
                                                    tCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCCGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTc                                                                                            <  1:263210/149‑1 (MQ=255)
                                                             acacCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTg                                                                                   <  2:266693/149‑1 (MQ=255)
                                                                  aTCGAAGACTCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGGACCTGCAATGCGAGCTGGTCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGa                                                                              <  1:241374/149‑1 (MQ=255)
                                                                     gAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACa                                                                            <  5:202991/148‑1 (MQ=255)
                                                                           gCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCAc                                                                     >  1:222898/1‑149 (MQ=255)
                                                                                cGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCGGCGAAGAAGGGCGTTTCGAGTTCGATg                                                                                        <  4:107646/125‑1 (MQ=255)
                                                                                cGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATg                                                                                        >  3:107646/1‑125 (MQ=255)
                                                                                       gCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCAt                                                         >  7:143600/1‑149 (MQ=255)
                                                                                           acacCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATccg                                                      <  2:222898/148‑1 (MQ=255)
                                                                                            caacATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGCGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATacgcc                                                    >  3:240894/4‑149 (MQ=255)
                                                                                             acCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATccgccg                                                   >  2:74836/1‑149 (MQ=255)
                                                                                                       gcggcgATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCAc                                         <  5:158466/149‑1 (MQ=255)
                                                                                                         ggcgATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCAcc                                                                    >  5:148616/1‑120 (MQ=255)
                                                                                                         ggcgATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCAcc                                                                    <  6:148616/120‑1 (MQ=255)
                                                                                                              tttCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGATGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTc                                  >  2:216824/1‑149 (MQ=255)
                                                                                                                       ccTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGc                                                                 >  7:198240/1‑109 (MQ=255)
                                                                                                                       ccTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGc                                                                 <  8:198240/109‑1 (MQ=255)
                                                                                                                        cTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGc                          >  3:106195/1‑147 (MQ=255)
                                                                                                                                    cTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCTGTACgcgc             <  8:143600/148‑1 (MQ=255)
                                                                                                                                        gCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCTGTACGCGCCGGTg        <  1:74836/149‑1 (MQ=255)
                                                                                                                                           aCCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCTGTACGCGCCGGTGGac     <  4:106195/149‑1 (MQ=255)
                                                                                                                                              tgctgTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCt                     >  2:239732/1‑130 (MQ=255)
                                                                                                                                              tgctgTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCt                     <  1:239732/130‑1 (MQ=255)
                                                                                                                                              tgctgTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCTGTACGCGCCGGTGAACACg  >  7:170692/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CGCGACCCGCAGTACGGCTGCCCGTGGGACCTAAAGCAGAACTACGCCAGCATCGTCCCGCACACCATCGAAGAGGCCTACGAGGTCGCTGACACCATTGAGCGCGGCGATTTCGAGCACCTGCAAGGCGAGCTGGGCGACCTGCTGTTCCAGGTGGTCTACTACAGCCAGCTGGCCCGCGAAGAAGGGCGTTTCGAGTTCGATGGCGTGGTCGACAGCATCACCCGCAAGCTCATCCGCCGCCACCCGCACGTGTTTCCGACTGGCGAGCTGTACGCGCCGGTGGACACG  >  NC_002947/1853324‑1853614

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