Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,743,385 C→T 13.3% R161C (CGC→TGC)  PP_2397 → EF hand domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,743,3850CT13.3% 71.1 / 3.0 30R161C (CGC→TGC) PP_2397EF hand domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (15/11);  new base T (2/2);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

GCAAGGCTTCCCAAAGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCAGGCCAATGGCAAAGCCATGCCGGAT  >  NC_002947/2743239‑2743528
                                                                                                                                                   |                                                                                                                                               
gCAAGGCTTCCCAAAGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTcgcgc                                                                                                                                               >  3:99217/1‑149 (MQ=255)
      tgatccAAAGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAg                                                                                                                                         <  6:165641/145‑1 (MQ=255)
       ttCCCAAAGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTcgcgcg                                                                                                                                              >  8:43412/1‑143 (MQ=255)
       ttCCCAAAGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTcgcgcg                                                                                                                                              <  7:43412/143‑1 (MQ=255)
          ccAAAGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTc                                                                                                                                     >  4:332586/1‑149 (MQ=255)
            aaaGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTg                                                                                                                                   <  4:99217/149‑1 (MQ=255)
                     gCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGcc                                                                                                                          <  6:125008/149‑1 (MQ=255)
                       atgTGGGGTTAGGTTATTTGGGGCGCCGGGTCCAAGTCCGGGGCTGA‑GGGCATGTGTGGTGAATGCTAGTGTTGTGACGCCTGGTTCGGCGAATCCGACAGCATTCAGGCTGGAAAGGTCTCGCGCGATGAGTTTCTGGGGGTGGCCaa                                                                                                                        <  4:37835/147‑1 (MQ=255)
                           gCGGTGACGGTAACTCCGGCGCCGGCGCCAAGGGCGGGGCGGA‑GGTCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGAc                                                                                                                    <  8:223180/149‑1 (MQ=255)
                                    gTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGAGGGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCgg                                                                                                            >  1:63587/1‑149 (MQ=255)
                                    gTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGt                                                                                                           <  3:332586/149‑1 (MQ=255)
                                                                   tGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACg                                                                            >  8:266118/1‑149 (MQ=255)
                                                                          cAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCa                                                                      >  4:87243/1‑149 (MQ=255)
                                                                                  gtgAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGa                                                               <  2:63587/148‑1 (MQ=255)
                                                                                   tgtatgCAAGTGTGGGGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGGCGCGCGCGATGAGTTTCTGGCGGTGGGCAAGGACCGTGCCGGTGATTTCGTCACGATCGATCGCAACCATGACGGGTTCATCTCCGAGc                                                             <  5:195642/144‑1 (MQ=255)
                                                                                        gCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGGGATTTCGACACGATCGATCGCAACCCTGACGGTTTCATCTCCGAGCAgg                                                          >  8:245438/1‑147 (MQ=255)
                                                                                           aGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCg                                                     >  8:63816/1‑149 (MQ=255)
                                                                                                   gACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGTGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGc                                                   <  3:326543/143‑1 (MQ=255)
                                                                                                   gACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGTGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGc                                                   <  7:188882/143‑1 (MQ=255)
                                                                                                   gACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGTGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGc                                                   >  8:188882/1‑143 (MQ=255)
                                                                                                   gACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGTGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGc                                                   >  4:326543/1‑143 (MQ=255)
                                                                                                      gCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCt                                          >  8:150535/1‑149 (MQ=255)
                                                                                                      gCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGACGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGc                                                   >  2:35877/1‑140 (MQ=255)
                                                                                                      gCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGACGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGc                                                   <  1:35877/140‑1 (MQ=255)
                                                                                                                ccGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGAtt                                >  8:189199/1‑149 (MQ=255)
                                                                                                                     aCCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCa                           <  7:266118/149‑1 (MQ=255)
                                                                                                                        gACAGCAATCACGATAGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCgg                                                    >  2:293436/1‑121 (MQ=255)
                                                                                                                             cAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCAGGCCAATg                   >  4:15767/1‑149 (MQ=255)
                                                                                                                              aaTCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCAGGCCAATgg                  >  2:327749/1‑149 (MQ=255)
                                                                                                                                   cGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCAGGCCAATGGCAAAg             >  4:116868/1‑149 (MQ=255)
                                                                                                                                          aaGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCAGGCCAATGGCAAAGCCATGcc      >  4:172228/1‑149 (MQ=255)
                                                                                                                                          aaGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGACCGATCACCTCAACGCGATTTACCAGGCCAATGGAAAAGCCATg        >  6:246143/1‑147 (MQ=255)
                                                                                                                                              tctcGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCAGGCCAATGGCAAAGCCATGCCGGAt  <  7:150535/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCAAGGCTTCCCAAAGCGAAGGCAAGTGCGGTGAGGGTAAATGCGGCGCCGGCGCCAAGGCCGGGGCTGA‑GGGCAAGTGTGGTGAAGGCAAGTGTGGTGACGCCTCGTTCGCCGAAACCGACAGCAATCACGATGGCAAGGTCTCGCGCGATGAGTTTCTGGCGGTGGCCAAGGACCGTGCCGGTGATTTCGACACGATCGATCGCAACCATGACGGTTTCATCTCCGAGCAGGAAGCGGCCGATCACCTCAAGGCGATTTACCAGGCCAATGGCAAAGCCATGCCGGAT  >  NC_002947/2743239‑2743528

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑