Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,998,988 G→A 7.9% intergenic (+97/‑208) PP_4406 → / → PP_4407 hypothetical protein/gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,998,9880GA 7.9% 140.0 / 2.0 50intergenic (+97/‑208)PP_4406/PP_4407hypothetical protein/gamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base A (2/2);  total (24/26)
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

CCAAGGTGGACGGCGGCCTGGATGAAGAGCAACAAGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTATTTGCCACGCGTGGACGCAGGC  >  NC_002947/4998839‑4999131
                                                                                                                                                     |                                                                                                                                               
ccAAGGTGGACGGCGGCCTGGATGAAGAGCAACAAGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATgg                                                                                                                                                 <  1:301837/150‑1 (MQ=255)
      tGGACGGCGGCCTGGATGAAGAGCAACAAGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGt                                                                                                                                           <  1:137932/150‑1 (MQ=255)
        gACGGCGGCCTGGATGAAGAGCAACAAGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCa                                                                                                                                         <  2:116328/150‑1 (MQ=255)
           ggcggcCTGAATGAAGAGAAACAAGCAGCAAGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCAAACTTGCCGGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACt                                                                                                                                      <  2:200788/150‑1 (MQ=255)
             cggcCTGGATGAAGAGCAACAAGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGc                                                                                                                                    >  1:613204/1‑150 (MQ=255)
             cggcCTGGATGAAGAGCAACAAGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGc                                                                                                                                    <  2:231224/150‑1 (MQ=255)
                               acaaGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATc                                                                                                                  >  1:1101424/1‑150 (MQ=255)
                                            ccGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGa                                                                                                     >  2:316336/1‑150 (MQ=255)
                                            ccGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGa                                                                                                     >  1:84469/1‑150 (MQ=255)
                                             cGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGACGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGaa                                                                                                    <  1:1302736/150‑1 (MQ=255)
                                                  tGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGACTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATc                                                                                               >  2:1277564/1‑150 (MQ=255)
                                                  tGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGACTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATc                                                                                               >  2:1277511/1‑150 (MQ=255)
                                                       gtgtGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAg                                                                                          >  1:230872/1‑150 (MQ=255)
                                                            gCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGAtgt                                                                                     >  2:1121252/1‑150 (MQ=255)
                                                               gCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGAtgttgt                                                                                  <  2:367089/150‑1 (MQ=255)
                                                               gCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTAAGGAAAAATCGGCAGGAtgttgt                                                                                  <  2:366845/150‑1 (MQ=255)
                                                                       tAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTAcc                                                                                                                           <  2:595696/101‑1 (MQ=255)
                                                                       tAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTAcc                                                                                                                           >  1:595696/1‑101 (MQ=255)
                                                                          gAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTc                                                                       <  2:1101424/150‑1 (MQ=255)
                                                                           aaTAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCg                                                                      >  2:675609/1‑150 (MQ=255)
                                                                                 gccgcGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAg                                                                <  2:230872/150‑1 (MQ=255)
                                                                                       gTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCtt                                                          <  2:84469/150‑1 (MQ=255)
                                                                                              cAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTa                                                   <  1:915348/150‑1 (MQ=255)
                                                                                              cAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTa                                                   <  2:613204/150‑1 (MQ=255)
                                                                                               aaGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCa                                                                                           >  1:118536/1‑109 (MQ=255)
                                                                                               aaGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCa                                                                                           <  2:118536/109‑1 (MQ=255)
                                                                                                  gcgcCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGt                                               >  1:979942/1‑150 (MQ=255)
                                                                                                      cccAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCgg                                                                                                                         >  1:1305723/1‑72 (MQ=255)
                                                                                                      cccAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCgg                                                                                                                         <  2:1305723/72‑1 (MQ=255)
                                                                                                       ccAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTc                                              >  2:1192179/1‑146 (MQ=255)
                                                                                                       ccAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTc                                              <  1:1192179/146‑1 (MQ=255)
                                                                                                       ccAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAgg                                          >  1:464845/1‑150 (MQ=255)
                                                                                                       ccAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAgg                                          >  1:210667/1‑150 (MQ=255)
                                                                                                        cAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAggg                                         >  1:236118/1‑150 (MQ=255)
                                                                                                             tGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGACTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCa                                    <  1:1277511/150‑1 (MQ=255)
                                                                                                             tGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGACTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCa                                    <  1:1277564/150‑1 (MQ=255)
                                                                                                              gATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAg                                   <  1:1121252/150‑1 (MQ=255)
                                                                                                                    cgtcCGGTGGAGTCTACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCACGCATCATCGGCAGGATGTTCTGCACCGAGTACGGTTGAGGGCGTTGGGTTTATTGTCCAGGGTGGCAGGCTGGt                             <  1:248859/150‑1 (MQ=21)
                                                                                                                      tcCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAGCGGCTATAGGATGTTGTGCAGAGAATGCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAg                           >  1:1095990/1‑150 (MQ=38)
                                                                                                                      tcCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAg                                                                >  2:608611/1‑113 (MQ=255)
                                                                                                                      tcCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAg                                                                <  1:608611/113‑1 (MQ=255)
                                                                                                                              cGTCGACGTTCAGGTTGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGTCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTATTTGc                   <  2:574281/150‑1 (MQ=255)
                                                                                                                               gTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCt                                <  2:1151428/136‑1 (MQ=255)
                                                                                                                               gTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCt                                >  1:1151428/1‑136 (MQ=255)
                                                                                                                                      ttCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTATTTGCCACGCGTg           >  2:254288/1‑150 (MQ=255)
                                                                                                                                          ggtggtGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTAttt                     >  1:1208919/1‑136 (MQ=255)
                                                                                                                                          ggtggtGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTAttt                     <  2:1208919/136‑1 (MQ=255)
                                                                                                                                          ggtggtGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTATTTGCCACGCGTGGACg       >  2:423139/1‑150 (MQ=255)
                                                                                                                                               tgCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTATTTGCCACGCGTGGACGCAGgc  <  1:675609/150‑1 (MQ=255)
                                                                                                                                                   tGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTa                                                   <  2:632283/97‑1 (MQ=255)
                                                                                                                                                   tGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTa                                                   >  1:632283/1‑97 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CCAAGGTGGACGGCGGCCTGGATGAAGAGCAACAAGCAGCACGCCCGGTATGAGGGTGTGGCAGCGGCCGGTAAGAATAGCGCCGCGGTTCGTCCAAGGCGCCCCAGCCTGATCGTCGTCCGGTGGCGTCGACGTTCAGGTGGTGCATGGCTTCGTCAACTGCAACGCTACCGGCGCCATCACCGGTCTCAGGAAAAATCGGCAGGATGTTGTGCAGAGAATTCGTTTGAGGTCCTTGGCCTTATTGTCCAGGGTGGCAGGCTGGTAGCTATTTGCCACGCGTGGACGCAGGC  >  NC_002947/4998839‑4999131

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: