Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,191,109 A→T 10.0% intergenic (‑441/‑716) PP_5612 ← / → PP_5613 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,191,1090AT10.0% 132.0 / 4.0 47intergenic (‑441/‑716)PP_5612/PP_5613hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/17);  new base T (2/3);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCGAGCTTATCCAACATGGCGCTCCCCACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGTTCTTCAATCATTATGTG  >  NC_002947/4190969‑4191257
                                                                                                                                             |                                                                                                                                                    
cTGCCGAGCTTATCCAACATGGCGCTCCCCACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCTAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGttt                                                                                                                                              <  4:136091/149‑1 (MQ=255)
cTGCCGAGCTTATCCAACATGGCGCTCCCCACGTTATTTTTGGAAACGTTATCTTCGGATACAATGAAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGttt                                                                                                                                              <  1:41529/149‑1 (MQ=255)
    cGAGCTTATCCAACATGGCGCTCCCCACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTTTCAAAGtttttga                                                                                                                                          <  3:138406/149‑5 (MQ=255)
         ttATCCAACATGGCGCTCCCCACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAaat                                                                                                                                      >  7:96746/1‑148 (MQ=255)
            tCCAACATGCCGCTCCCCACGTTTTTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTTAAAAATTTTGGGCAAATAAtt                                                                                                                                  <  4:36077/149‑1 (MQ=255)
                     gcgcTCCCCACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATaaa                                                                                                                         >  4:68780/1‑149 (MQ=255)
                           cccACGTTATTTTTGGAAACTTTATCTTCTGATACATTTCAAGATCAACCCACACCAAAAACAACTTAGAGCAACGCTACTGCGAGAGCAAGTGTTTTACATCA‑ATAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAATTTGAc                                                                                                                    <  4:210007/148‑1 (MQ=255)
                           cccACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTTTCAAAGtttttga                                                                                                                                          >  8:41203/1‑122 (MQ=255)
                           cccACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTTTCAAAGtttttga                                                                                                                                          <  7:41203/126‑5 (MQ=255)
                           cccACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTaa                                                                                                                                >  5:218984/1‑136 (MQ=255)
                           cccACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTaa                                                                                                                                <  6:218984/136‑1 (MQ=255)
                           cccACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGAc                                                                                                                    <  3:171/148‑1 (MQ=255)
                           cccACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACg                                                                                                                   <  8:96746/149‑1 (MQ=255)
                                gTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATAAACCCACAACAAAAACAACTTAGACCAACGCTACTGCGAGAGAAAGTGTTTTCAATAA‑AAAAATTTGAAAAAGTTTGGGCAAATAATTAAATAAAAATTTACGGataa                                                                                                              >  5:10361/1‑149 (MQ=255)
                                  tATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACAAAAGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGataata                                                                                                            >  3:246476/1‑149 (MQ=255)
                                      tttGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATcc                                                                                                         <  6:158488/148‑1 (MQ=255)
                                          gAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACt                                                                                                    >  6:277472/1‑149 (MQ=255)
                                                    cTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCAc                                                                                          >  4:150189/1‑149 (MQ=255)
                                                     ttCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCAct                                                                                         >  1:200506/1‑149 (MQ=255)
                                                          aTACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCActctctc                                                                                    <  3:68780/149‑1 (MQ=255)
                                                              aaTGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATATTTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGa                                                                                 >  2:232505/1‑148 (MQ=255)
                                                              aaTGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCGATACTCATTTACCACTCTCTCTGat                                                                                >  8:129096/1‑149 (MQ=255)
                                                                      aTCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAAt                                                                        >  2:90146/1‑149 (MQ=255)
                                                                       tCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAAtt                                                                       <  5:277472/149‑1 (MQ=255)
                                                                        cAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAAttt                                                                      >  1:135936/1‑149 (MQ=255)
                                                                          aCCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCATAAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTc                                                                     >  1:242094/1‑149 (MQ=255)
                                                                                  cAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTaa                                                              >  5:291190/1‑147 (MQ=255)
                                                                                       acaacTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGTAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAgcg                                                         >  7:64612/1‑147 (MQ=255)
                                                                                       acaacTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCgca                                                       <  6:15160/149‑1 (MQ=255)
                                                                                                  cAACGCTACTTCTAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACTGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGc                                            <  6:10361/149‑1 (MQ=255)
                                                                                                  cAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGc                                            >  8:149897/1‑149 (MQ=255)
                                                                                                    aCGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATACTTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCt                                          <  1:232505/149‑1 (MQ=255)
                                                                                                      gCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTa                                        <  3:150189/149‑1 (MQ=255)
                                                                                                        tACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCgcagc                                                     >  6:181174/1‑134 (MQ=255)
                                                                                                        tACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGa                                      >  3:229477/1‑149 (MQ=255)
                                                                                                             cGAGAGAAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTTACGGATAATATACATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGAcccc                                  >  3:281852/1‑148 (MQ=255)
                                                                                                              gagaGCAAGTGTTTTCCATCA‑AAAACTTTGTAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCaa                                <  8:64612/149‑1 (MQ=255)
                                                                                                                         ttttCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGcc                      >  5:139484/1‑148 (MQ=255)
                                                                                                                          tttCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCt                     <  6:139484/148‑1 (MQ=255)
                                                                                                                           ttCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGt                   >  2:139470/1‑149 (MQ=255)
                                                                                                                           ttCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGt                   >  8:175681/1‑149 (MQ=255)
                                                                                                                               aTCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGttctt               <  2:200506/149‑1 (MQ=255)
                                                                                                                                  a‑aaaaCTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGttctt               <  6:47373/146‑1 (MQ=255)
                                                                                                                                  a‑aaaaCTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGttctt               >  5:47373/1‑146 (MQ=255)
                                                                                                                                    aaaaCTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGTGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGTTCTTCAAt           >  2:173950/1‑149 (MQ=255)
                                                                                                                                    aaaaCTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGTTCTTCaa            >  3:244505/1‑148 (MQ=255)
                                                                                                                                     aaaCTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGTTCTTCAAt           >  6:54658/1‑148 (MQ=255)
                                                                                                                                             aaaaaGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTTATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGTTCTTCAATCATTAtgtg  >  8:144148/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CTGCCGAGCTTATCCAACATGGCGCTCCCCACGTTATTTTTGGAAACGTTATCTTCGGATACAATGCAAGATCAACCCACACCAAAAACAACTTAGACCAACGCTACTGCGAGAGCAAGTGTTTTCCATCA‑AAAACTTTGAAAAAGTTTGGGCAAATAATTAACTATAAATTGACGGATAATATCCATACTCATTTACCACTCTCTCTGATATTTGAATTTCGAGTTAACAGCGCAGCATCAACTGCCTTAGACCCCAAGAAATCCGCCTGTTCTTCAATCATTATGTG  >  NC_002947/4190969‑4191257

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