Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,798,380 C→A 12.4% D495Y (GAC→TAC) ‡ PP_5077 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,798,3800CA12.4% 135.6 / 8.3 48D495Y (GAC→TAC) ‡PP_5077hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base A (2/4);  total (20/28)
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.09e-01

CTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGC  >  NC_002947/5798238‑5798516
                                                                                                                                              |                                                                                                                                        
cTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGtt                                                                                                                                    >  3:172212/1‑149 (MQ=255)
  aCTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGTACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCGACGCAGCGGTTgc                                                                                                                                  <  1:238056/149‑2 (MQ=255)
   cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                                                                                                 >  6:311805/1‑149 (MQ=255)
   cTGGGCGGCGGGGGGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCGAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCATCCTGGAGCTTCGCTGGCAATTTCTTGATTGCCGCAACCGCGACGCAGCGGTTgcg                                                                                                                                 <  1:8893/149‑3 (MQ=255)
      ggcggcGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCGAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGGCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                                                                                                 <  4:72219/146‑1 (MQ=255)
         ggcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACt                                                                                                                           <  4:12471/149‑1 (MQ=255)
             gggTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCg                                                                                                                       <  4:172212/149‑1 (MQ=255)
                  aTGCCGCCCGGAAGGTATCAGCGGGCAGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACGAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGt                                                                                                                  <  7:289960/149‑1 (MQ=255)
                      cgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                                                                                                               >  8:256139/1‑148 (MQ=255)
                       gccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCacga                                                                                                             <  1:215066/149‑1 (MQ=255)
                       gccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCacga                                                                                                             <  6:20857/149‑1 (MQ=255)
                          cGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACTCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCTTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGt                                                                                                          <  2:189030/149‑1 (MQ=255)
                                  aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGACACTGCCGTAAGTCACGACGTACAGAGGc                                                                                                  >  2:228825/1‑149 (MQ=255)
                                  aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGc                                                                                                  >  4:275995/1‑149 (MQ=255)
                                     aGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTg                                                                                               <  5:311805/149‑1 (MQ=255)
                                                  cAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGtttt                                                                                  >  4:113931/1‑149 (MQ=255)
                                                          gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                             <  3:275995/146‑1 (MQ=255)
                                                          gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGa                                                                            <  7:256139/147‑1 (MQ=255)
                                                              gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                           <  8:184986/144‑1 (MQ=255)
                                                              gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                           >  7:184986/1‑144 (MQ=255)
                                                              gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGc                                                                      <  3:113931/149‑1 (MQ=255)
                                                                     ccGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCg                                                                     >  4:223518/1‑143 (MQ=255)
                                                                     ccGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCg                                                                     <  3:223518/143‑1 (MQ=255)
                                                                       gACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGc                                                             <  5:208650/149‑1 (MQ=255)
                                                                             aCGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCt                                                       >  6:260502/1‑149 (MQ=255)
                                                                              cGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTg                                                      <  5:191612/149‑1 (MQ=255)
                                                                                tGGCCAAGGTTTACCAGACGGGACGTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGTCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTggg                                                    <  6:127397/149‑1 (MQ=255)
                                                                                  gCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTc                                                               >  4:279973/1‑136 (MQ=255)
                                                                                  gCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTc                                                               <  3:279973/136‑1 (MQ=255)
                                                                                     aaGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCa                                                                                        >  1:32057/1‑108 (MQ=255)
                                                                                     aaGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCa                                                                                        <  2:32057/108‑1 (MQ=255)
                                                                                     aaGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                >  4:41061/1‑148 (MQ=255)
                                                                                      aGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCAGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGa                                              <  1:228825/149‑1 (MQ=255)
                                                                                          ttACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGt                                           >  7:302639/1‑148 (MQ=255)
                                                                                            aacAGCCTGGACGTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTa                                          <  4:81200/145‑1 (MQ=255)
                                                                                            aCCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTaaa                                        >  3:81200/1‑148 (MQ=255)
                                                                                             ccAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGc                                       <  8:302639/149‑1 (MQ=255)
                                                                                                               ggCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGCCGTAAGCCTGGGCCGACGCCTCg                       >  5:255476/1‑147 (MQ=255)
                                                                                                                   aGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTg                 <  5:260502/149‑1 (MQ=255)
                                                                                                                       cTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTa              <  3:41061/148‑1 (MQ=255)
                                                                                                                           aTTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAg         >  6:268537/1‑149 (MQ=255)
                                                                                                                               ccGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGc                                                             <  5:276816/93‑1 (MQ=255)
                                                                                                                               ccGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGc                                                             >  6:276816/1‑93 (MQ=255)
                                                                                                                                 gCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTg   >  2:48763/1‑149 (MQ=255)
                                                                                                                                  cAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTTCAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTg   >  1:225921/1‑148 (MQ=255)
                                                                                                                                  cAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCg                                                              >  4:316503/1‑89 (MQ=255)
                                                                                                                                  cAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCg                                                              <  3:316503/89‑1 (MQ=255)
                                                                                                                                  cAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGc  <  2:225921/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                        
CTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGC  >  NC_002947/5798238‑5798516

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