Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,063,193 A→C 16.7% intergenic (+50/‑157) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,063,1930AC16.7% 83.3 / 12.6 36intergenic (+50/‑157)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/19);  new base C (3/3);  total (14/22)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGG  >  NC_002947/2063050‑2063341
                                                                                                                                                 |                                                                                                                                                    
ggAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACAATCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGc                                                                                                                                                  <  1:102675/149‑1 (MQ=255)
        tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCT‑CCCACAGGCATTGCACAGGGt                                                                                                                                           >  4:4991/1‑147 (MQ=255)
            gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGAc                                                                                                                                     >  8:315773/1‑149 (MQ=255)
               gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CGCACAGGACCTGCACAAGGTCCGAAccct                                                                                                                                  <  8:248108/149‑3 (MQ=37)
                cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTACCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                  <  1:31365/149‑1 (MQ=255)
                cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTtg                                                                                                                                 <  3:172047/149‑1 (MQ=255)
                   gAAGCTTATATCCAGCCCGCCACCAAGAGCCTGCTGGATAACACCGCCGTCGTGTTGCACGGCAC‑ACCGAAGTAAACCGGGGACGCTTTGCGGCCCATTCGCTGTCTTGCCCGCT‑ACCACATGCATTTCACAGTGTTCGTACCTtttg                                                                                                                               <  4:44669/148‑1 (MQ=255)
                   gAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGTGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                              <  3:173503/149‑1 (MQ=255)
                      gCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAAGGTCCGAACCCTGTGCaat                                                                                                                           <  2:303238/149‑3 (MQ=255)
                               cAGACCGAGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGc                                                                                                                   <  5:192137/148‑1 (MQ=255)
                                aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACGGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                 <  2:328573/149‑1 (MQ=21)
                                aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGa                                                                                                                 <  1:309951/149‑1 (MQ=21)
                                 gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTTGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTTCAGGTCCTGTGCGAg                                                                                                                <  6:140582/149‑1 (MQ=255)
                                    cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                             <  6:54799/149‑1 (MQ=255)
                                     gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGtt                                                                                                                                          <  3:25200/119‑1 (MQ=255)
                                     gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGtt                                                                                                                                          >  4:25200/1‑119 (MQ=255)
                                     gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCggg                                                                                                            <  3:137011/149‑1 (MQ=255)
                                         ccAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCggg                                                                                                            >  3:265745/1‑145 (MQ=255)
                                         ccAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCggg                                                                                                            <  4:265745/145‑1 (MQ=255)
                                         ccAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAg                                                                                                        <  7:315773/149‑1 (MQ=255)
                                              agCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCcgc                                                                                                   <  3:4991/149‑1 (MQ=255)
                                                 ctgctgGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGaa                                                                                                <  4:230671/149‑1 (MQ=255)
                                                                                                 cGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACCGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAAc                                                                  <  5:304092/132‑1 (MQ=255)
                                                                                                                  ccATTCGCGGGCTTGCCCGCT‑CGCACAGGACCTGCACCAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCACGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                >  7:249916/1‑149 (MQ=14)
                                                                                                                     ttCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATca                              >  1:12685/1‑148 (MQ=255)
                                                                                                                       cgcgGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATcaca                            <  2:12685/148‑1 (MQ=255)
                                                                                                                                  ccGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                >  3:203771/1‑149 (MQ=255)
                                                                                                                                   cGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGTTGTTCg               >  6:275843/1‑149 (MQ=255)
                                                                                                                                    gCT‑CGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt              >  2:331572/1‑149 (MQ=18)
                                                                                                                                     cT‑CGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt             >  4:216163/1‑149 (MQ=18)
                                                                                                                                         ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTgg                     <  7:13485/138‑1 (MQ=255)
                                                                                                                                         ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTgg                     >  8:13485/1‑138 (MQ=255)
                                                                                                                                          cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                    >  2:32916/1‑138 (MQ=255)
                                                                                                                                          cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgt                    <  1:32916/138‑1 (MQ=255)
                                                                                                                                          cacaGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt         >  1:310981/1‑149 (MQ=255)
                                                                                                                                                cATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAgg   >  1:179097/1‑149 (MQ=255)
                                                                                                                                                 aTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAggg  >  8:204868/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCAC‑ACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCT‑CCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGG  >  NC_002947/2063050‑2063341

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