Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,377,518 C→T 9.3% G33D (GGC→GAC)  PP_3852 ← BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,377,5180CT 9.3% 106.7 / 2.6 43G33D (GGC→GAC) PP_3852BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base T (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GAATAAAGATAGTCTGAACGACCTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAACGCTGGGAGGCTGATG  >  NC_002947/4377380‑4377663
                                                                                                                                           |                                                                                                                                                 
gAATAAAGATAGTCTGAACGACCTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCg                                                                                                                                        >  1:1192059/1‑150 (MQ=255)
  aTAAAGATAGTCTGAACGACCTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCt                                                                                                                                      >  1:1088955/1‑150 (MQ=255)
                aaCGACCTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGt                                                                                                                        >  2:129867/1‑150 (MQ=255)
                   gACCTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTAc                                                                                                                     >  1:875257/1‑150 (MQ=255)
                     ccTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTTAGATTCCACATTCAAGACCCAGTGGTTGAGTTCTTTATTGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATTCCTCTGGTTCCGCTTTCTAGACCTTAGTTACAg                                                                                                                   >  2:535298/1‑150 (MQ=255)
                      cTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTTTTTTTCCATATTCAAAACCCATTGTTTTATTTTTTTTTCGAATTTCGCAAGCTCCTTTCCGCCCACTTTAACATCTGATGCCTTTGTTGCCGCTTTCGGAACCTTATTTTCAgg                                                                                                                  >  2:70070/1‑150 (MQ=255)
                          ggCCCCTGGGGAAATATCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAgt                                                                                                              <  2:544795/149‑1 (MQ=255)
                                gggggAAATACATGCTTTTTTTATTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCACTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAgtccgtc                                                                                                         <  2:544408/148‑1 (MQ=255)
                                          tCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACGGGCAGTCCGTCCATCGATTttc                                                                                              >  1:511235/1‑150 (MQ=255)
                                             tgctgtttt‑tttGAGCTTCCACATCCAAGCCCCAGTGGTGGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCACCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTttcttc                                                                                           <  1:94133/145‑1 (MQ=255)
                                             gCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTttcttc                                                                                           <  1:70070/150‑1 (MQ=255)
                                               ttttttt‑tttGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTttcttc                                                                                           <  1:1359494/148‑1 (MQ=255)
                                               ttttttt‑tttGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTttcttc                                                                                           >  2:1359494/1‑148 (MQ=255)
                                                   ttt‑tttGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCGCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTtcttcctcttc                                                                                     <  1:535298/150‑1 (MQ=255)
                                                     t‑tttGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACGGGCAGTCCGTCCATCGATTTTcttcctcttcct                                                                                   <  2:511235/150‑1 (MQ=255)
                                                              ttCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGGTCCCTTCCGCCCACGATAACAGCTAATGCCTCTGTTTCCGCTATCTAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATcgg                                                                           >  2:462794/1‑148 (MQ=255)
                                                               tCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCgg                                                                          <  2:1192059/150‑1 (MQ=255)
                                                                        aaGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCa                                                                 >  2:302055/1‑150 (MQ=255)
                                                                             ccAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATcc                                                                            <  2:539188/134‑1 (MQ=255)
                                                                             ccAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATcc                                                                            >  1:539188/1‑134 (MQ=255)
                                                                             ccAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTATTCTTGTTCGTGGTATCCTGGCATTCCCaaaaaa                                                            >  1:525251/1‑146 (MQ=255)
                                                                              cAGTGGTTGAGCTATGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATg                                                           >  2:874518/1‑150 (MQ=255)
                                                                                     tGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAAc                                                    <  1:302055/150‑1 (MQ=255)
                                                                                               aTCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCa                                                                       <  2:635086/121‑1 (MQ=255)
                                                                                               aTCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCa                                                                       >  1:635086/1‑121 (MQ=255)
                                                                                                  gAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAgtccgt                                                                                                          >  2:996770/1‑83 (MQ=255)
                                                                                                  gAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAgtccgt                                                                                                          <  1:996770/83‑1 (MQ=255)
                                                                                                    aTCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATg                                              >  2:753436/1‑141 (MQ=255)
                                                                                                    aTCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATg                                              <  1:753436/141‑1 (MQ=255)
                                                                                                              cTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGc                           <  1:874518/150‑1 (MQ=255)
                                                                                                               tCCCTTCCGCCCACGATAACAGCTGATGTCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAAcatcat                                               >  1:1295464/1‑129 (MQ=255)
                                                                                                               tCCCTTCCGCCCACGATAACAGCTGATGTCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAAcatcat                                               >  1:1295805/1‑129 (MQ=255)
                                                                                                               tCCCTTCCGCCCACGATAACAGCTGATGTCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAAcatcat                                               <  2:1295464/129‑1 (MQ=255)
                                                                                                               tCCCTTCCGCCCACGATAACAGCTGATGTCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAAcatcat                                               <  2:1295805/129‑1 (MQ=255)
                                                                                                                   ttCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGaaa                      >  2:262114/1‑150 (MQ=255)
                                                                                                                        cccACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAAc                 >  2:61198/1‑150 (MQ=255)
                                                                                                                         ccACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAACg                >  1:1202096/1‑150 (MQ=255)
                                                                                                                              aTAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCt                                         <  2:850337/120‑1 (MQ=255)
                                                                                                                              aTAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCt                                         >  1:850337/1‑120 (MQ=255)
                                                                                                                                  cAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAACGCTGGGAGGc       <  2:1088955/150‑1 (MQ=255)
                                                                                                                                     cTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAACGCTGGGAGGCtga    >  2:849070/1‑150 (MQ=255)
                                                                                                                                      tGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAACGCTGGGAGGCtgat   <  1:129867/150‑1 (MQ=255)
                                                                                                                                       gATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAACGCTGGGAGGCtgatg  >  2:1353827/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GAATAAAGATAGTCTGAACGACCTGCAGCCCCTTGGGAAATATCTGCTTTTTTT‑TTTGAGATTCCACATTCAAGACCCAGTGGTTGAGCTCTGTATCGAATCGCGCAAGCTCCCTTCCGCCCACGATAACAGCTGATGCCTCTGGTGCCGCTGTCGAGACCTGAGTTACAGGCAGTCCGTCCATCGATTTTCTTCCTCTTCCTGGTATCCGGGCATTTCCAAGGATGACGTAACATCATGCTCCTGCGTAATCTCATGCGGAAACTAACGCTGGGAGGCTGATG  >  NC_002947/4377380‑4377663

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