Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,298,936 T→A 9.7% intergenic (+47/+44) PP_3770 → / ← PP_3771 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,298,9360TA 9.7% 121.3 / 3.7 41intergenic (+47/+44)PP_3770/PP_3771hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (2/2);  total (21/20)
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.79e-01

CCTGGTTCGCAACGGCGAAGCCGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCTTGGCGAGCA  >  NC_002947/4298790‑4299080
                                                                                                                                                  |                                                                                                                                                
ccTGGTTCGCAACGGCGAAGCCGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTg                                                                                                                                               >  1:225943/1‑150 (MQ=255)
  tGGTTCGCAACGGCGAAGCCGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGc                                                                                                                                             <  2:516212/150‑1 (MQ=255)
          aaCGGCGAAGCCGTGAGCGGCCCGATGCCTACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCtttt                                                                                                                                     >  1:254228/1‑150 (MQ=255)
          aaCGGCGAAGCCGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCCGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCtttt                                                                                                                                     >  1:651683/1‑150 (MQ=255)
          aaCGGCGAAGCCGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCtttt                                                                                                                                     >  2:141992/1‑150 (MQ=255)
           aCGGCGAAGCCGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAAATTGGCGGGcgacaa                                                                                                                                    <  1:280459/150‑6 (MQ=255)
                    ccGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTcgc                                                                                                                           >  2:282486/1‑150 (MQ=255)
                            ggCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGcc                                                                                                                   >  2:966831/1‑150 (MQ=255)
                              cccGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCg                                                                                                                 >  1:1398743/1‑150 (MQ=255)
                               ccGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCgg                                                                                                                <  1:321418/150‑1 (MQ=255)
                                   tGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGAtt                                                                                                                      <  1:1155200/140‑1 (MQ=255)
                                   tGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGAtt                                                                                                                      >  2:1155200/1‑140 (MQ=255)
                                       gACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGcccc                                                                                                        >  1:282053/1‑150 (MQ=255)
                                         cccGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGcccccc                                                                                                      <  1:966831/150‑1 (MQ=255)
                                                    gCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGcc                                                                                           >  1:1313474/1‑150 (MQ=255)
                                                     cccAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCt                                                                                                     <  2:1311083/139‑1 (MQ=255)
                                                     cccAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCt                                                                                                     >  1:1311083/1‑139 (MQ=255)
                                                            aCTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGc                                                                                   <  1:141992/150‑1 (MQ=255)
                                                               ggCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCa                                                                                <  2:1398743/150‑1 (MQ=255)
                                                                      cAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACt                                                                         <  1:282486/150‑1 (MQ=255)
                                                                                       tttGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGAcct                                                        >  1:719626/1‑150 (MQ=255)
                                                                                          ggTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGAccttgc                                                     <  2:282053/150‑1 (MQ=255)
                                                                                           gTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGAccttgcc                                                    >  1:1080424/1‑150 (MQ=255)
                                                                                              gCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGttt                                                                                               >  1:539894/1‑104 (MQ=255)
                                                                                              gCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGttt                                                                                               <  2:539894/104‑1 (MQ=255)
                                                                                               ccTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACcttgccttgc                                                <  1:934623/150‑1 (MQ=255)
                                                                                                              ttgtGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGtc                                 <  2:1313474/150‑1 (MQ=255)
                                                                                                              ttgtGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGtc                                 >  1:491585/1‑150 (MQ=255)
                                                                                                              ttgtGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGtc                                 >  1:1141660/1‑150 (MQ=255)
                                                                                                                  gCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCATACGGtcatca                             >  1:1287145/1‑150 (MQ=255)
                                                                                                                            cTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGttcgtt                   <  2:1080424/150‑1 (MQ=255)
                                                                                                                             tACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGttcgttt                  >  1:1268135/1‑150 (MQ=255)
                                                                                                                                     tGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCtt          <  2:1268135/150‑1 (MQ=255)
                                                                                                                                     tGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCtt          <  2:225943/150‑1 (MQ=255)
                                                                                                                                        cgcccgccAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCTTGGc       <  2:491585/150‑1 (MQ=255)
                                                                                                                                        agcccgctAAAGTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCTTGGc       >  2:415824/2‑150 (MQ=255)
                                                                                                                                         gcccgccAAATTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGAc                                                          >  1:439097/1‑98 (MQ=255)
                                                                                                                                         gcccgccAAATTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGAc                                                          <  2:439097/98‑1 (MQ=255)
                                                                                                                                          cccgccAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCTTGGCgt     <  2:651683/150‑2 (MQ=255)
                                                                                                                                            cgccAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCTTGGCGAGc   >  2:869722/1‑150 (MQ=255)
                                                                                                                                             gccAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCTTGGCGAGCa  <  2:254228/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCTGGTTCGCAACGGCGAAGCCGTGAGCGGCCCGATGCCGACCCGCGCCTCAGCCCAGGCACTGGCCGAACAGTTGCAACCAGGCACTTTGGTTGCCTGAACCGAATTTGTTGTGCGTTGCTCCCTACGCCCTTGTCGCCCGCCAATTTGGCGGGCTTTTTTTTCTTCGCGGATTGCCCGGGCTGCCCCCCTACGTTTGGCCGCTACTGCCCAGCTTACTGGTCAGCACCATGACCTTGCCTTGCAGCCGGCCTACGGTCATCACCAGTTCGTTTCGTTCCTTGGCGAGCA  >  NC_002947/4298790‑4299080

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