Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,491,518 C→T 9.3% D301N (GAT→AAT)  PP_3096 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,491,5180CT 9.3% 115.8 / 2.7 43D301N (GAT→AAT) PP_3096hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/22);  new base T (2/2);  total (19/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 = 5.61e-01

TGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGAC  >  NC_002947/3491376‑3491653
                                                                                                                                              |                                                                                                                                       
tgGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTc                                                                                                                                   <  5:108625/149‑1 (MQ=255)
       cGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCt                                                                                                                            <  4:149547/149‑1 (MQ=255)
         gcgcAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCtat                                                                                                                          <  7:304132/149‑1 (MQ=255)
          cgcAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCtata                                                                                                                         <  6:339685/149‑1 (MQ=255)
                   aaGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAgcg                                                                                                                 >  3:299873/1‑148 (MQ=255)
                    aGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCa                                                                                                               <  1:356408/149‑1 (MQ=255)
                         ccGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAAc                                                                                                          >  8:59876/1‑149 (MQ=255)
                               cACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGAc                                                                                                     >  2:177661/1‑148 (MQ=255)
                                aCGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGAc                                                                                                     >  8:342214/1‑147 (MQ=255)
                                      aGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGc                                                                                             <  4:153178/149‑1 (MQ=255)
                                        ggCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCgg                                                                                           >  2:243137/1‑149 (MQ=255)
                                           gTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCGCGCGCGTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAg                                                                                        <  4:183127/149‑1 (MQ=255)
                                            tttGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGc                                                                                       <  1:226085/149‑1 (MQ=255)
                                                aGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGc                                                                                   >  5:374698/1‑149 (MQ=255)
                                                 gCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATTCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCg                                                                                  >  7:131215/1‑149 (MQ=255)
                                                 accAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATTCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCg                                                                                  >  7:131216/2‑149 (MQ=255)
                                                   cAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGc                                                                                       >  3:216260/1‑142 (MQ=255)
                                                   cAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGc                                                                                       <  4:216260/142‑1 (MQ=255)
                                                   cAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGt                                                                                >  2:136732/1‑149 (MQ=255)
                                                       ccagccagCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTc                                                                                                        <  8:65513/121‑1 (MQ=255)
                                                               cccAGGGCGACATTGCCGCTGCGGCCCAGTTCGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGaa                                                                    <  3:82635/149‑1 (MQ=255)
                                                                    ggCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCGCCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCt                                                               >  2:101561/1‑149 (MQ=255)
                                                                      cGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTgg                                                             <  1:177661/149‑1 (MQ=255)
                                                                      cGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTgg                                                             <  4:187437/149‑1 (MQ=255)
                                                                        tcATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTa                                                           <  1:319808/148‑1 (MQ=255)
                                                                          aTTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTa                                                           >  2:319808/1‑147 (MQ=255)
                                                                                 cTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAg                                                  >  4:353671/1‑149 (MQ=255)
                                                                                      gCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCGGGTAGTCGGCAAGGCt                                               >  7:152732/1‑147 (MQ=255)
                                                                                          aGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCgg                                         <  4:299873/149‑1 (MQ=255)
                                                                                           gTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCggg                                        <  2:115554/149‑1 (MQ=255)
                                                                                                  gcgTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATTCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCt                                 <  8:131215/149‑1 (MQ=255)
                                                                                                  gcgTGGTACCGCCTCGCGCGTGAGGCGCGGGGCAATTTCGAAATTCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCt                                 <  8:131216/149‑1 (MQ=255)
                                                                                                     tGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGa                              <  1:136732/149‑1 (MQ=255)
                                                                                                        tACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAATTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGc                           >  1:238938/1‑149 (MQ=255)
                                                                                                           cGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACAAGCTTGCGGTAGCAGGCGGTGCCCGGCTGTAAGCGGTGGTAGTCGGAAAGGCTGCGCGGGCCCAGCTCGATGCGGa                        >  7:175800/1‑149 (MQ=255)
                                                                                                                 gcgcTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCAc                   <  7:342214/148‑1 (MQ=255)
                                                                                                                 gcgcTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCAc                   <  7:59876/148‑1 (MQ=255)
                                                                                                                 gcgcTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCAc                   <  3:353671/148‑1 (MQ=255)
                                                                                                                 gcgcTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCAc                   <  6:374698/148‑1 (MQ=255)
                                                                                                                   gcTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGAc                >  5:269540/1‑149 (MQ=255)
                                                                                                                     ttGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTg              >  7:274569/1‑149 (MQ=255)
                                                                                                                          gcgcGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGt         <  8:333484/149‑1 (MQ=255)
                                                                                                                                 gCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGAc  >  7:328256/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                       
TGGAAGACGGCGCAACGCGAAGGCTCCGCGTCACGACCAGGGCGTTTGAGCCAACCCAGCCAGCCCAGGGCGACATTGCCGCTGCGGCCCAGTTGGGCGCGTGGTACCGCCTCGCGCTTGAGGCGCGGGGCAATTTCGAAATCCAGTTCGGCACCTATATAGAAGCGCACCAACTCGACCAGCTTGCGGTAGCAGGCGGTGCCCGGCTGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGAC  >  NC_002947/3491376‑3491653

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