Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,644,017 G→A 8.1% L25L (CTG→CTA)  PP_2314 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,644,0170GA 8.1% 139.3 / 3.2 49L25L (CTG→CTA) PP_2314hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (2/2);  total (28/21)
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.86e-01

TAAAGCGCTGGTTCGGTGACGCCTCCCACCGCTGCCCAGCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGACCGTTTCGACCTGGCTGT  >  NC_002947/2643876‑2644163
                                                                                                                                               |                                                                                                                                                  
tAAAGCGCTGGTTCGGTGACGCCTCCCACCGCTGCCCAGCGGGATTCGCGAGGCGGAGAAGCAAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTcgg                                                                                                                                               >  1:290727/1‑149 (MQ=255)
 aaaGCGCTGGTTCGGTGACGCCTCCCACCGCTGCCCAGCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTcggcgg                                                                                                                                            <  3:83822/149‑1 (MQ=255)
               gTGACGCCTCCCACCGCTGCCCAGCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCa                                                                                                                              <  3:357455/149‑1 (MQ=255)
               gTGACGCCTCCCACCGCTGCCCAGCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCa                                                                                                                              <  8:246705/149‑1 (MQ=255)
                      cTCCCACCGCTGCCCAGCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCtttt                                                                                                                       <  2:641/149‑1 (MQ=255)
                           aCCGCTGCCCAGCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGt                                                                                                                   >  3:10742/1‑148 (MQ=255)
                                      gCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGc                                                                                                       <  7:98130/149‑1 (MQ=255)
                                               gcgAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGGTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACgg                                                                                              >  1:267940/1‑149 (MQ=255)
                                                      ggAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAgg                                                                                       <  5:108431/149‑1 (MQ=255)
                                                      ggAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAgg                                                                                       >  3:355522/1‑149 (MQ=255)
                                                       gag‑ag‑aAGCTCCATGAGTTGCCCACGTACCAGTGGCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAggg                                                                                      <  7:292414/149‑1 (MQ=255)
                                                              aaGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAgg                                                                                 >  5:159362/1‑149 (MQ=255)
                                                                    cATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCACGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTg                                                                           >  1:278631/1‑149 (MQ=255)
                                                                    cATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCACGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTg                                                                           >  1:273813/1‑149 (MQ=255)
                                                                         gTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCt                                                                                                                          <  1:34926/97‑1 (MQ=255)
                                                                         gTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCt                                                                                                                          >  2:34926/1‑97 (MQ=255)
                                                                          ttGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTc                                                                     <  6:159362/149‑1 (MQ=255)
                                                                          gtGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTc                                                                     <  5:30079/148‑1 (MQ=255)
                                                                            gCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGaa                                                                                                          <  1:180380/110‑1 (MQ=255)
                                                                            gCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGaa                                                                                                          >  2:180380/1‑110 (MQ=255)
                                                                            gCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCACGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCCAGGGGCAGGCGCGTGCCTTTCt                                                                    >  4:192626/1‑148 (MQ=255)
                                                                              ccACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGaa                                                                                                          <  6:257079/108‑1 (MQ=255)
                                                                              ccACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGaa                                                                                                          >  5:257079/1‑108 (MQ=255)
                                                                              ccACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGAt                                                                 <  6:138262/149‑1 (MQ=255)
                                                                                    aCCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGtt                                                           >  5:194941/1‑149 (MQ=255)
                                                                                      cAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCg                                                         >  2:299172/1‑149 (MQ=255)
                                                                                      cAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCg                                                         >  1:278602/1‑149 (MQ=255)
                                                                                       aGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGc                                                        <  2:267940/149‑1 (MQ=255)
                                                                                        gtgtCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCATGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGcc                                                       >  7:329385/1‑149 (MQ=255)
                                                                                               tttGAGCCCTTTGCCTGCCAGCCAGGCGTCNCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCg                                                >  1:34902/1‑149 (MQ=255)
                                                                                               tttGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAg                                                  <  4:10742/147‑1 (MQ=255)
                                                                                               tttGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCg                                                >  6:31138/1‑149 (MQ=255)
                                                                                                   aGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCt                                             <  2:290727/148‑1 (MQ=255)
                                                                                                       cTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCAc                                        <  6:194941/149‑1 (MQ=255)
                                                                                                           gcctgccAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCAc                                        >  8:27203/1‑145 (MQ=255)
                                                                                                           gcctgccAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCAc                                        <  7:27203/145‑1 (MQ=255)
                                                                                                                ccagccagGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAg                               >  8:131602/1‑149 (MQ=255)
                                                                                                                 cagccagGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAg                               >  4:175757/1‑148 (MQ=255)
                                                                                                                      agGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATgg                          >  1:210254/1‑148 (MQ=255)
                                                                                                                          gTCCCGCACACCGCGCATCCTACTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAAt                      >  8:267376/1‑148 (MQ=255)
                                                                                                                          gTCCCGCAAACCGCGCATCCTACTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAAt                      >  8:262435/1‑148 (MQ=255)
                                                                                                                           tCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGa                    >  5:283189/1‑149 (MQ=255)
                                                                                                                             ccGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGt                                           >  7:237599/1‑124 (MQ=255)
                                                                                                                             ccGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGt                                           <  8:237599/124‑1 (MQ=255)
                                                                                                                             ccGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGAcc                  >  7:146472/1‑149 (MQ=255)
                                                                                                                                  gacCGCGCATCCTACTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGACCGTTTc             <  7:267376/148‑1 (MQ=255)
                                                                                                                                  aacCGCGCATCCTACTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGACCGTTTc             <  7:262435/148‑1 (MQ=255)
                                                                                                                                   acCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGACCGTTTc             >  6:370255/1‑148 (MQ=255)
                                                                                                                                             cTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGACCGTTTCGACCTGGCtgt  >  6:158053/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TAAAGCGCTGGTTCGGTGACGCCTCCCACCGCTGCCCAGCGGGATTCGCGAGGCGGAG‑AG‑AAGCTCCATGAGTTGCCCACGTACCAGTGTCTGTTTGAGCCCTTTGCCTGCCAGCCAGGCGTCCCGCACACCGCGCATCCTGCTCGGCGGCCGGCATCAACCCAGCCTTTTGCGTCATCTTGAAGGCCTGTCACGGCGCAAGGGGCAGGCGCGTGCCTTTCTGATCCAGTTCGCCGATAGCGGCTGTCACCTCGCGCAGTATGGCAATGACCGTTTCGACCTGGCTGT  >  NC_002947/2643876‑2644163

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