Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,230,075 G→C 10.7% A63G (GCG→GGG)  PP_4607 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,230,0750GC10.7% 131.5 / 5.7 44A63G (GCG→GGG) PP_4607transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base C (2/3);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATT  >  NC_002947/5229929‑5230213
                                                                                                                                                  |                                                                                                                                          
gCCAACCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAgccgcgc                                                                                                                                           <  1:446058/148‑1 (MQ=255)
gCCAACCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAgccgcgc                                                                                                                                           >  2:446058/1‑148 (MQ=255)
                  ccACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCt                                                                                                                       >  2:16599/1‑150 (MQ=255)
                  aaacatccAGCCCCAGGCAGGCTGCCAACGCCAGGCCTGCACCCCGACCCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGAGCCGCGGGAGCGGCGCGTGTCt                                                                                                                       >  2:334924/7‑150 (MQ=255)
                    aCGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGc                                                                                                                     <  2:383088/150‑1 (MQ=255)
                    aCGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGc                                                                                                                     >  2:427081/1‑150 (MQ=255)
                        ccAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCaaa                                                                                                                 <  1:144081/150‑1 (MQ=255)
                        ccAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCaaa                                                                                                                 >  2:416669/1‑150 (MQ=255)
                        caaGCACCAGTCAGGCTGCCAACGCCAGGCCTGCCCCCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCAGCGATGTAGCGCTGCAGGGCGCCCTGCTCTTCTTCCGCCAGCCGCCCCGCCGGGGCGGCGCTTTTCTGCCaaa                                                                                                                 <  2:279188/148‑1 (MQ=255)
                                      gcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTgg                                                                                                   <  2:513057/150‑1 (MQ=255)
                                            aaCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAg                                                                                             >  1:222616/1‑150 (MQ=255)
                                            aaCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAg                                                                                             >  1:729784/1‑150 (MQ=255)
                                                ccATGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCCCCGgcggcgcggctgg                                                                                                                            <  2:309706/115‑4 (MQ=255)
                                                ccAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTg                                                                                         >  2:414998/1‑150 (MQ=255)
                                                ccAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCCCCGgcggcgcggctgg                                                                                                                            >  1:309972/1‑112 (MQ=255)
                                                ccAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCCCCGgcggcgcggctgg                                                                                                                            >  1:309706/1‑112 (MQ=255)
                                                ccAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCCCCGgcggcgcggctgg                                                                                                                            <  2:309972/115‑4 (MQ=255)
                                                  aGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTgcg                                                                                       >  1:428094/1‑150 (MQ=255)
                                                        gCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgtgt                                                                                 <  1:416669/150‑1 (MQ=255)
                                                        gCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgtgt                                                                                 <  2:428094/150‑1 (MQ=255)
                                                           aCCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCg                                                                              <  1:455676/150‑1 (MQ=255)
                                                           aCCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAAGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCAGCCGGGGCCGAGCTTTTCTGCCACAGCACATGCGCATGGTCGCAGGCCGCGCggcggccg                                                                              >  2:82623/1‑144 (MQ=255)
                                                             ccgccgCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCGTCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCggg                                                                            <  1:16599/150‑1 (MQ=255)
                                                             ccgccgCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCGTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCggg                                                                            >  2:573248/1‑150 (MQ=255)
                                                              cgccgcCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCgggg                                                                           <  2:212429/150‑1 (MQ=255)
                                                               gccgccAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGt                                                                          >  1:406774/1‑150 (MQ=255)
                                                                   ccagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGc                                                                      <  1:573248/150‑1 (MQ=255)
                                                                             cgcgCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAg                                                            >  1:720964/1‑150 (MQ=255)
                                                                                 ccGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCt                                                        <  2:406774/150‑1 (MQ=255)
                                                                                      ggtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAAt                                                   <  1:427081/150‑1 (MQ=255)
                                                                                       gtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATg                                                  >  1:64204/1‑150 (MQ=255)
                                                                                         cTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgcc                                                <  2:64204/150‑1 (MQ=255)
                                                                                          tGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccg                                               >  1:827679/1‑150 (MQ=255)
                                                                                          tGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGATGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccg                                               >  1:42464/1‑150 (MQ=255)
                                                                                           ggCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccac                                              <  2:827999/150‑3 (MQ=255)
                                                                                            gCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccgcc                                             <  2:118756/150‑1 (MQ=255)
                                                                                              cATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGcccc                                           >  2:542487/1‑150 (MQ=255)
                                                                                                    aTCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGc                                     >  1:237835/1‑150 (MQ=255)
                                                                                                               gcTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCggg                                                                            <  2:341863/100‑1 (MQ=255)
                                                                                                               gcTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCggg                                                                            >  1:341863/1‑100 (MQ=255)
                                                                                                               gcTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGtt                          <  1:542487/150‑1 (MQ=255)
                                                                                                                          tCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATc               <  2:720964/150‑1 (MQ=255)
                                                                                                                               ctctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGtt          <  2:137236/150‑1 (MQ=255)
                                                                                                                                  ttCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGc                                                                      <  1:656200/87‑1 (MQ=255)
                                                                                                                                  ttCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGc                                                                      >  2:656200/1‑87 (MQ=255)
                                                                                                                                       ccGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCAtt  <  2:222616/150‑1 (MQ=255)
                                                                                                                                               agcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCg             <  1:671886/130‑1 (MQ=255)
                                                                                                                                               agcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCg             >  2:671886/2‑131 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GCCAACCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATT  >  NC_002947/5229929‑5230213

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