Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,905,502 A→G 7.3% intergenic (‑96/+32) PP_4313 ← / ← PP_4314 peptidylprolyl isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,905,5020AG 7.3% 153.0 / 2.6 55intergenic (‑96/+32)PP_4313/PP_4314peptidylprolyl isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/30);  new base G (2/2);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

GCAGCTTTGCCGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCAGGCTTGG  >  NC_002947/4905354‑4905648
                                                                                                                                                       |                                                                                                                                                  
gCAGCTTTGCCGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGa                                                                                                                                                    >  8:318592/1‑149 (MQ=255)
 caaagctGCCGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCAAGCAGAATTGTACGGAGGCATTCGCCGAc                                                                                                                                                   >  6:102608/7‑149 (MQ=255)
          cGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAATAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCt                                                                                                                                          <  4:60184/149‑1 (MQ=255)
            tctcCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTcc                                                                                                                                        <  3:9033/149‑1 (MQ=255)
             ctcCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCg                                                                                                                                       >  1:98201/1‑149 (MQ=255)
                ccTCGACCATCGCGAAGATCTGCAGTTCTTGGCTCTCGCCCCATACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCCGTCGGCGAATGCCTcc                                                                                                                                        <  1:316103/146‑1 (MQ=255)
                ccTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTa                                                                                                                                     >  4:82896/1‑148 (MQ=255)
                ccTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTa                                                                                                                                     >  2:98282/1‑148 (MQ=255)
                ccTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTa                                                                                                                                     >  3:358196/1‑148 (MQ=255)
                 cTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACa                                                                                                                                   <  3:192223/149‑1 (MQ=255)
                    gACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGa                                                                                                                                 <  4:67822/148‑1 (MQ=255)
                    gACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGAc                                                                                                                                <  5:354117/149‑1 (MQ=255)
                      ccAGGTCGATGATCTGCAGTTATTGGCTCACTGGTAAAATCCAAACGCTTAAAAACAAAACCGACCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGAAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACAt                                                                                                                              <  1:238174/149‑1 (MQ=255)
                             gatgatCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTg                                                                                                                       <  4:137042/149‑1 (MQ=255)
                              attatCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAATCCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGa                                                                                                                      <  7:329935/146‑1 (MQ=255)
                              atgatCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTg                                                                                                                       >  3:193397/1‑148 (MQ=255)
                              atgatCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGa                                                                                                                      <  1:98282/149‑1 (MQ=255)
                                      cAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCa                                                                                                                  <  5:102608/145‑1 (MQ=255)
                                       aGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCGCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGtt                                                                                                                        <  3:95097/138‑1 (MQ=255)
                                       aGTTCTTGGCTCACCGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGtt                                                                                                                        >  4:95097/1‑138 (MQ=255)
                                             tGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGc                                                                                                        <  3:121990/148‑1 (MQ=255)
                                              ggCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGc                                                                                                        <  8:322566/147‑1 (MQ=255)
                                                     gggTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACAt                                                                                                                              <  5:1712/118‑1 (MQ=255)
                                                     gggTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACAt                                                                                                                              >  6:1712/1‑118 (MQ=255)
                                                        gTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGATGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACctc                                                                                             >  4:260800/1‑149 (MQ=255)
                                                           cacTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTg                                                                                          <  4:193397/149‑1 (MQ=255)
                                                             cTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGc                                                                                        <  1:206016/149‑1 (MQ=255)
                                                             cTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGc                                                                                        <  5:271642/149‑1 (MQ=255)
                                                             cTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGc                                                                                        >  3:147982/1‑149 (MQ=255)
                                                                cAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTGACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGCCGTTGATTCAGCTGTTCTGCTCGAACACCTCTTGGCCg                                                                                      >  8:192018/1‑148 (MQ=255)
                                                                cAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGcc                                                                                       <  2:98201/147‑1 (MQ=255)
                                                                 aaaCGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGa                                                                                    <  7:318592/149‑1 (MQ=255)
                                                                    cGCTGAAAAACAAACCCGCCATTTTCTCAGGAATGCAGCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAg                                                                                 <  4:229919/149‑1 (MQ=255)
                                                                       tGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCCGTCGGCGAATGCCTcc                                                                                                                                        >  5:282528/1‑91 (MQ=37)
                                                                       tGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCCGTCGGCGAATGCCTcc                                                                                                                                        <  6:282528/91‑1 (MQ=37)
                                                                         aaaaaCAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTg                                                                                          <  8:216004/135‑1 (MQ=255)
                                                                         aaaaaCAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTg                                                                                          >  7:216004/1‑135 (MQ=255)
                                                                          aaaacaaaacCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTc                                                                                                                                         <  6:358592/87‑1 (MQ=255)
                                                                              caaaacCGCCCTTTTCTCAGGAATGCACCGGTTTGCGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTgg                                                                         >  1:56906/1‑148 (MQ=255)
                                                                              caaaacCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGATGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTc                                                                       <  3:260800/149‑1 (MQ=255)
                                                                                 aacCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCtca                                                                    >  2:238757/1‑147 (MQ=255)
                                                                                         ttttCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCg                                                              >  4:335438/1‑147 (MQ=255)
                                                                                           ttCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAg                                                           <  4:147982/148‑1 (MQ=255)
                                                                                           ttCTCAGGAATGCACCGGTTTCGCACGCTCACGCAGAATTGTACGCAGGCATTCGCCGACGGAAGGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGTCAGCTCCTTGGTCTGATTGCCGCAg                                                           <  2:323264/148‑1 (MQ=255)
                                                                                              tcAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCTTTGGTCTGATTGGCGAAGCAg                                                        >  1:162259/1‑148 (MQ=255)
                                                                                                aGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTAc                                                     >  4:321173/1‑149 (MQ=255)
                                                                                                     tGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGGATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCAcgc                                                >  5:360338/1‑149 (MQ=255)
                                                                                                        aCCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGc                                                               <  2:44921/131‑1 (MQ=255)
                                                                                                        aCCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGc                                                               >  1:44921/1‑131 (MQ=255)
                                                                                                            gTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGACTCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTtc                                         <  5:174129/149‑1 (MQ=255)
                                                                                                                    aaGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTcga                                  <  3:82896/149‑1 (MQ=255)
                                                                                                                         cACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGaa                              <  4:358196/148‑1 (MQ=255)
                                                                                                                                     tACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTc                  >  6:127035/1‑148 (MQ=255)
                                                                                                                                        ggAGGCATTCCGTCGGCGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGaa               >  7:34391/1‑148 (MQ=255)
                                                                                                                                              aTTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCTTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCAg        <  2:162259/149‑1 (MQ=255)
                                                                                                                                                     cGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCGTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTAGACGAAGCTCTGCAGGTCGAAATCCCAggaatgg  >  2:334391/1‑144 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
GCAGCTTTGCCGTCTCCCTCGACCAGGTCGATGATCTGCAGTTCTTGGCTCACGGGTACACTCCAAACGCTGAAAAACAAAACCGCCCTTTTCTCAGGAATGCACCGGTTTGGCAAGTTCACGCAGAATTGTACGGAGGCATTCGCCGACGGAATGCCTCCGTACAGACATGACGTTGAATCAGCTGTTCTGCTCGAACACCTCTTGGCCGGACAGCTCCTTGGTCTGATTGGCGAAGCAGTACCACGCGGGTTTCTCTTCGACGAAGATCTGCAGGTCGAAATCCCAGGCTTGG  >  NC_002947/4905354‑4905648

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: