Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,538,569 A→T 16.4% intergenic (+59/+62) PP_4871 → / ← PP_4872 lysine decarboxylase family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,538,5690AT16.4% 80.6 / 8.2 41intergenic (+59/+62)PP_4871/PP_4872lysine decarboxylase family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/18);  new base T (3/4);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGT  >  NC_002947/5538426‑5538716
                                                                                                                                               |                                                                                                                                                    
gcagcaGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCaa                                                                                                                                                 >  3:172224/1‑149 (MQ=255)
 cagcagCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAt                                                                                                                                                <  4:274646/149‑1 (MQ=255)
    cagcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGGCACAAGGCCGCTACCACAAAGGCAATtgg                                                                                                                                             >  5:283521/1‑149 (MQ=255)
     agcCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATgggg                                                                                                                                            <  8:35101/149‑1 (MQ=255)
                   ctgGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTtgt                                                                                                                                <  5:184816/147‑1 (MQ=255)
                      gAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcgc                                                                                                                           <  6:108454/149‑1 (MQ=255)
                      gAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTcgc                                                                                                                           <  4:172224/149‑1 (MQ=255)
                         gggATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGAt                                                                                                                        <  8:216566/149‑1 (MQ=255)
                           gagggACTGGTTTGGGGCGCGGGTCGCGCCTATGTGGGTCTACAAGCACCCTCACTGAAGCTGAGACTTTTGGGACCGCTTCGGGGTCTTTTCGCGACACAAGGCCGCCCCCACAACGGCAATGGGGGAGCGGCTTTGTGTCGCGATgg                                                                                                                      <  4:155937/146‑1 (MQ=255)
                           gATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATgg                                                                                                                      >  3:155937/1‑149 (MQ=255)
                                    cTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAATGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAag                                                                                                             >  5:274032/1‑149 (MQ=255)
                                     tttGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGa                                                                                                                         >  3:184177/1‑136 (MQ=255)
                                     tttGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGa                                                                                                                         <  4:184177/136‑1 (MQ=255)
                                             gcCGGTCGCGCCTATGTGGGTCTACAAGCAGTGTGACTGATGCTGATACTTTTGGTACGGGTTCGCGGTCCTTTCGCGACACTAGGCCGCTCCCACATAGTCAATGGGGGAGGGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCaa                                                                                                    <  8:63931/149‑1 (MQ=37)
                                              ccGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCaa                                                                                                    >  7:156173/1‑148 (MQ=255)
                                                        cTAAGTGGGTCGTCAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAAttg                                                                                                                                              <  5:165276/96‑3 (MQ=255)
                                                        cTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAg                                                                                         >  7:252434/1‑149 (MQ=255)
                                                        cTAAGTCGGTCGACAAGCAGCCTGACTGAACCTGAGACTTTTGGGACCGCTTCGCGGTCCTTACGCCACACAAGGCCGCTCCCACAAAGGCAATgg                                                                                                                                              >  6:165276/1‑96 (MQ=38)
                                                         tAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                                                                  <  6:303150/107‑1 (MQ=255)
                                                         tAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCtt                                                                                                                                  >  5:303150/1‑107 (MQ=255)
                                                               ggTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAAcgc                                                                                  <  8:206438/149‑1 (MQ=255)
                                                                         ctgcattacttaAGGTTAGACTTTTGGGACCGGTTGGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGGGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAg                                                                        <  6:283521/142‑1 (MQ=38)
                                                                         ccgcCTGACGGGAGCGGATCCTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAg                                                                        <  6:132979/147‑1 (MQ=38)
                                                                          agcCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAgg                                                                       <  2:270785/149‑1 (MQ=255)
                                                                           gcCTGACTGAAGCTGAGACTTTTGGGACAGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGc                                                                      >  6:163830/1‑149 (MQ=255)
                                                                                actgaAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAATGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGAc                                                                  <  6:274032/148‑1 (MQ=255)
                                                                                                 tGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGccc                                                <  5:100402/148‑1 (MQ=255)
                                                                                                   gcaCCGATTCGCGGTCCTTTCGCCACACCCGCCCGCTCCCACAATGGCATGGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCtt                                              <  8:156173/146‑1 (MQ=38)
                                                                                                      ccGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCa                                          >  8:116746/1‑149 (MQ=255)
                                                                                                         cTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGCCGCGGCT‑TCTTGTTGCCCTTGCCAccc                                       >  2:241594/1‑149 (MQ=255)
                                                                                                         cTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCAccc                                       <  7:116746/149‑1 (MQ=255)
                                                                                                               ggTCCTTTCGCGACACAAGGCCGCACCCAAAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTc                                 <  7:267776/149‑1 (MQ=255)
                                                                                                               ggTCCTTTCGCGAAACACGGCCGCTCCCACAATGGCAAGGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTc                                 <  8:252434/149‑1 (MQ=255)
                                                                                                                    tttCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCatc                            >  5:23211/1‑149 (MQ=255)
                                                                                                                             acaAGGCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTT‑CCCTTGCCACCCTTGCTCGCCTCTTTGCGCTcc                  >  2:170492/1‑149 (MQ=37)
                                                                                                                                  gCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGct‑tctt                                                       <  8:255529/108‑1 (MQ=25)
                                                                                                                                  gCCGCTCCCCCATTGCCTTTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGct‑tctt                                                       >  7:255529/1‑108 (MQ=25)
                                                                                                                                    cGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGCCAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTgg              >  3:264636/1‑147 (MQ=255)
                                                                                                                                          cacaAAGGCAATGGGGGAGCGGCGTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgc      >  6:157167/1‑149 (MQ=255)
                                                                                                                                           acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCTGTCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgc      <  1:241594/149‑1 (MQ=255)
                                                                                                                                           acaAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCctgct     >  1:214584/1‑149 (MQ=255)
                                                                                                                                             aaaGGCAATGGGGGTGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACTCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTgg   <  4:166854/149‑1 (MQ=255)
                                                                                                                                             aaaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTgg   >  8:234296/1‑149 (MQ=255)
                                                                                                                                              aaGGCAATGGGGGAGCGGCGTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt  <  5:157167/149‑1 (MQ=255)
                                                                                                                                              aaGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGt  <  5:245613/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCAGCAGCCGGACGAGCTGCTGGAAGGGATGGACAGCTTTGTGGCGCCGGTCGCGCCTAAGTGGGTCGACAAGCAGCCTGACTGAAGCTGAGACTTTTGGGACCGCTTCGCGGTCCTTTCGCGACACAAGGCCGCTCCCACAAAGGCAATGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGCCCCAATAGTCGATCAGTTAACGCGGAATGACAGGCTGACGCGGCT‑TCTTGTTGCCCTTGCCACCCTTGCTCGCCTCTTTGCGCTCCTTGGCAGCCTGCTGGT  >  NC_002947/5538426‑5538716

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