Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,956 C→A 7.8% R85L (CGC→CTC) ‡ PP_5039 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,9560CA 7.8% 162.0 / 3.4 51R85L (CGC→CTC) ‡PP_5039hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base A (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 = 8.08e-01

ACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGCGTCGGC  >  NC_002947/5741813‑5742101
                                                                                                                                               |                                                                                                                                                 
acaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcggcg                                                                                                                                             >  1:1134985/1‑150 (MQ=255)
 caGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGGCGACCGGCgt                                                                                                                                            <  2:1088953/150‑1 (MQ=255)
     ttCTGCCAGTCCAGCACTTCGGCATGCTCCAGTCCCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTcg                                                                                                                                        >  2:674568/1‑150 (MQ=255)
        tGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                                                     <  1:655084/150‑1 (MQ=255)
            ggTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGa                                                                                                                                 <  2:514635/149‑1 (MQ=255)
                cAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCAt                                                                                                                             >  2:1005921/1‑150 (MQ=255)
                  gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                                           <  2:1104483/150‑1 (MQ=255)
                    aCTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGt                                                                                                                         <  1:876539/150‑1 (MQ=255)
                                 ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTc                                                                                                              <  1:562442/148‑1 (MQ=255)
                                 ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTc                                                                                                              >  2:562442/1‑148 (MQ=255)
                                 ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                            <  2:1145924/150‑1 (MQ=255)
                                    gTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCa                                                                                                         >  1:817834/1‑150 (MQ=255)
                                     tACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaa                                                                                                        >  1:641861/1‑150 (MQ=255)
                                       ccAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaa                                                                                                      >  2:864077/1‑150 (MQ=255)
                                       ccAGCAGCACGAACACGCGCTGGATAGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGACACGCCGCGCGCCTTCGAGGGCGACCGGCgtgcccag                                                                                                                                      >  2:1192612/1‑113 (MQ=21)
                                       ccAGCAGCACGAACACGCGCTGGATAGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGACACGCCGCGCGCCTTCGAGGGCGACCGGCgtgcccag                                                                                                                                      <  1:1192612/118‑6 (MQ=21)
                                        cagcagCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAAt                                                                                                     >  1:291231/1‑150 (MQ=255)
                                          gcagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATca                                                                                                   <  1:864077/150‑1 (MQ=255)
                                            agcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagc                                                                                                 <  2:291231/150‑1 (MQ=255)
                                                     acGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCccagc                                                                                        >  2:449433/1‑150 (MQ=255)
                                                             gATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                            >  1:1044298/1‑122 (MQ=255)
                                                             gATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                            <  2:1044298/122‑1 (MQ=255)
                                                                cGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgg                                                                             >  2:637037/1‑150 (MQ=255)
                                                                 ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                            <  2:641861/150‑1 (MQ=255)
                                                                        gcTGTTAGTGCGAGCCCTTCTGCAGGCCCTGCATGGCGCCCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                                     <  2:498749/150‑1 (MQ=255)
                                                                            gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCa                                                                 >  1:924969/1‑150 (MQ=255)
                                                                                 cccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCtt                                                            <  2:1112120/150‑1 (MQ=255)
                                                                                  ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                                           <  1:637037/150‑1 (MQ=255)
                                                                                    aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagc                                                                                                 >  2:431263/1‑110 (MQ=255)
                                                                                    aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagc                                                                                                 <  1:431263/110‑1 (MQ=255)
                                                                                              cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGc                                               >  1:998223/1‑150 (MQ=255)
                                                                                                ggCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCa                                             <  2:817834/150‑1 (MQ=255)
                                                                                                ggCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGGAGATCATGCGTGTCAGGTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCa                                             >  1:549713/1‑150 (MQ=255)
                                                                                                     tGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCg                                        >  1:234286/1‑150 (MQ=255)
                                                                                                       gagggCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCg                                                                                                                          >  2:415559/4‑66 (MQ=255)
                                                                                                       gagggCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCg                                                                                                                          <  1:415559/63‑1 (MQ=255)
                                                                                                          ggCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTc                                   <  2:998223/150‑1 (MQ=255)
                                                                                                             gACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGAc                                >  2:294800/1‑150 (MQ=255)
                                                                                                               ccaccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTgg                                                                                                                       <  2:10015/61‑1 (MQ=255)
                                                                                                               ccaccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTgg                                                                                                                       >  1:10015/1‑61 (MQ=255)
                                                                                                                caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAg                             >  1:660642/1‑150 (MQ=255)
                                                                                                                caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAg                             <  2:660642/150‑1 (MQ=255)
                                                                                                                   cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGccc                                                     <  2:1152151/123‑1 (MQ=255)
                                                                                                                   cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGAGGGGTTTGCTGCCACTCTTCGAGccc                                                     >  1:1152151/1‑123 (MQ=255)
                                                                                                                      ggCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGCCCAGgccgcc                       <  2:203961/150‑1 (MQ=255)
                                                                                                                        cACGCCGCGCGCCTTCGAGGGCGACCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGGTCCCGCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCaa                     >  2:96949/1‑150 (MQ=255)
                                                                                                                         aCGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTc                                   >  1:413725/1‑135 (MQ=255)
                                                                                                                         aCGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTc                                   <  2:413725/135‑1 (MQ=255)
                                                                                                                             cGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAgc                >  2:393055/1‑150 (MQ=255)
                                                                                                                                   cccTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATgg          <  2:924969/150‑1 (MQ=255)
                                                                                                                                           ggCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGCGTCGGc  >  2:688733/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
ACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGGCGTCGGC  >  NC_002947/5741813‑5742101

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