Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 3,254,711 C→A G232G (GGG→GGT)  PP_2851 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,254,7110CA98.0% 165.3 / ‑3.1 49G232G (GGG→GGT) PP_2851hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (21/27);  total (21/28)
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.33e-01

GCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGCCCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAACCCGCGGGCGTAGGTTGAGTAGGT  >  NC_002947/3254588‑3254839
                                                                                                                           |                                                                                                                                
gcgGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCgg                                                                                                                      >  2:681976/1‑136 (MQ=255)
gcgGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCgg                                                                                                                      >  1:435086/1‑136 (MQ=255)
  gggTAAACCCGCTCCTACGGGGCGACCGCGCATGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTc                                                                                                                    <  1:1362012/136‑1 (MQ=255)
           cGCTCCTACGGGGCTACCGCGCTGGCCTTGTAATTTTGCGCAGTACCGTTG‑CCCCCAGAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATTGCGGTCTTGGTGCCGc                                                                                                          <  1:618489/136‑1 (MQ=255)
            gCTCCTACGGGGTGACCGCGGAGGCCGGCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGGAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGc                                                                                                          <  1:537731/136‑1 (MQ=255)
            gCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGc                                                                                                          <  1:632573/136‑1 (MQ=255)
            gCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGc                                                                                                          <  1:904660/136‑1 (MQ=255)
                 tACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCa                                                                                                     >  1:975713/1‑136 (MQ=255)
                 tACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGACGCACCCa                                                                                                     >  1:405817/1‑136 (MQ=255)
                   cGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGt                                                                                                   <  2:838336/136‑1 (MQ=255)
                    ggggCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGCCCTTCATAGCGGTCATGGTGCCGCACCCAGTc                                                                                                  <  1:1204043/136‑1 (MQ=255)
                         gACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATgg                                                                                             <  2:926818/136‑1 (MQ=255)
                         gACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATgg                                                                                             <  2:144601/136‑1 (MQ=255)
                         gACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATgg                                                                                             <  2:85468/136‑1 (MQ=255)
                             gcgcAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTccc                                                                                         >  1:99204/1‑136 (MQ=255)
                             gcgcAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTccc                                                                                         >  1:155032/1‑136 (MQ=255)
                                    ccGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCg                                                                                  <  1:681976/136‑1 (MQ=255)
                                      gTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGtt                                                                                <  2:1268822/136‑1 (MQ=255)
                                         gAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACg                                                                             <  2:264330/136‑1 (MQ=255)
                                                   cAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGgc                                                                   <  2:50755/136‑1 (MQ=255)
                                                    aGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGgcg                                                                  <  2:1087258/136‑1 (MQ=255)
                                                    aGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGgcc                                                                  >  1:1309732/1‑135 (MQ=255)
                                                         aGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATg                                                             >  1:1306529/1‑136 (MQ=255)
                                                             gCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGa                                                         <  1:262134/136‑1 (MQ=255)
                                                             gCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGa                                                         <  2:297842/136‑1 (MQ=255)
                                                              ccccccACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAg                                                        >  1:331221/1‑136 (MQ=255)
                                                              ccccccACAGGTGTGCGTGTTGCCAGCATCATCTGAGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAg                                                        <  2:948477/136‑1 (MQ=255)
                                                                   cacaGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGt                                                   <  2:975713/136‑1 (MQ=255)
                                                                   cacaGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGACATGGTGCCGCACCCAGTCCATGGTCCCTCCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGatgg                                                   >  1:16267/1‑133 (MQ=255)
                                                                      aGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGt                                                >  1:583130/1‑136 (MQ=255)
                                                                         tgtgCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAgg                                             <  1:655352/136‑1 (MQ=255)
                                                                           tgCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTCGCGCAATGTCGAGGAAGTTGTAGgtg                                           >  2:510870/1‑136 (MQ=255)
                                                                              gtgtTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTg                                        >  1:767255/1‑136 (MQ=255)
                                                                              gtgtTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTg                                        >  2:1334739/1‑136 (MQ=255)
                                                                                gtTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCa                                                                                                     >  1:22976/1‑73 (MQ=255)
                                                                                gtTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCa                                                                                                     <  2:22976/73‑1 (MQ=255)
                                                                                gtTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGAc                                      >  1:183550/1‑136 (MQ=255)
                                                                                gtTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGAc                                      >  1:1332228/1‑136 (MQ=255)
                                                                                gtTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGAc                                      >  2:892506/1‑136 (MQ=255)
                                                                                             tcTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGTGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAAc                         <  1:434280/136‑1 (MQ=255)
                                                                                             tcTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAAc                         <  2:99204/136‑1 (MQ=255)
                                                                                             tcTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAAc                         <  2:155032/136‑1 (MQ=255)
                                                                                             tcTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAAc                         <  1:767393/136‑1 (MQ=255)
                                                                                             tcTGTGACAGCAGTGAGGCTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAAc                         <  1:653841/136‑1 (MQ=255)
                                                                                                           gAGGCTTGTCCGGTTGACCTTCATAGCGGTCATCGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTAGTAGGTGTTGACCAGAATGTCCAACCCGCGGGCGTAGGt           >  1:95131/1‑136 (MQ=255)
                                                                                                               cTTGTCCGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAACCCGCGGGCGTAGGTTGAg       <  1:861940/136‑1 (MQ=255)
                                                                                                                    ccGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAACCCGCGGGCGTAGGt           >  2:776701/1‑127 (MQ=255)
                                                                                                                    ccGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAACCCGCGGGCGTAGGt           <  1:776701/127‑1 (MQ=255)
                                                                                                                    ccGGTTGACCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAACCCGCGGGCGTAGGTTGAGTAGgt  >  2:765624/1‑136 (MQ=255)
                                                                                                                           |                                                                                                                                
GCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGCGTGTTGCCAGCATCATCTGTGACAGCAGTGAGGCTTGTCCGGTTGCCCTTCATAGCGGTCATGGTGCCGCACCCAGTCCATGGTCCCACCTTCGTTACGCCCCAGTGGCGCAATGTCGAGGAAGTTGTAGGTGTTGACCAGAATGTCCAACCCGCGGGCGTAGGTTGAGTAGGT  >  NC_002947/3254588‑3254839

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