Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,471,656 A→T 7.2% intergenic (‑467/‑427) PP_3961 ← / → PP_3962 transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,471,6560AT 7.2% 157.6 / 3.3 56intergenic (‑467/‑427)PP_3961/PP_3962transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base T (2/2);  total (30/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 = 4.11e-01

GACTTGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACA  >  NC_002947/4471514‑4471803
                                                                                                                                              |                                                                                                                                                   
gaCTTGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCa                                                                                                                                               <  7:241837/149‑1 (MQ=255)
   ttGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAggg                                                                                                                                            <  7:259361/149‑1 (MQ=255)
        gCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTc                                                                                                                                       >  2:90883/1‑149 (MQ=255)
        gCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTc                                                                                                                                       <  2:300091/149‑1 (MQ=255)
              ggcgAAGGTTGAGGATTTCGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGTAATGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATg                                                                                                                                 <  1:273151/148‑1 (MQ=255)
                     gTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTgg                                                                                                                          >  1:154262/1‑149 (MQ=255)
                     gTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTgg                                                                                                                          <  2:154262/149‑1 (MQ=255)
                      ttGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTgg                                                                                                                          <  6:117206/148‑1 (MQ=255)
                           gATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCAtt                                                                                                                    <  3:177/149‑1 (MQ=255)
                                 gAATCGCTATTGCCGTGGAAGTGAAAGACGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCg                                                                                                              >  8:257980/1‑149 (MQ=255)
                                   aTCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGgc                                                                                                            >  5:79864/1‑149 (MQ=255)
                                    tCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCAt                                                                                                                                  >  2:83559/1‑126 (MQ=255)
                                    tCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCAt                                                                                                                                  >  7:141086/1‑126 (MQ=255)
                                    tCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCAt                                                                                                                                  <  1:83559/126‑1 (MQ=255)
                                    tCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCAt                                                                                                                                  <  8:141086/126‑1 (MQ=255)
                                              cGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCt                                                                                                  >  3:171335/1‑148 (MQ=255)
                                                   aaGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCCGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAt                                                                                            <  7:250622/149‑1 (MQ=255)
                                                    aGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAtt                                                                                           >  3:250442/1‑149 (MQ=255)
                                                          aGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTa                                                                                                <  3:217290/138‑1 (MQ=255)
                                                          aGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTa                                                                                                >  4:217290/1‑138 (MQ=255)
                                                           gTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACAcc                                                                                    <  8:250898/149‑1 (MQ=255)
                                                               gggAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACAcc                                                                                    >  6:115286/1‑145 (MQ=255)
                                                               gggAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACAcc                                                                                    >  7:150076/1‑145 (MQ=255)
                                                               gggAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGa                                                                                <  6:79864/149‑1 (MQ=255)
                                                                    gCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAATCATTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACAGCAAGATTCaa                                                                           <  3:181275/149‑1 (MQ=255)
                                                                       aCGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGa                                                                        >  5:196785/1‑149 (MQ=255)
                                                                       aCGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTACGCTTACCATTTGACACCAAGATTCAATGa                                                                        >  3:202984/1‑149 (MQ=255)
                                                                            aaaCCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTg                                                                    >  5:74795/1‑148 (MQ=255)
                                                                                  aGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTc                                                             <  5:115286/149‑1 (MQ=255)
                                                                                    cccTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTAc                                                                                               >  8:40322/1‑113 (MQ=255)
                                                                                    cccTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTAc                                                                                               <  7:40322/113‑1 (MQ=255)
                                                                                    cccTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCg                                                            >  7:9172/1‑148 (MQ=255)
                                                                                     ccTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGc                                                           >  8:65475/1‑148 (MQ=255)
                                                                                                ggCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCg                                                   >  1:53008/1‑145 (MQ=255)
                                                                                                       aGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCAcgc                                         <  4:202984/148‑1 (MQ=255)
                                                                                                                 gATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTg                               <  7:257980/148‑1 (MQ=255)
                                                                                                                      gTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTg                         >  4:296091/1‑149 (MQ=255)
                                                                                                                      gTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTg                         >  3:180308/1‑149 (MQ=255)
                                                                                                                       ttAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGc                        >  6:125453/1‑149 (MQ=255)
                                                                                                                       ttAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTACGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATAAGTCGCAGATCACGCGCATGGGCTTGACCTGc                        >  5:7169/1‑149 (MQ=255)
                                                                                                                         aaGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCg                       <  7:65475/148‑1 (MQ=255)
                                                                                                                          aGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGTCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGtt                     <  6:196785/149‑1 (MQ=255)
                                                                                                                            aaaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTgg                   >  5:46429/1‑149 (MQ=255)
                                                                                                                            aaaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTgg                   <  4:171335/149‑1 (MQ=255)
                                                                                                                            aaaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTgg                   <  3:244492/149‑1 (MQ=255)
                                                                                                                            aaaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTgg                   >  1:32604/1‑149 (MQ=255)
                                                                                                                                 cAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAAc              <  8:9172/149‑1 (MQ=255)
                                                                                                                                 aggaaaGCCCTGATTATTCAAGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGcc                            >  7:257713/4‑135 (MQ=255)
                                                                                                                                 aggaaaGCCCTGATTATTCAAGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGcc                            <  8:257713/132‑1 (MQ=255)
                                                                                                                                  aaaaaGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCg                 >  4:187498/1‑145 (MQ=255)
                                                                                                                                  aaaaaGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCg                 <  3:187498/145‑1 (MQ=255)
                                                                                                                                      aGCCCTGATTATTCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCat                                                         <  7:276834/101‑1 (MQ=255)
                                                                                                                                      aGCCCTGATTATTCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCat                                                         >  8:276834/1‑101 (MQ=255)
                                                                                                                                          ttGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGAt      >  5:292653/1‑148 (MQ=255)
                                                                                                                                           tGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCAAGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGa    >  5:5678/1‑149 (MQ=255)
                                                                                                                                            gAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGa    >  4:96470/1‑148 (MQ=255)
                                                                                                                                             aataatCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACa  >  2:26742/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GACTTGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGATGACA  >  NC_002947/4471514‑4471803

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