Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,466,841 T→G 16.6% intergenic (‑16/+108) PP_3081 ← / ← PP_5560 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,466,8410TG16.6% 76.4 / 8.9 47intergenic (‑16/+108)PP_3081/PP_5560hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (9/27);  new base G (6/2);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 1.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

TGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTC  >  NC_002947/3466700‑3466980
                                                                                                                                             |                                                                                                                                            
tGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGcggcgg                                                                                                                                       >  2:81970/1‑149 (MQ=255)
   gggCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGCAGGCGGCCaa                                                                                                                                   <  3:221117/149‑1 (MQ=255)
      cAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACg                                                                                                                                 <  5:134782/149‑1 (MQ=255)
          ttGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCtgt                                                                                                                             >  6:139611/1‑149 (MQ=255)
              ccccGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTGGCAGGGCAATGACGCCGCCATCAGTGAAAGCCGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGGGCAg                                                                                                                        <  3:203542/149‑1 (MQ=255)
                 aGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCAGCGGCGGCCAACGCTGTGGACGgc                                                                                                                     >  7:60450/1‑149 (MQ=255)
                  gCTCAGTTGCGGGGTCCGCTCTTCACCGGAACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGTAACTGACCTAAATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGgca                                                                                                                    <  6:205232/149‑1 (MQ=255)
                   cTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGgcag                                                                                                                   <  8:6180/149‑1 (MQ=255)
                      aGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGCAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGgcagca                                                                                                                 <  4:51768/148‑1 (MQ=255)
                        ttGGGGGGTCCGCTCTTCACCGGCACCGGATGCCGGGTCCTGGCCTACGTGTTGCAGGGAAATTACGCCGCCAACTTTGAAAGACGAGTAATTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAAAGCTGTGCAGGGCAGCAAag                                                                                                              <  2:260755/149‑1 (MQ=255)
                         tGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCCCCTGGGGCGGGCAACGCTGTGGAGGGCCCCaaacc                                                                                                             >  3:147090/1‑147 (MQ=255)
                          gCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAagc                                                                                                             <  5:139611/148‑1 (MQ=255)
                                tCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAAc                                                                                                      <  4:110570/149‑1 (MQ=255)
                                  cGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGCACAg                                                                                                    >  8:194799/1‑149 (MQ=255)
                                             ggCACCGCATGCCGGTTCCTGGCCCACCTGTTGCAGGCCAATGACGCCGCCATCATTTAAAGACTAGTAAGTCACCTACATCGCCCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGa                                                                                         <  4:147090/149‑1 (MQ=255)
                                             ggCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGa                                                                                         <  1:81970/149‑1 (MQ=255)
                                                 tgcggtGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGt                                                                                      <  1:53394/144‑1 (MQ=255)
                                                            gTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGCCGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATc                                                                          <  7:287351/149‑1 (MQ=255)
                                                               cTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGGCCACGCGGTGCCGGGCAGCACACCACCGAACACCCAGAGTCGGTGGTGCATGGTGAATACtt                                                                       >  8:123683/1‑149 (MQ=255)
                                                                  gCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAg                                                                                                >  5:322340/1‑121 (MQ=255)
                                                                  gCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAg                                                                                                <  6:322340/121‑1 (MQ=255)
                                                                   ccGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCa                                                                   <  4:59948/149‑1 (MQ=255)
                                                                          gTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCCCCAGCGGCGGCCNACCCGGTGCAGGGGAGCAAAGCAGCGGACAGCCCGTGGGGGAGGTGGAAGGTAAATGCTTGGGAGGGCGt                                                             >  1:104454/1‑148 (MQ=255)
                                                                           ttGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGGCGCCGGCGGGGGCCCCCGCCGTGCCGGGCCGGCAAGAAGCGCACCCCCAGTGTGTGGGGTGGATGGTGAATGCGTGGGAggggggga                                                           >  2:308842/1‑144 (MQ=255)
                                                                             gCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGGGGCGGCCAAAGCTGTGCACGGCAGCAAAGCAGCGCACAGCCAGGGGCGGAGGTGCAAGGTCAATCCATGGCAGGGCGTCAcc                                                         >  7:193921/1‑148 (MQ=255)
                                                                              cAGGGCAATGACGCGGCCATCAGTCAACCAGGAGGACCGGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCATAGCAGCGAACAGCCAGTGTCTGAGGTGCTTGGTCAATCCTTGGCAGGGCGTCACgg                                                        <  7:222524/149‑1 (MQ=255)
                                                                                gggCAATGACGCCGCCATCAGTGAAACACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGGCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCt                                                      <  8:90852/149‑1 (MQ=255)
                                                                                  gCAATGACGCCGCCATCAGTGAATGACGAGGAACTCACCGACATCGCGCCTGGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCa                                                    <  5:198908/149‑1 (MQ=255)
                                                                                            cgccATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGCCAACCCAGTGCAGGGCCGCAAAGCAGCGAACAGCCAGTGTGTGAGGTGCATGGTCAATCCATGGCAGGGCGGGACGGCGCAGTCTAGCCCa                                          >  6:7012/1‑149 (MQ=255)
                                                                                              ccATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAg                                         <  1:170806/148‑1 (MQ=255)
                                                                                               cATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAggg                                                                <  7:273369/124‑1 (MQ=255)
                                                                                                    gTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCCCCAGCGGCGGCCCACGCGGTGGAGGGCAGCAAAGCCGCGCACAGCCAGGGTCTGCGGTGGATGGGCAAACCTTGGCCGGGGGTCCCGGCTGCGGGTAGCCCCGGTGTGca                                  >  6:62245/1‑148 (MQ=255)
                                                                                                         aGACGAGGAACTGACCGACATCGCGCCTTGGCGGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAACGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCTTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGc                             <  1:209075/149‑1 (MQ=255)
                                                                                                                  aCTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAg                                                   >  7:150066/1‑118 (MQ=255)
                                                                                                                  aCTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAg                                                   <  8:150066/118‑1 (MQ=255)
                                                                                                                  aCTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGCCAACCCTGTGCATGGCAGCAAAGCACCGAACAGCCCGTGTCTGTGGTGCATGGTGAATCCGTGGCAGGGCGTCACGGGGCAGTCTACCCCAGGGGAGGTTGACCGGCTGTccg                    >  2:30786/1‑148 (MQ=255)
                                                                                                                      aCCGACATCCCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTc                <  2:104454/149‑1 (MQ=255)
                                                                                                                        cGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCt              <  4:307360/149‑1 (MQ=255)
                                                                                                                           cATCGGGCCTTGGCCGCGTGGGGCGGCCCACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGGCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCgg           <  5:188926/149‑1 (MQ=255)
                                                                                                                               gcgcGTTGGCCGCCTGCGGGGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGGCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAg        <  7:123683/148‑1 (MQ=255)
                                                                                                                                 gcgttGGCCGCCGCAGGCGGCCAAGGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGc     >  8:13537/4‑149 (MQ=255)
                                                                                                                                  ccTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGt       >  3:69403/1‑146 (MQ=255)
                                                                                                                                  ccTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGt       <  4:69403/146‑1 (MQ=255)
                                                                                                                                    ttGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc  <  8:193921/149‑1 (MQ=255)
                                                                                                                                    ttGGCCGCCGCAGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc  >  3:315733/1‑149 (MQ=255)
                                                                                                                                    ttGGCCGCCGCAGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc  >  8:3434/1‑149 (MQ=255)
                                                                                                                                      ggCCGCCGGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTc  <  7:201061/147‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
TGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTC  >  NC_002947/3466700‑3466980

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