Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,104,110 T→A 10.9% intergenic (+32/‑82) PP_5353 → / → PP_5354 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,104,1100TA10.9% 130.7 / 7.6 47intergenic (+32/‑82)PP_5353/PP_5354hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base A (2/3);  total (22/25)
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.71e-01

GAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCTTCCCG  >  NC_002947/6103967‑6104258
                                                                                                                                               |                                                                                                                                                    
gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                                                                                                                >  2:163839/1‑150 (MQ=255)
gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                                                                                                                <  1:661551/150‑1 (MQ=255)
 aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                                                                                                               <  1:101765/150‑1 (MQ=255)
 aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                                                                                                               <  2:133374/150‑1 (MQ=255)
 aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                                                                                                               <  1:589306/150‑1 (MQ=255)
       tCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCGGGGttctg                                                                                                                                         <  1:666116/150‑4 (MQ=255)
         gcgcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttctgt                                                                                                                                        <  2:641719/149‑5 (MQ=255)
         gcgcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttctgt                                                                                                                                        >  1:641719/1‑145 (MQ=255)
           gcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCGGGGt                                                                                                                                             <  2:76580/142‑1 (MQ=255)
           gcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCGGGGt                                                                                                                                             >  1:76580/1‑142 (MQ=255)
           gcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGAGTGCTAAACAGAGAACCCCGGGCTTGCCCGGGGGttttgctt                                                                                                                                     >  1:507295/1‑150 (MQ=255)
               cGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttttgc                                                                                                                                 <  1:999778/150‑1 (MQ=255)
                  ccAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCaa                                                                                                                              >  2:960318/1‑150 (MQ=255)
                  ccAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGGTTTTGCTTTTGCCaa                                                                                                                              >  1:941200/1‑150 (MQ=255)
                  ccAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCAGGCTTGCCCGGGGTTTTTGCTTTTGCCaa                                                                                                                              >  2:90805/1‑150 (MQ=255)
                   cAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCaaa                                                                                                                             >  2:939213/1‑150 (MQ=255)
                   cAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGAGAGAAAAAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCaaa                                                                                                                             >  2:1122153/1‑150 (MQ=255)
                     gggCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGt                                                                                                                           <  2:815378/150‑1 (MQ=255)
                     gggCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGt                                                                                                                           >  2:1067759/1‑150 (MQ=255)
                              gtTGAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGACCCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctacta                                                                                                                  <  2:130502/150‑1 (MQ=255)
                               ttCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGTGGTTTTTGCTTTTGCCAAAGTGGGCTACTAg                                                                                                                 >  2:580432/1‑150 (MQ=255)
                                tCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGAGTGCTAAACACAGATCCCCGGTCTTGTCCGGGGTTGTTGCTTTTGCCAAAGTGGGCTACGAgg                                                                                                                <  2:710087/150‑1 (MQ=255)
                                tCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAAGCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTGGCCAAAGTGGGCTACTAgg                                                                                                                <  2:159980/150‑1 (MQ=255)
                                                 aTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGc                                                                                                                       <  1:465573/126‑1 (MQ=255)
                                                 aTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGc                                                                                                                       >  2:465573/1‑126 (MQ=255)
                                                 aTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt                                                                                               <  2:994967/150‑1 (MQ=255)
                                                       cgcgCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCTGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATc                                                                                         >  2:249111/1‑150 (MQ=255)
                                                        gcgcCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCt                                                                                        <  1:939213/150‑1 (MQ=255)
                                                                         cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGaa                                                                       >  2:944781/1‑150 (MQ=255)
                                                                                ggCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGGACCCCGGGCTTGCCCGGGGtttt                                                                                                                                          <  2:521661/76‑1 (MQ=255)
                                                                                ggCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGGACCCCGGGCTTGCCCGGGGtttt                                                                                                                                          >  1:521661/1‑76 (MQ=255)
                                                                                   tGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGgaaga                                                             >  2:1114807/1‑150 (MQ=255)
                                                                                     gTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGaagaag                                                           <  1:684262/150‑1 (MQ=255)
                                                                                      tAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGt                                                          <  1:960318/150‑1 (MQ=255)
                                                                                        gCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGa                                                        <  1:1122153/150‑1 (MQ=255)
                                                                                         ccGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATCAGGCGCAGTGGGGTGATCTATGTCAGAGGAAGTGAACCGTGGAAGAAGTGAt                                                       >  2:68036/1‑150 (MQ=255)
                                                                                               tCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTCCGGCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaaca                                                 <  1:68036/150‑1 (MQ=255)
                                                                                                             tAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGa                                   >  2:286689/1‑150 (MQ=255)
                                                                                                              gaTCGTTGCGTGCTAAAGACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGGGATCGAACAACTCCGTGAAGCGaa                                  <  2:371304/149‑1 (MQ=255)
                                                                                                               aTCGCTGCGTGCTAAACACAGAACCCCAGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAAt                                 <  1:90805/150‑1 (MQ=255)
                                                                                                                   cTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGc                             <  1:1114807/150‑1 (MQ=255)
                                                                                                                          cTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGa                                                        >  1:711039/1‑116 (MQ=255)
                                                                                                                          cTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGa                                                        <  2:711039/116‑1 (MQ=255)
                                                                                                                                acaGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGcc                >  1:405373/1‑150 (MQ=255)
                                                                                                                                       cccGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAg         >  1:1095087/1‑150 (MQ=255)
                                                                                                                                         cGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCt       >  1:663128/1‑150 (MQ=255)
                                                                                                                                             cTTGCTCGGGGTTTTTGCTTTTGACAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGATGACGTGACCCGTGGAAGAAGTGATCGAACAACTCCGTGTATCGAATGAGCCGGTACCGGTGCCTCTGGAGCTTccc   <  2:196018/150‑1 (MQ=255)
                                                                                                                                              ttGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGATTCCcg  >  1:1078045/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCTTCCCG  >  NC_002947/6103967‑6104258

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