Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,044,739 A→G 7.8% intergenic (+45/+223) PP_4448 → / ← PP_5681 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,044,7390AG 7.8% 141.3 / 2.7 51intergenic (+45/+223)PP_4448/PP_5681hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base G (2/2);  total (23/28)
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.12e-01

CGATCCAGCACCATCGTCTGCCGACATTCCATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCAAGCTCTAGATCCAGCCTTCTTCACA  >  NC_002947/5044592‑5044885
                                                                                                                                                   |                                                                                                                                                  
cGATCCAGCACCATCGTCTGCCGACATTCCATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCAc                                                                                                                                                   <  7:86461/149‑1 (MQ=255)
 gATCCAGCACCATCGTCTGCCGACATTCCATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCGGCGcc                                                                                                                                                  >  5:353701/1‑149 (MQ=255)
     cAGCACCATCGTCTGCCGACATTCCATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTaa                                                                                                                                               <  1:171796/148‑1 (MQ=255)
          ccATCGTCTGCCGACATTCCATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAg                                                                                                                                         <  5:298083/149‑1 (MQ=255)
                   gCCGACATTCCATACTGTGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGAGcc                                                                                                                                >  4:235836/1‑149 (MQ=255)
                          ttCCATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGt                                                                                                                         <  1:12354/149‑1 (MQ=255)
                            ccATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGtt                                                                                                                        <  4:177571/148‑1 (MQ=255)
                             cATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCGGCGCCTTCCTGAAGTTAGGTGCGGGGCCGttgg                                                                                                                      <  6:353701/149‑3 (MQ=37)
                               tACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAgc                                                                                                                    >  8:107550/1‑149 (MQ=255)
                                 cTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAgcgc                                                                                                                  >  7:250722/1‑149 (MQ=255)
                                       aTTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTc                                                                                                                >  5:297060/1‑145 (MQ=255)
                                          gtgGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATc                                                                                                          <  4:265205/148‑1 (MQ=255)
                                              aaggaagAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGACGCCGGGCCGTTCAGCGCTCAAGATcagaag                                                                                                     >  7:102166/1‑149 (MQ=255)
                                                 gaagaaTGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATcagaagcca                                                                                                  <  8:250722/149‑1 (MQ=255)
                                                      aTGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCGGCGCCTTCCTGAAGTTAGGTGCGGGGCCGttg                                                                                                                       >  8:27117/1‑122 (MQ=25)
                                                      aTGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCGGCGCCTTCCTGAAGTTAGGTGCGGGGCCGttg                                                                                                                       <  7:27117/123‑2 (MQ=25)
                                                       tGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCACACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCa                                                                                            <  8:130800/149‑1 (MQ=255)
                                                             gATCTTAGCTTTGATCAGATCACGACTTAGTGGATTAAATAAGCACTTCTGCTCATCTGATAAACTAATCGCACCAACGGAACCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCAtttttt                                                                                      <  5:353662/149‑1 (MQ=255)
                                                               tCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTg                                                                                     <  7:126867/148‑1 (MQ=255)
                                                                    gCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATcagaagccagaa                                                                                               >  8:145913/1‑133 (MQ=255)
                                                                    gCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATcagaagccagaa                                                                                               <  7:145913/133‑1 (MQ=255)
                                                                     cTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGtt                                                                                                                        <  6:5604/107‑1 (MQ=255)
                                                                     cTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGtt                                                                                                                        >  5:5604/1‑107 (MQ=255)
                                                                       ttGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATcagaagccag                                                                                                 >  1:206109/1‑128 (MQ=255)
                                                                       ttGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATcagaagccag                                                                                                 <  2:206109/128‑1 (MQ=255)
                                                                        tGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGACGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAg                                                                             <  8:102166/147‑1 (MQ=255)
                                                                            cagatcaCGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTc                                                                       >  3:2163/1‑149 (MQ=255)
                                                                               atcaCGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCg                                                                     >  6:207868/1‑148 (MQ=255)
                                                                               atcaCGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCg                                                                     >  4:92533/1‑148 (MQ=255)
                                                                                         aGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGc                                                            <  5:207868/147‑1 (MQ=255)
                                                                                             gCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGc                                                      >  1:345851/1‑149 (MQ=255)
                                                                                                    tAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGt                                                                    <  1:347739/128‑1 (MQ=255)
                                                                                                    tAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGt                                                                    >  2:347739/1‑128 (MQ=255)
                                                                                                        cACTGCTGCTGATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTGAGGAAGGCGCCGGGCCGTTCAGCGCTCAATATCATAAGGCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGTGCTCCTCCGTcgc                                           <  6:297060/149‑1 (MQ=255)
                                                                                                        cACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTcgc                                           <  3:92533/149‑1 (MQ=255)
                                                                                                                 tCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTAAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGATTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCt                                  >  1:85825/1‑149 (MQ=255)
                                                                                                                   aTCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGa                                <  2:345851/149‑1 (MQ=255)
                                                                                                                        agaAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCa                           <  4:2163/149‑1 (MQ=255)
                                                                                                                         gaAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTc                                     >  6:367368/1‑138 (MQ=255)
                                                                                                                         gaAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTc                                     <  5:367368/138‑1 (MQ=255)
                                                                                                                         gaAACGAAGCCCCCCTACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGAATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCaa                          <  4:38004/149‑1 (MQ=255)
                                                                                                                          aaaCGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCGGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGCGCTCCTCCGTCGCGCATTCCCTGATTCCAAg                         >  1:252878/1‑149 (MQ=255)
                                                                                                                            aCGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCAAGct                       >  5:81803/1‑149 (MQ=255)
                                                                                                                                    cACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGc                                                      >  6:53209/1‑110 (MQ=255)
                                                                                                                                    cACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGc                                                      <  5:53209/110‑1 (MQ=255)
                                                                                                                                       cAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGGTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGCTACCAAGCTCTAGATCCAGc            >  6:328295/1‑149 (MQ=255)
                                                                                                                                           ggCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCAAGCTCTAGATCCAGCcttc        >  3:342748/1‑149 (MQ=255)
                                                                                                                                              cccGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCAAGCTCTAGATCCAGCCttcttc     <  3:235836/149‑1 (MQ=255)
                                                                                                                                               ccGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTc                                     >  7:190994/1‑116 (MQ=255)
                                                                                                                                               ccGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTc                                     <  8:190994/116‑1 (MQ=255)
                                                                                                                                                 gCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCAAGCTCTAGATCCAGCCTTCTTcaca  <  7:107550/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGATCCAGCACCATCGTCTGCCGACATTCCATACTGGGAATTGTGGAAGGAAGAATGGAAAGATCTTAGCTTTGATCAGATCACGACTGAGTGGCTTAAATAAGCACTGCTGCTCATCTGAGAAACGAAGCGCACCAACGGCCCCGCACCTAACTTCAGGAAGGCGCCGGGCCGTTCAGCGCTCAAGATCAGAAGCCAGAAGCATTTTTTGGAAGGTAGGGCTTCCGTAATGTTGCCAGAGCTCCTCCGTCGCGCATTCCCTGATACCAAGCTCTAGATCCAGCCTTCTTCACA  >  NC_002947/5044592‑5044885

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