Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,995,018 T→A 14.6% intergenic (+146/+31) PP_5251 → / ← PP_5252 transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,995,0180TA14.6% 110.6 / 11.8 48intergenic (+146/+31)PP_5251/PP_5252transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/24);  new base A (4/3);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.34e-01

CCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACA  >  NC_002947/5994875‑5995156
                                                                                                                                               |                                                                                                                                          
ccGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTg                                                                                                                                        <  7:77039/148‑1 (MQ=255)
    cccTTACGGCTCATTCGCGGGTGAACCCGCTCCTACCCAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGc                                                                                                                                      <  8:147452/146‑1 (MQ=255)
    cccTTACGGCTCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGc                                                                                                                                      >  7:147452/1‑146 (MQ=255)
    cccTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGcccc                                                                                                                                   <  4:122473/149‑1 (MQ=255)
    cccTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGcccc                                                                                                                                   >  2:151353/1‑149 (MQ=255)
     ccTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccct                                                                                                                                  >  6:129120/1‑148 (MQ=255)
          cGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccctttctt                                                                                                                             <  5:129120/149‑7 (MQ=37)
           ggCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGATTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGt                                                                                                                            >  8:92022/1‑149 (MQ=255)
                     cgGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACAtcagc                                                                                                                   >  6:130466/1‑148 (MQ=255)
                       ggTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGccc                                                                                                                                    <  7:150278/129‑1 (MQ=255)
                       ggTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGccc                                                                                                                                    >  8:150278/1‑129 (MQ=255)
                            aCCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGCAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTa                                                                                                           <  1:65039/149‑1 (MQ=255)
                            aCCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTa                                                                                                           <  3:68398/149‑1 (MQ=255)
                              ccGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGtt                                                                                                                                            <  3:44671/114‑1 (MQ=255)
                              ccGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGtt                                                                                                                                            >  4:44671/1‑114 (MQ=255)
                                    cTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGc                                                                                                     <  4:274548/147‑1 (MQ=255)
                                    cTACACAGACCTGCCGAACACACACACCGGCCGTTCGAACCAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGa                                                                                                   <  4:251595/149‑1 (MQ=38)
                                      acacAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAg                                                                                                 >  3:105384/1‑149 (MQ=255)
                                             tCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGaa                                                                                           <  1:151353/148‑1 (MQ=255)
                                                  ggTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGcc                                                                                      >  8:167193/1‑148 (MQ=255)
                                                      gACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAgc                                                                                  <  5:33981/148‑1 (MQ=255)
                                                         acacacCTGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTa                                                                                                           <  4:162339/120‑1 (MQ=255)
                                                             acCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAccc                                                                            >  7:227513/1‑147 (MQ=255)
                                                                 gACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAg                                                                      >  7:205182/1‑149 (MQ=255)
                                                                   cGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGcc                                                                    <  7:92022/149‑1 (MQ=255)
                                                                    gTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACa                                                                                         <  1:121586/127‑1 (MQ=255)
                                                                    gTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACa                                                                                         >  2:121586/1‑127 (MQ=255)
                                                                     ttGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCCCATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCt                                                                  >  6:14937/1‑149 (MQ=255)
                                                                               gCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCAtt                                                        >  3:145243/1‑149 (MQ=255)
                                                                                     cctttcgCGGGTAAACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCgg                                                  >  2:272082/7‑149 (MQ=255)
                                                                                         gAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAACGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCAAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACg                                              >  5:305260/1‑149 (MQ=255)
                                                                                            cGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTc                                           <  3:269104/149‑1 (MQ=255)
                                                                                            cGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTc                                           <  4:105384/149‑1 (MQ=255)
                                                                                                tttACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCTTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAAt                                       <  7:167193/149‑1 (MQ=255)
                                                                                                    cccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGcc                                    <  5:130466/148‑1 (MQ=255)
                                                                                                      cgcgAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCcgc                                  <  3:255992/148‑1 (MQ=255)
                                                                                                      cgcgAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCcgc                                  >  4:255992/1‑148 (MQ=255)
                                                                                                               gCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGt                        <  4:145243/149‑1 (MQ=255)
                                                                                                                ccGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCCGCTTAGAAAGCGAAGAACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTg                       <  1:260937/149‑1 (MQ=255)
                                                                                                                ccGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTg                       >  5:286938/1‑149 (MQ=255)
                                                                                                                ccGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTg                       <  6:286938/149‑1 (MQ=255)
                                                                                                                         aaGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTgc                <  5:14937/147‑1 (MQ=255)
                                                                                                                              acaAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGca         >  4:292935/1‑149 (MQ=255)
                                                                                                                                 aGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTc                          <  6:286720/129‑1 (MQ=255)
                                                                                                                                 aGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTc                          >  5:286720/1‑129 (MQ=255)
                                                                                                                                  gAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCAcacc     <  8:227513/149‑1 (MQ=255)
                                                                                                                                  aacgggggCCAATAACTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCa                                                          <  6:264640/92‑1 (MQ=255)
                                                                                                                                  aacgggggCCAATAACTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCa                                                          >  5:264640/5‑96 (MQ=255)
                                                                                                                                    cgggggCCAATAACTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCAcacca    >  6:171461/3‑148 (MQ=255)
                                                                                                                                     aGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACcaca  <  6:305260/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                          
CCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCACA  >  NC_002947/5994875‑5995156

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: ‑