Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,586,057 +C intergenic (+140/+75) PP_4061 → / ← PP_4063 hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,586,0561.C87.1% 91.2 / 8.2 31intergenic (+139/+76)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (10/17);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 3.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTT  >  NC_002947/4585910‑4586197
                                                                                                                                                       |                                                                                                                                             
aGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGa                                                                                                                                                <  5:380404/148‑1 (MQ=255)
 ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAcc                                                                                                                                             <  8:115060/149‑2 (MQ=255)
        gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATCTCCGGGTGCTTGGCTAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCgg                                                                                                                                       <  1:59157/149‑1 (MQ=38)
         ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGc                                                                                                                                      <  2:176916/149‑1 (MQ=39)
            aGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGACGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCt                                                                                                                                    <  6:383279/148‑1 (MQ=38)
                    gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcgg                                                                                                                           <  2:230148/149‑1 (MQ=39)
                     cACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGAGAAGGTGAAGATCTTGCATTGATCGGGCCCGCCTCTTCGCGGGTAAACCCGCTCATACCAGGGACCCTGCAGGCCTTGAATgcggc                                                                                                                          >  2:156937/1‑149 (MQ=255)
                          gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGc                                                                                                                     <  8:337655/149‑1 (MQ=39)
                              gAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGTGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCAATTAGCGGGTAAAACCCCTCCAAACAGGGGCCGGGCCGGCCGTGAATGGGGGGCGGcccc                                                                                                                  >  1:277725/1‑149 (MQ=255)
                                 cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGCCCGTGCCGGCCGTGAATTCGGCGGTGCCCCggg                                                                                                               >  6:300217/1‑148 (MQ=37)
                                 cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg                                                                                                               <  4:101802/148‑1 (MQ=39)
                                 cAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                                                                                                              <  1:156937/149‑1 (MQ=39)
                                 cAGCTGATCTGCGCTGCTCCATTGAAACGTTGAGGATGACCGGGTGCGTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                                                                                                              <  2:108801/149‑1 (MQ=38)
                                  aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                                                                                                              <  6:100147/148‑1 (MQ=39)
                                    cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGaa                                                                                                                                >  8:263607/1‑128 (MQ=38)
                                    cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGaa                                                                                                                                <  7:263607/128‑1 (MQ=38)
                                    cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAgg                                                                                                           <  2:75219/149‑1 (MQ=39)
                                    cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAgg                                                                                                           <  8:335719/149‑1 (MQ=39)
                                       aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGaa                                                                                                                                <  2:283993/125‑1 (MQ=38)
                                       aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGaa                                                                                                                                >  1:283993/1‑125 (MQ=38)
                                                  tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCACTTAGCGGGTAAAACCCCAACAACCAGGGGCCCGGCCGGCCATGAATTCGGCGCGGGCCCGGGTAGGCGCGCGTTTTCCcg                                                                                             >  8:79076/1‑149 (MQ=255)
                                                            ggTTGAGGATGACCGGGTGCTTGGCATGGTGAATATCTTGCATTGAGGGGGCCGGCCTCTTCGCGGGTAACCCCGCGCCTACCAGGGACCGGGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGGGGGTTTACCCGCGAAGAggcc                                                                                   <  2:277725/149‑1 (MQ=37)
                                                                         ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg                                                                                                      >  7:145915/1‑118 (MQ=37)
                                                                         ccGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg                                                                                                      <  8:145915/118‑1 (MQ=37)
                                                                             gTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg                                                                                                               >  4:325/1‑105 (MQ=37)
                                                                              tGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt                                                                                                              <  3:325/105‑1 (MQ=37)
                                                                                                                 cggccTCTTCGCGGGTAAACCCGCTCCTAACAGGGGCCCGGGCGGCCTTGAATGCGGCGCTGCCCCGGGGAGGTGCGGGTTTACCCGCGGATCGGCCGGCCGCGCCCGGCCAACACCCCGCCGATGCAGCCGCAATCTGCGCCCCCCTg                               >  2:7588/1‑149 (MQ=255)
                                                                                                                                           cTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCa      >  3:273565/1‑148 (MQ=255)
                                                                                                                                            tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc    >  7:380769/1‑149 (MQ=255)
                                                                                                                                             aCCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGa                              >  6:191413/1‑122 (MQ=255)
                                                                                                                                              ccAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGAt                             <  5:191413/122‑1 (MQ=255)
                                                                                                                                              ccAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGTATGCGGAACTTCTGCACCtt  >  2:102532/1‑149 (MQ=255)
                                                                                                                                              ccAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCtt  >  7:9311/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
AGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTT  >  NC_002947/4585910‑4586197

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: