Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 353,067 A→G intergenic (+34/+84) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,0670AG97.6% 135.6 / ‑6.0 42intergenic (+34/+84)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/24);  minor base C (0/1);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352926‑353208
                                                                                                                                             |                                                                                                                                               
gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGc                                                                                                                                          >  7:290685/1‑149 (MQ=255)
 aaGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGcc                                                                                                                                         >  5:345916/1‑149 (MQ=255)
       aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCCAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCCTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTcc                                                                                                                                   <  5:336913/149‑1 (MQ=255)
          ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTAGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCcac                                                                                                                                <  3:81347/149‑1 (MQ=255)
          ccAGCGCGGGGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGTCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                                 <  2:253674/148‑1 (MQ=255)
            aGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAg                                                                                                                              <  3:325686/149‑1 (MQ=255)
                cgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                                              >  8:133636/1‑129 (MQ=255)
                cgGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGc                                                                                                                                              <  7:133636/129‑1 (MQ=255)
                    gCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTg                                                                                                                      <  7:324686/149‑1 (MQ=255)
                        gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑At                                                                                                                 <  7:189719/149‑1 (MQ=255)
                            ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCTCTCCCACAGGTTTGCTGAG‑ATctct                                                                                                             <  5:7754/149‑1 (MQ=255)
                             cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                                        >  1:47443/1‑122 (MQ=255)
                             cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGccc                                                                                                                                        <  2:47443/122‑1 (MQ=255)
                             cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCtg                                                                                                            >  6:118729/1‑149 (MQ=255)
                                  tCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTgg                                                                                                         >  7:68616/1‑147 (MQ=255)
                                    ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                       >  8:187344/1‑132 (MQ=255)
                                    ttCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                       <  7:187344/132‑1 (MQ=255)
                                       cACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTCTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                                       >  2:325552/1‑129 (MQ=255)
                                        aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGc                                                                                                <  8:68616/149‑1 (MQ=255)
                                                   tACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCcgcg                                                                                       >  1:310653/1‑147 (MQ=255)
                                                          tGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCcgc                                                                              <  3:284877/149‑1 (MQ=255)
                                                            ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                                    <  3:376813/111‑1 (MQ=255)
                                                            ccGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                                    >  4:376813/1‑111 (MQ=255)
                                                             cGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGAGAT‑CTCAGCAACCCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCTGCaaa                                                                           <  5:181787/149‑1 (MQ=17)
                                                                  gCCAAGGCCTACATGGCTGCGCGCGGGATCGTCGTTCGCTGAGCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGc                                                                      <  8:81130/149‑1 (MQ=255)
                                                                    cAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                    >  5:165829/1‑149 (MQ=255)
                                                                         ccTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCTTTGACCGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAAtt                                                               <  5:162666/149‑1 (MQ=255)
                                                                               tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGCG‑CGGGCGTGCCCGCGCAGGGCCGCGAAGCGGccc                                                                    >  2:51413/1‑138 (MQ=255)
                                                                               tGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAGCCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGccc                                                                    <  1:51413/139‑1 (MQ=255)
                                                                                ggCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTTCTGTG‑ATGTGTGTTGGCT‑CGGGCGGTCCCCCGa                                                                                      >  7:344770/1‑119 (MQ=21)
                                                                                       cgcgGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGGTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCCCGAAGCGGCCCCAATTGCTTAAACCGTGGt                                                 >  7:52366/1‑149 (MQ=255)
                                                                                            gATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTg                                            <  2:310653/149‑1 (MQ=255)
                                                                                             aTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGTAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                           <  6:345916/149‑1 (MQ=255)
                                                                                                   tttCGCTTATCCATTGCCGGGGGCTGCTTTGCAGCCCTTCGGGGCCACGCCCGCTCCCACTGGTTTGCTCAT‑ATCTCTGGGGTTG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagca                                     <  8:52366/148‑1 (MQ=37)
                                                                                                           aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTgg                             <  5:118729/149‑1 (MQ=255)
                                                                                                           aTCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTg                              <  6:165829/148‑1 (MQ=255)
                                                                                                                  gACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCt                      >  4:252522/1‑149 (MQ=255)
                                                                                                                       gggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTa                 <  3:252522/149‑1 (MQ=255)
                                                                                                                             ctTTGCAGCCCTTCGCGGGCACGCCCCCACCCACACGGTTGGTTAG‑ATGTCTGTGGGAG‑CCGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTTGTTCAAACGTGGTTCATTGCGAGCAGCGCGTGGTAAAGCTCGGGAtacac            >  6:140742/1‑144 (MQ=21)
                                                                                                                                      ccTTCGCGGGCACGCCCGGTCCCACAGGTTTGCTGAG‑ATCTATTTGGTAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATcccc  <  5:140742/149‑1 (MQ=39)
                                                                                                                                       cTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa          <  7:73746/140‑1 (MQ=255)
                                                                                                                                       cTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGa          >  8:73746/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAG‑ATCTCTGTGGGAG‑CGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAGTATCCCC  >  NC_002947/352926‑353208

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