Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 718,639 G→T 7.8% F426L (TTC→TTA)  PP_0613 ← amidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947718,6390GT 7.8% 157.0 / 3.7 51F426L (TTC→TTA) PP_0613amidase family protein
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base T (2/2);  total (26/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 = 7.16e-01

CAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATAC  >  NC_002947/718501‑718785
                                                                                                                                          |                                                                                                                                                   
cAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCCTGCGGGTTCACCCGcgat                                                                                                                                           <  2:288158/149‑4 (MQ=255)
 aGATTTCAGACAGTAATCGGTGCTTG‑AGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACAGGCAACCCCTCGCGGCTGAACCCGC‑Acgg                                                                                                                                         >  8:123643/1‑148 (MQ=255)
         gCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAgg                                                                                                                                 >  3:146635/1‑149 (MQ=255)
           cAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGgcc                                                                                                                               >  2:313948/1‑149 (MQ=255)
                aTCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGccgccg                                                                                                                           <  4:259957/148‑1 (MQ=255)
                    gTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGt                                                                                                                      >  7:311082/1‑149 (MQ=255)
                     tGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGtt                                                                                                                     <  3:327505/149‑1 (MQ=255)
                      gCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTg                                                                                                                    >  3:148195/1‑149 (MQ=255)
                       cTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTgg                                                                                                                   <  4:253433/149‑1 (MQ=255)
                               tGCCGGCTCGGATAATGCTGCTCATACACCTTGCACCCCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGGCAGGtt                                                                                                           <  8:311082/149‑1 (MQ=255)
                                      tCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGgcc                                                                                                                               <  3:102400/122‑1 (MQ=255)
                                      tCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGgcc                                                                                                                               >  4:102400/1‑122 (MQ=255)
                                        ggATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATa                                                                                                  <  4:178900/149‑1 (MQ=255)
                                              tgctgcTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTaa                                                                                            >  7:30466/1‑149 (MQ=255)
                                                     cATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTc                                                                                     >  1:34169/1‑149 (MQ=255)
                                                      aTACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACg                                                                                                      >  7:47714/1‑131 (MQ=255)
                                                      aTACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACg                                                                                                      <  8:47714/131‑1 (MQ=255)
                                                         cacCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACt                                                                                 <  7:18545/149‑1 (MQ=255)
                                                          acCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAAc                                                                                           >  4:308812/1‑138 (MQ=255)
                                                          acCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAAc                                                                                           <  3:308812/138‑1 (MQ=255)
                                                          acCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACtc                                                                                >  8:272288/1‑149 (MQ=255)
                                                           ccTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACtc                                                                                >  2:177418/1‑148 (MQ=255)
                                                             ttGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACtccat                                                                             >  4:292354/1‑149 (MQ=255)
                                                                 acacCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGAAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTccatccat                                                                         >  1:80558/1‑149 (MQ=255)
                                                                   acCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTg                                                                       <  4:146635/149‑1 (MQ=255)
                                                                   acCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTg                                                                       >  8:204584/1‑149 (MQ=255)
                                                                     ccGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTgcg                                                                     <  1:177418/149‑1 (MQ=255)
                                                                            aCTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATccc                                                              <  7:272288/149‑1 (MQ=255)
                                                                             cTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATcccc                                                             <  3:292354/149‑1 (MQ=255)
                                                                                       aGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAggcg                                                   >  8:128547/1‑149 (MQ=255)
                                                                                       aGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAggcg                                                   >  5:367307/1‑149 (MQ=255)
                                                                                         cAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACtccatc                                                                            <  7:93848/122‑1 (MQ=255)
                                                                                         cAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACtccatc                                                                            >  8:93848/1‑122 (MQ=255)
                                                                                          aaaGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCa                                                >  7:148937/1‑149 (MQ=255)
                                                                                           aaGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCAc                                               <  1:266392/149‑1 (MQ=255)
                                                                                              cgcCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGtt                                            <  6:367307/149‑1 (MQ=255)
                                                                                                  cGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTgtgt                                        >  2:359497/1‑149 (MQ=255)
                                                                                                        ccaccTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCt                                  <  1:359497/149‑1 (MQ=255)
                                                                                                          accTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTc                                <  2:34169/149‑1 (MQ=255)
                                                                                                              gCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACTGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAAc                            >  8:132795/1‑149 (MQ=255)
                                                                                                                  gCCCACCGGTAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGtt                                            <  8:295302/129‑1 (MQ=255)
                                                                                                                  gCCCACCGGAAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGtt                                            >  7:295302/1‑129 (MQ=255)
                                                                                                                    ccACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGccacc                      <  7:128547/149‑1 (MQ=255)
                                                                                                                    ccACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGccac                       >  8:16604/1‑148 (MQ=255)
                                                                                                                              ccccTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGgcagca            >  1:228443/1‑149 (MQ=255)
                                                                                                                                 cctgcgGGTTCACCCGCGAGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGc                                   <  3:56654/121‑1 (MQ=255)
                                                                                                                                 cctgcgGGTTCACCCGCGAGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGc                                   >  4:56654/4‑124 (MQ=255)
                                                                                                                                 cctgcgGGTTCACCCGCGAGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATc          >  6:368412/4‑149 (MQ=255)
                                                                                                                                    gcgGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGc       <  4:148195/148‑1 (MQ=255)
                                                                                                                                     cgGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGTGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGa     <  7:123643/149‑1 (MQ=255)
                                                                                                                                        gTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATAc  <  1:313948/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
CAGATTTCAGCCAGTAATCGGTGCTTGCAGGTGCCGGCTCGGATAATGCTGCTCATACACCTTGCACACCCGCAGCACTTGCTCATCAGCAAAGCGCCCGGCCACCACCTGCAGGCCCACCGGCAACCCCTCGCGGGTGAACCCGC‑AGGGTACCGAGGCCGCCGGCTGTTGGGTCAGGTTGAACGGATAGCTTAACGGCGTCCACTCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATAC  >  NC_002947/718501‑718785

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: ‑