Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,843,918 C→A 8.3% L241M (CTG→ATG) ‡ PP_4259 → iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,843,9180CA 8.3% 142.2 / 2.4 48L241M (CTG→ATG) ‡PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base A (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCT  >  NC_002947/4843777‑4844066
                                                                                                                                             |                                                                                                                                                    
cAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgtg                                                                                                                                              >  1:652865/1‑150 (MQ=255)
 aGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCgtgt                                                                                                                                             >  1:591871/1‑150 (MQ=255)
   gCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTcc                                                                                                                                           >  2:303603/1‑150 (MQ=255)
    ccGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCt                                                                                                                                          <  2:867270/150‑1 (MQ=255)
      gATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCttg                                                                                                                                        <  1:303603/150‑3 (MQ=255)
            tGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCa                                                                                                                                  >  1:1060130/1‑150 (MQ=255)
             gggCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCAc                                                                                                                                 >  1:937429/1‑150 (MQ=255)
                 cTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCg                                                                                                                             >  1:848323/1‑150 (MQ=255)
                     ttCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGc                                                                                                                         <  1:442737/150‑1 (MQ=255)
                            tCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAAccc                                                                                                                  >  1:242892/1‑150 (MQ=255)
                               gCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCcgcg                                                                                                               <  2:389194/150‑1 (MQ=255)
                                 ttcttcGCCCTCGCCACTTACGGCAACGCCGGCTGGCTAAGCGAACGGGTGTGCGTGTACATGTGCCCCTATGTCCGCTTCCAGAGCGTGATTTTGGACAAGGACAAGCTGTTCGTGTCGTACGACCCACGCCGTGGCGAACCCCGCGGc                                                                                                             >  1:416758/1‑150 (MQ=255)
                                 ttcttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGc                                                                                                             >  2:404372/1‑150 (MQ=255)
                                  tcttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGCGTCCTACGACCCACGCCGTGGCGAAACCCGCGGcc                                                                                                            <  2:591871/150‑1 (MQ=255)
                                          cTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaag                                                                                                    <  2:1060130/150‑1 (MQ=255)
                                             gCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGAGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaagaaa                                                                                                 <  1:1092577/150‑1 (MQ=255)
                                                aCTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGAt                                                                                              <  1:404372/150‑1 (MQ=255)
                                                aCTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGAt                                                                                              >  1:327106/1‑150 (MQ=255)
                                                  ttACGGCAACGCCGGCTGGCTCCGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTCACAACGGTCAAGGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATAc                                                                                            <  1:748115/150‑1 (MQ=255)
                                                   tACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATAcc                                                                                           <  2:191483/150‑1 (MQ=255)
                                                                   ggctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAgg                                                                           <  2:848323/150‑1 (MQ=255)
                                                                     ctACGCGAACAGGTGTGCGTGTACATGTGCCCCTAGTCTCGCTTGCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGAGGCCCGCGCAAGAGAGATACAGAGTACAAGGCCCAGGGc                                                                         >  1:420916/1‑150 (MQ=37)
                                                                       aCGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCt                                                                       <  2:327106/150‑1 (MQ=255)
                                                                           gAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGt                                                                   <  2:937429/150‑1 (MQ=255)
                                                                            aaCAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATa                                                                                             >  2:602821/1‑123 (MQ=255)
                                                                            aaCAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATa                                                                                             <  1:602821/123‑1 (MQ=255)
                                                                                   gtgCGTGTACATGTGCCCCTATGCCCGCTCCTAGTGCATTTTGTTCAACAAGGACTCGCTGATCGTGTCCTACGACCCGCGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACCAGGCCCAGGGCCTGGGTGACTGCAt                                                           <  2:420916/150‑1 (MQ=255)
                                                                                      cGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGa                                                        >  1:769918/1‑150 (MQ=255)
                                                                                      cGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGa                                                        <  1:565132/150‑1 (MQ=255)
                                                                                                     cTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgcgt                                         >  1:994388/1‑150 (MQ=255)
                                                                                                         gCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGTACGCGCAAGAAAGATACCATTTTTGATCCCCATTGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                     <  1:748944/150‑1 (MQ=37)
                                                                                                          cccGCTTCCAGAGCGTGATGGTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACGGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAgg                                    >  1:562022/1‑150 (MQ=255)
                                                                                                           ccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGt                                   >  2:593577/1‑150 (MQ=255)
                                                                                                           ccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGt                                   >  2:484794/1‑150 (MQ=255)
                                                                                                                    gagCGTGATGTTCGACAAGGACACGCTGATCGTGTACTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGGCTGCACCATGTGCGTGCAGGTCTGCCCGAc                          >  2:594310/1‑150 (MQ=255)
                                                                                                                          gATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCGCGCGGAGACCAAACCGCCGTCCCGCTCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGTTCTGCCCGACCGGCAt                    <  2:242892/150‑1 (MQ=38)
                                                                                                                          gATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                     >  1:167121/1‑101 (MQ=255)
                                                                                                                          gATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                     <  2:167121/101‑1 (MQ=255)
                                                                                                                          gATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCAt                    <  1:594620/150‑1 (MQ=255)
                                                                                                                            tGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCg                  <  2:769918/150‑1 (MQ=255)
                                                                                                                                     aGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTg         >  1:874095/1‑150 (MQ=255)
                                                                                                                                      ggACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGAc                                                       <  2:762428/103‑1 (MQ=255)
                                                                                                                                      ggACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGAc                                                       >  1:762428/1‑103 (MQ=255)
                                                                                                                                        acacGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACACGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACg      <  2:689708/150‑1 (MQ=255)
                                                                                                                                            gCTGATCGTGTCCTACGACCCACGGCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGACCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCt  >  2:428133/1‑150 (MQ=255)
                                                                                                                                            gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGc                               <  1:690592/121‑1 (MQ=255)
                                                                                                                                            gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGc                               >  2:690592/1‑121 (MQ=255)
                                                                                                                                            gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCt  >  2:373869/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCT  >  NC_002947/4843777‑4844066

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑