Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,269 G→C 15.1% intergenic (+45/+97) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,2690GC15.1% 138.8 / 15.2 52intergenic (+45/+97)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (19/25);  new base C (6/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TGCTGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTT  >  NC_002947/1707121‑1707410
                                                                                                                                                    |                                                                                                                                             
tgctgTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCcgcgc                                                                                                                                              >  1:572468/1‑150 (MQ=255)
    gTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgggc                                                                                                                                          >  2:104878/1‑150 (MQ=255)
     ttCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCAATATTGGGGCCTGCGAGATCGGGCGCCGCgcgggcg                                                                                                                                         >  2:977125/1‑150 (MQ=255)
           acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                   >  1:149787/1‑150 (MQ=255)
           acCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                   >  2:737085/1‑150 (MQ=255)
             cgaggagAAAGGCCTGGACGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCTCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 <  2:788277/145‑1 (MQ=255)
             ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 <  2:1051150/150‑1 (MQ=255)
             ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                 >  1:647231/1‑150 (MQ=255)
              cagcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                >  2:89691/1‑150 (MQ=255)
              cagcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                                <  1:737085/150‑1 (MQ=255)
                  agAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                            >  1:307122/1‑150 (MQ=255)
                               cGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaac                                                                                                               <  1:444472/150‑1 (MQ=255)
                                  cGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACAtt                                                                                                            <  1:859542/150‑1 (MQ=255)
                                     gCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATc                                                                                                         <  2:591391/150‑1 (MQ=255)
                                      cTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCg                                                                                                        <  1:1212308/150‑1 (MQ=255)
                                        ggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCTGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                      <  1:104878/150‑1 (MQ=255)
                                        ggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGc                                                                                                      <  1:188964/150‑1 (MQ=255)
                                         gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGTCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCg                                                                                                     >  1:992402/1‑150 (MQ=255)
                                         gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                            <  1:1180752/127‑1 (MQ=255)
                                         gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                            >  2:1180752/1‑127 (MQ=255)
                                                 ccGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGACCGAACGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGa                                                                                             >  1:1253595/1‑150 (MQ=255)
                                                  cGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGaa                                                                                            >  1:1039848/1‑150 (MQ=255)
                                                  cGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGaa                                                                                            >  1:307030/1‑150 (MQ=255)
                                                                 cgcTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGAt                                                                             <  2:647231/150‑1 (MQ=255)
                                                                      cATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                        <  2:894282/150‑1 (MQ=255)
                                                                            cGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTg                                                                  <  1:202689/150‑1 (MQ=255)
                                                                             gAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGt                                                                 >  2:319619/1‑150 (MQ=255)
                                                                              aGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTc                                                                <  1:319619/150‑1 (MQ=255)
                                                                              aGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTc                                                                <  1:143771/150‑1 (MQ=255)
                                                                               ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaac                                                                                                               >  1:218950/1‑102 (MQ=255)
                                                                               ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaac                                                                                                               <  2:218950/102‑1 (MQ=255)
                                                                                  aaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCttt                                                            <  2:572468/150‑1 (MQ=255)
                                                                                  aaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCttt                                                            <  2:1039848/150‑1 (MQ=255)
                                                                                   acacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                           >  2:114997/1‑150 (MQ=255)
                                                                                     acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                         >  1:692436/1‑150 (MQ=255)
                                                                                             aGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccag                                                 >  1:491908/1‑150 (MQ=255)
                                                                                                 gCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCtt                                             <  2:692436/150‑1 (MQ=255)
                                                                                                   aGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAg                                           <  2:697060/150‑1 (MQ=255)
                                                                                                              ggCGCAGCATTTTTGGCGCCCGCGAGATGGAACGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                >  1:33106/1‑150 (MQ=255)
                                                                                                              ggCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgggg                                                                      >  2:169115/1‑112 (MQ=255)
                                                                                                              ggCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgggg                                                                      <  1:169115/112‑1 (MQ=255)
                                                                                                                    gcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAc                          >  2:589196/1‑150 (MQ=255)
                                                                                                                         tttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATATAACAATCAACACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGc                     <  2:33106/150‑1 (MQ=255)
                                                                                                                           tgcagaccGCGAGATCGAGCGCCGCGCGAGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgt                   <  1:1231421/144‑1 (MQ=255)
                                                                                                                              cgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGCTGCCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                <  1:977125/150‑1 (MQ=255)
                                                                                                                                   gcgAGATCTAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGGGCCAc           <  2:491908/150‑1 (MQ=255)
                                                                                                                                       gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGc       >  1:1079570/1‑150 (MQ=255)
                                                                                                                                        aTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCt      <  1:89691/150‑1 (MQ=255)
                                                                                                                                        aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCt      >  2:716574/1‑150 (MQ=39)
                                                                                                                                         tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCTGcgg     >  2:579270/1‑148 (MQ=255)
                                                                                                                                         tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTg     >  1:795640/1‑150 (MQ=255)
                                                                                                                                          cGAGCGCCGCCCGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTgg    <  1:716574/150‑1 (MQ=39)
                                                                                                                                           gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt   >  2:624659/1‑150 (MQ=255)
                                                                                                                                           gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt   >  1:440670/1‑150 (MQ=255)
                                                                                                                                            aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt  >  2:624975/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TGCTGTTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTT  >  NC_002947/1707121‑1707410

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: