Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,210,278 C→A 7.7% V160V (GTC→GTA)  PP_1953 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,210,2780CA 7.7% 165.3 / 3.4 52V160V (GTC→GTA) PP_1953short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (26/22);  new base A (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

AAGAGCGGGCGGGCGCCCACTTTCGACTTCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCTTACAACATTCGTGTGAATGCCAT  >  NC_002947/2210129‑2210412
                                                                                                                                                     |                                                                                                                                      
aaGAGCGGGCGGGCGCCCACTTTCGACTTCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTc                                                                                                                                        <  1:357539/150‑1 (MQ=255)
                ccACTTTCGACTTCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTgg                                                                                                                        <  1:465971/150‑1 (MQ=255)
                ccACTTTCGACTTCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTgg                                                                                                                        >  1:625210/1‑150 (MQ=255)
                   cTTTCGACTTCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGt                                                                                                                          <  2:1278575/145‑1 (MQ=255)
                   cTTTCGACTTCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGt                                                                                                                          >  1:1278575/1‑145 (MQ=255)
                         aCTTCCTGGCCGCTACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTcc                                                                                                               >  2:419730/1‑150 (MQ=255)
                           ttCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTg                                                                                                             >  1:938103/1‑150 (MQ=255)
                            tCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTgg                                                                                                            <  2:625210/150‑1 (MQ=255)
                                     aTACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCaa                                                                                                   >  2:143346/1‑150 (MQ=255)
                                      tACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAAt                                                                                                  <  1:1294527/150‑1 (MQ=255)
                                        ccTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGc                                                                                                >  1:812717/1‑150 (MQ=255)
                                         cTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGcc                                                                                               >  1:278514/1‑150 (MQ=255)
                                                 aGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCATTGCCAcgcccag                                                                                       <  1:868932/150‑7 (MQ=255)
                                                 aGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTtact                                                                                       >  1:132719/1‑150 (MQ=255)
                                                   gATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTtactac                                                                                     <  2:132719/150‑1 (MQ=255)
                                                    aTTTCCGGCTCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCaa                                                                                                   >  2:1090622/1‑135 (MQ=255)
                                                    aTTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCaa                                                                                                   <  1:1090622/135‑1 (MQ=255)
                                                            cgcgTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTcc                                                                            <  1:682440/150‑1 (MQ=255)
                                                              cgTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCaa                                                                          <  2:938103/150‑1 (MQ=255)
                                                              cgTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCaa                                                                          <  1:419730/150‑1 (MQ=255)
                                                               gTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAg                                                                         <  2:278514/150‑1 (MQ=255)
                                                               gTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAg                                                                         >  2:1131633/1‑150 (MQ=255)
                                                               gTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAg                                                                         >  2:452452/1‑150 (MQ=255)
                                                               gTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAg                                                                         >  2:723586/1‑150 (MQ=255)
                                                                 aaTGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAggc                                                                       >  2:1199708/1‑150 (MQ=255)
                                                                    ggAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAggcggc                                                                    >  1:918537/1‑150 (MQ=255)
                                                                     gAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGgcggcg                                                                   <  1:143346/150‑1 (MQ=255)
                                                                         tCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGc                                                               <  2:812717/150‑1 (MQ=255)
                                                                           aaTTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTg                                                             <  1:452452/150‑1 (MQ=255)
                                                                            aTTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACCTTACTACCCCGCCTCCAAGGCGGCGGTGCTgg                                                            >  1:1311300/1‑150 (MQ=255)
                                                                                         gTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGc                                               >  1:706173/1‑150 (MQ=255)
                                                                                          ttttACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACAtctc                                                                                                                                >  1:403088/1‑68 (MQ=255)
                                                                                          ttttACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACAtctc                                                                                                                                <  2:403088/68‑1 (MQ=255)
                                                                                                   aTTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCCGCt                                     >  2:785360/1‑150 (MQ=255)
                                                                                                       gagaGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTAAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGcggcgg                                                                  <  1:1051952/117‑1 (MQ=255)
                                                                                                       gagaGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTAAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGcggcgg                                                                  <  1:1051704/117‑1 (MQ=255)
                                                                                                       gagaGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTAAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGcggcgg                                                                  >  2:1051704/1‑117 (MQ=255)
                                                                                                       gagaGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTAAACATCTCCAGCGTGGAAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGcggcgg                                                                  >  2:1051952/1‑117 (MQ=255)
                                                                                                        agagGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGtt                                <  1:931013/150‑1 (MQ=255)
                                                                                                          agGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAgcgc                                                                                                                  >  1:1117576/1‑66 (MQ=255)
                                                                                                          agGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAgcgc                                                                                                                  <  2:1117576/66‑1 (MQ=255)
                                                                                                          agGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTgg                              >  2:453128/1‑150 (MQ=255)
                                                                                                               gtTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCt                         >  2:333079/1‑150 (MQ=255)
                                                                                                               gtTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCt                         >  2:663211/1‑150 (MQ=255)
                                                                                                                     cTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCTTacaac                   >  2:153019/1‑150 (MQ=255)
                                                                                                                         tGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCTTACAACATTc               <  1:1199708/150‑1 (MQ=255)
                                                                                                                           cgcCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCTTACAACATTCgt             <  2:706173/150‑1 (MQ=255)
                                                                                                                            gcCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTgg                                                            >  1:1404906/1‑102 (MQ=255)
                                                                                                                            gcCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTgg                                                            <  2:1404906/102‑1 (MQ=255)
                                                                                                                               gAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCTTACAACATTCGTGTGa         >  2:690898/1‑150 (MQ=255)
                                                                                                                                      tGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCTTACAACATTCGTGTGAATGCCAt  <  1:153019/150‑1 (MQ=255)
                                                                                                                                                   gTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTAccc                                                                                   <  1:1407097/56‑1 (MQ=255)
                                                                                                                                                   gTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTAccc                                                                                   >  2:1407097/1‑56 (MQ=255)
                                                                                                                                                     |                                                                                                                                      
AAGAGCGGGCGGGCGCCCACTTTCGACTTCCTGGCCGATACCTCGGATCAGGATTTCCGGCGCGTAATGGAAGTCAATTTGTTCAGCCAGTTTTACTGCATTCGAGAGGGTGTTCCGCTGATGCGCCGAGCGGGTGGCGGCAGCATCGTCAACATCTCCAGCGTGGCAGCGCTCCTGGGCGTGGCAATGCCACTTTACTACCCCGCCTCCAAGGCGGCGGTGCTGGGCCTCACCCGTGCAGCGGCAGCTGAGTTGGCACCTTACAACATTCGTGTGAATGCCAT  >  NC_002947/2210129‑2210412

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