Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 919,675 G→C 10.2% intergenic (‑143/+48) PP_5458 ← / ← PP_0800 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947919,6750GC10.2% 136.8 / 5.9 46intergenic (‑143/+48)PP_5458/PP_0800hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base C (3/2);  total (25/21)
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.83e-01

ATTTGAAACTACCTTCTTGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACATGCAAG  >  NC_002947/919531‑919822
                                                                                                                                                |                                                                                                                                                   
aTTTGAAACTACCTTCTTGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGaaa                                                                                                                                                <  1:756363/150‑1 (MQ=255)
     aaaCTACCTTCTTGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGCCGAAAGcccc                                                                                                                                           <  2:313065/150‑1 (MQ=255)
     aaaCTACCTTCTTGCGGTGCAAGTGATAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGCCGAAAGcccc                                                                                                                                           <  1:394909/150‑1 (MQ=255)
                  gCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGcc                                                                                                                              <  2:310022/150‑1 (MQ=255)
                               tAATGAGGACAGGGGTAAAGTTAATCCTTTCGGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGcgcc                                                                                                                 >  2:490102/1‑150 (MQ=255)
                                 aTGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGAACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGAGCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGcgccgc                                                                                                               >  1:63601/1‑150 (MQ=255)
                                        cAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCa                                                                                                        <  2:141931/150‑1 (MQ=255)
                                         aGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAg                                                                                                       >  2:437577/1‑150 (MQ=255)
                                            ggTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCt                                                                                                    >  2:635366/1‑150 (MQ=255)
                                             gTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTa                                                                                                   >  1:470352/1‑150 (MQ=255)
                                                   ttAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAAcc                                                                                             >  2:605291/1‑150 (MQ=255)
                                                     aGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGa                                                                                           <  1:635366/150‑1 (MQ=255)
                                                       tCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATa                                                                                         >  1:374541/1‑150 (MQ=255)
                                                          tttCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGGAGGCTAGAAACCGATAAAg                                                                                      >  2:524737/1‑150 (MQ=255)
                                                                 gCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGtaat                                                                               <  2:616075/150‑1 (MQ=255)
                                                                   aCTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGtaataa                                                                             >  1:590249/1‑150 (MQ=255)
                                                                    cTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATaaa                                                                            <  2:583081/150‑1 (MQ=255)
                                                                     tCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAAc                                                                           >  1:602445/1‑150 (MQ=255)
                                                                        aaaCAGAGCGAAGTGCGATAAGCACTGTTGCGCGAACGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCgg                                                                        <  2:63601/150‑1 (MQ=255)
                                                                            agtgCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAg                                                                                      <  2:596429/129‑1 (MQ=255)
                                                                            agagCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAg                                                                                      >  1:596429/1‑132 (MQ=255)
                                                                                 gAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCCTGGCAGGCTAGAAACCCGTAAAGGCGTAATAACCCGGGGGTTCATg                                                               <  2:232162/150‑1 (MQ=255)
                                                                                     tGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGtt                                                           >  1:533106/1‑150 (MQ=255)
                                                                                     tGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGtt                                                           <  1:605291/150‑1 (MQ=255)
                                                                                        gATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTtgtg                                                        >  1:90627/1‑150 (MQ=255)
                                                                                         aTAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTgc                                                       <  1:12391/150‑1 (MQ=255)
                                                                                         aTAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTgc                                                       >  1:340492/1‑150 (MQ=255)
                                                                                                      ctcAAAAGGTAACGATCAGTGGGTGGAACAGCGGGGGTTTGGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGc                                          <  2:564558/148‑1 (MQ=255)
                                                                                                         aaaaGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGccag                                       >  1:564558/1‑148 (MQ=255)
                                                                                                               tAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCCGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCt                                 >  2:507156/1‑150 (MQ=255)
                                                                                                                   gATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGgcgc                                                                                                                  >  1:425365/1‑65 (MQ=255)
                                                                                                                   gATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGgcgc                                                                                                                  <  2:425365/65‑1 (MQ=255)
                                                                                                                   gATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGcc                             >  2:190118/1‑150 (MQ=255)
                                                                                                                      cATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTc                          <  1:437577/150‑1 (MQ=255)
                                                                                                                       aTTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAATGCGTAATAAACCGGGGTTTCATGGTTTGTGCTGAAACCCATTGCCTGCAAGCTTGCCGTCa                         <  2:8890/150‑1 (MQ=255)
                                                                                                                            gTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCt                    <  1:190118/150‑1 (MQ=255)
                                                                                                                               gAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCa                 >  2:576621/1‑150 (MQ=255)
                                                                                                                                    acGGGGCTTTCGCCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTg            >  1:290676/1‑150 (MQ=255)
                                                                                                                                    acGGGGCTTTCGCCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTg            >  1:290477/1‑150 (MQ=255)
                                                                                                                                       gggCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGcac         >  2:79200/1‑150 (MQ=255)
                                                                                                                                       gggCTTTCGCCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGcac         >  2:734653/1‑150 (MQ=255)
                                                                                                                                        ggCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTa                                                                                                   >  1:647889/1‑59 (MQ=255)
                                                                                                                                        ggCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTa                                                                                                   <  2:647889/59‑1 (MQ=255)
                                                                                                                                         gCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACAt       <  2:470352/150‑1 (MQ=255)
                                                                                                                                         gCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACAt       >  1:757529/1‑150 (MQ=255)
                                                                                                                                              cGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACATGCAAg  <  2:590249/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATTTGAAACTACCTTCTTGCGGTGCAAGTGCTAATGAGGACAGGGGTAAAGTTAGTCCTTTCAGTGCACTCGAAACAGAGCGAAGTGCGATAAGCACTGTTGCGCAAAAGGTAACGATCATTGGGTGGAACAACGGGGCTTTCGGCGAAAGCCCCGTGCTAGGGCGCCTGCACAGGGCGCCGCCTTGGCAGGCTAGAAACCGATAAAGGCGTAATAAACCGGGGTTTCATGGGTTGTGCGCAAACCCAGTGCCTGCAAGCTTGCCGTCAGCGCTTCATCCTGCACATGCAAG  >  NC_002947/919531‑919822

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