Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,004,613 G→A 13.6% intergenic (+60/+90) PP_1787 → / ← PP_1788 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,004,6130GA13.6% 143.0 / 14.3 59intergenic (+60/+90)PP_1787/PP_1788hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base A (3/5);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

CGCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGG  >  NC_002947/2004468‑2004749
                                                                                                                                                 |                                                                                                                                        
cgcCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGAttt                                                                                                                                      <  1:200590/150‑1 (MQ=255)
      ttttgtCAAGTTTTCTACGTTCGCGGTCGAAAATTAATAGTGAGTTTTTTGACGTGACTGAAATGACGCTTCATGCTTAATAAGATATTTATATTCTACCTCCCCTTGACATCCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATa                                                                                                                                <  2:752498/145‑1 (MQ=255)
        ttttCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAg                                                                                                                              <  1:215205/150‑1 (MQ=255)
             aaGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGGTCGTTATATTCTACGTCCCCTTGACATGCCCCGAAAATGAAATAGGTTATTGCAGATTTGCAATAAGCTAtt                                                                                                                         >  1:848904/1‑150 (MQ=255)
             aaGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTAtt                                                                                                                         <  1:1498441/150‑1 (MQ=255)
              aGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGAAGTGCGTCAAATGACTCTTAATGCTTAAGAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTAttt                                                                                                                        >  1:691965/1‑150 (MQ=255)
                                  gAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGc                                                                                                    >  1:128796/1‑150 (MQ=255)
                                      tttAAAAGTGAGTTTTTTGACGTGACTGAAATTACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATc                                                                                                <  1:257916/150‑1 (MQ=255)
                                      tttAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATc                                                                                                >  1:169654/1‑150 (MQ=255)
                                          aaaGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATCGCTTATTGCAGATTTGCAATAAGCTATTTCTAAAttt                                                                                                                <  1:1135366/130‑1 (MQ=255)
                                          aaaGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAAttt                                                                                                                >  2:1135366/1‑130 (MQ=255)
                                          aaaGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGtt                                                                                            <  2:743599/150‑1 (MQ=255)
                                              tGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGtttttt                                                                                        <  2:205898/150‑1 (MQ=255)
                                                 gTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGa                                                                                     >  2:1444369/1‑150 (MQ=255)
                                                  ttttttGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAg                                                                                                                              >  2:1240991/1‑108 (MQ=255)
                                                  ttttttGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAATAAg                                                                                                                              <  1:1240991/108‑1 (MQ=255)
                                                    ttttGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGAtgt                                                                                  <  2:169654/150‑1 (MQ=255)
                                                      ttGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGa                                                                                <  1:1444369/150‑1 (MQ=255)
                                                          cGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCt                                                                                                                      <  2:642036/108‑1 (MQ=255)
                                                          cGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCt                                                                                                                      >  1:642036/1‑108 (MQ=255)
                                                           gTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCg                                                                           >  1:817457/1‑150 (MQ=255)
                                                              aCTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGTGGCGATCGGTTTTTTGGATGTGAAACCGCTc                                                                        <  2:1175861/150‑1 (MQ=255)
                                                               cTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTAttt                                                                                                                        <  1:606760/101‑1 (MQ=255)
                                                               cTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTAttt                                                                                                                        >  2:606760/1‑101 (MQ=255)
                                                               cTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAAttt                                                                                                                <  2:1541604/109‑1 (MQ=255)
                                                               cTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAAttt                                                                                                                >  1:1541604/1‑109 (MQ=255)
                                                                     tGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAg                                                                 <  2:128796/150‑1 (MQ=255)
                                                                             aatgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGAt                                                                                    <  1:1623264/123‑1 (MQ=255)
                                                                             aatgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGAt                                                                                    >  2:1623264/1‑123 (MQ=255)
                                                                               tgcttaatAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGAt                                                       <  2:1021187/150‑1 (MQ=255)
                                                                                    aataaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCa                                                                  <  1:1108792/134‑1 (MQ=255)
                                                                                    aataaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCa                                                                  >  2:1108792/1‑134 (MQ=255)
                                                                                      taaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAAt                                                                                                                                 <  2:584567/69‑1 (MQ=38)
                                                                                      taaGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAAATCTGCAAt                                                                                                                                 >  1:584567/1‑69 (MQ=38)
                                                                                         gATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTAt                                             <  2:691965/150‑1 (MQ=255)
                                                                                           tatTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGa                                           <  2:817457/150‑1 (MQ=255)
                                                                                                 taGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAgtgt                                     <  1:118213/150‑1 (MQ=255)
                                                                                                   gTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTa                                                                                                                           >  1:729147/1‑62 (MQ=255)
                                                                                                   gTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTa                                                                                                                           <  2:729147/62‑1 (MQ=255)
                                                                                                    tCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGa                                  >  1:1269341/1‑150 (MQ=255)
                                                                                                         cTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGttttt                                                                                         <  2:1074124/90‑1 (MQ=255)
                                                                                                         cTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGttttt                                                                                         >  1:1074124/1‑90 (MQ=255)
                                                                                                         cTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTaa                             >  2:1273008/1‑150 (MQ=255)
                                                                                                          tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAg                            <  2:1269341/150‑1 (MQ=255)
                                                                                                          tCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAg                            >  1:127549/1‑150 (MQ=255)
                                                                                                              gttGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTg                        >  1:536097/2‑150 (MQ=255)
                                                                                                               ttGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGtttttt                                                                                        >  1:726976/1‑85 (MQ=255)
                                                                                                               ttGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGtttttt                                                                                        <  2:726976/85‑1 (MQ=255)
                                                                                                                    aTGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTtgt                                                          <  2:1610790/110‑1 (MQ=255)
                                                                                                                    aTGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTtgt                                                          >  1:1610790/1‑110 (MQ=255)
                                                                                                                       ccccGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGAtgt                                                                                  >  2:660429/1‑83 (MQ=255)
                                                                                                                       ccccGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGAtgt                                                                                  <  1:660429/83‑1 (MQ=255)
                                                                                                                        cccGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGAt              >  1:975537/1‑150 (MQ=255)
                                                                                                                          cGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTg            <  2:848904/150‑1 (MQ=255)
                                                                                                                           gAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTgg           <  1:938499/150‑1 (MQ=255)
                                                                                                                                        cTTATTGCAAATCTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATgg  <  2:1457788/146‑1 (MQ=255)
                                                                                                                                        cTTATTGCAAATCTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCCCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATgg  >  1:1457788/1‑146 (MQ=255)
                                                                                                                                                aGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTAt                                             <  2:313096/95‑1 (MQ=255)
                                                                                                                                                aGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTAt                                             >  1:313096/1‑95 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
CGCCCGTATTTTCAAGGTTTCAAAGTTCGCGGTCGAAATTTAAAAGTGAGTTTTTTGACGTGACTGAAATGACGCTTAATGCTTAATAAGATATTTATAGTCTACCTCCCCTTGACATGCCCCGAAAATGAAATAGCTTATTGCAGATTTGCAATAAGCTATTTCTAAATTTGTTAGCTGAGGCGATCGGTTTTTTGGATGTGAAACCGCTCGCCTCAGGCTTTGTGATGTGTATTTATGATAGTGTTGATCTAAGCTTGCAATACGGATTGGTTTTAATGG  >  NC_002947/2004468‑2004749

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