Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,896,928 C→T 8.7% R328H (CGT→CAT)  PP_3440 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,896,9280CT 8.7% 127.2 / 2.4 47R328H (CGT→CAT) PP_3440hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/17);  new base T (2/2);  total (28/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 = 2.24e-01

ACTACACATTTATCGAGGGAACAGATCCAACGTGCGCTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACATCGGACTTCAGCGGGTTGA  >  NC_002947/3896783‑3897075
                                                                                                                                                  |                                                                                                                                                   
actacACATTTATCGAGGCAACAGATCCAACGTGCGCTCCACCTCGTCAATATCGAGCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAACGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACgg                                                                                                                                                   <  2:274827/148‑1 (MQ=255)
           aTCGAGGGAACAGATCCAACGTGCGCTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTg                                                                                                                                        >  3:106516/1‑148 (MQ=255)
                   aaCAGATACAACGTGCGCTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTcag                                                                                                                               >  3:183336/1‑149 (MQ=255)
                      aGATCCAACGTGCGCTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTcagcag                                                                                                                            <  5:181432/149‑1 (MQ=255)
                        aTCCAACGTGCGCTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAATGGTTCGGCGGCTGCAAAGTcagc                                                                                                                              <  1:288494/145‑1 (MQ=255)
                        aTCCAACGTGCGCTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAATGGTTCGGCGGCTGCAAAGTcagc                                                                                                                              >  2:288494/1‑145 (MQ=255)
                                  cgcTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGc                                                                                                                <  4:229670/149‑1 (MQ=255)
                                  cgcTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGc                                                                                                                <  4:106516/149‑1 (MQ=255)
                                      ccACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTccc                                                                                                                  >  5:277396/1‑143 (MQ=255)
                                      ccACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTccc                                                                                                                  <  6:277396/143‑1 (MQ=255)
                                         ccTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTgcag                                                                                                          >  8:320507/1‑148 (MQ=255)
                                          cTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAATGGTTcggcg                                                                                                                                            >  5:132239/1‑113 (MQ=255)
                                          cTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAATGGTTcggcg                                                                                                                                            <  6:132239/113‑1 (MQ=255)
                                           tCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACCGACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTgcagca                                                                                                        >  1:77030/1‑149 (MQ=255)
                                                aaTATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCa                                                                                                  <  2:123678/149‑1 (MQ=255)
                                                    tCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGc                                                                                               >  7:51282/1‑148 (MQ=255)
                                                       aTCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATg                                                                                           >  1:362505/1‑149 (MQ=255)
                                                            taaaCCATCGCATTCCGGTAGGCCTACA‑AAAAAGCCGAAATGCTGGTGATAACGATCCGTCAGGTATTTGTAGTAAATGACGAAAAGGTTCGGCGGCTGCAAAGTCATAAGACAGCCTCCCGACTCAAGCACATTAACGCCATGCAcc                                                                                       >  8:329440/2‑148 (MQ=255)
                                                                 cATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTgg                                                                                 >  8:307153/1‑149 (MQ=255)
                                                                  aTCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCCGCACCTTCACGCCATGCACCACCTGGc                                                                                >  8:273749/1‑149 (MQ=255)
                                                                      ccTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGATGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTAcc                                                                            >  4:107600/1‑149 (MQ=255)
                                                                           cGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTcc                                                                        <  7:320507/148‑1 (MQ=255)
                                                                                   cTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGAt                                                                >  5:62273/1‑148 (MQ=255)
                                                                                       c‑acaAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTAAGCAGACAGCCTACCGAATGCAGCACATTCACGCCATGCACCAACTGGCTACCCTACACCCCGATGaac                                                            >  3:300202/1‑146 (MQ=255)
                                                                                         acaAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCCCCCCGATGACCAc                                                          >  2:112918/1‑149 (MQ=255)
                                                                                         acaAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCAc                                                          >  2:108356/1‑149 (MQ=255)
                                                                                                 gAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGccc                                                  >  7:288035/1‑149 (MQ=255)
                                                                                                 gAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGccc                                                  <  2:362505/149‑1 (MQ=255)
                                                                                                          gTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCa                                           <  8:51282/147‑1 (MQ=255)
                                                                                                          gTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAgg                                         >  1:84586/1‑149 (MQ=255)
                                                                                                            gATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGcc                                       >  1:336165/1‑149 (MQ=255)
                                                                                                            gATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGcc                                       >  3:297197/1‑149 (MQ=255)
                                                                                                             aTCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCcg                                      >  7:260378/1‑149 (MQ=255)
                                                                                                              tCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCcgc                                     >  4:177614/1‑149 (MQ=255)
                                                                                                                aCGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGa                                   <  8:288035/149‑1 (MQ=255)
                                                                                                                        tCAGGTATTTGTATTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTTGCTACCCTCCACCCCGATGTCCACTTTAGCCCCGGTGCAGGCCGCGACAATACTc                           <  7:329440/149‑1 (MQ=255)
                                                                                                                        tCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTc                           <  3:107600/149‑1 (MQ=255)
                                                                                                                          aGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCg                          <  7:273749/148‑1 (MQ=255)
                                                                                                                            gTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGca                       >  4:301664/1‑149 (MQ=255)
                                                                                                                             tATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACACTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAAAACTCGGca                       >  4:254583/1‑148 (MQ=255)
                                                                                                                               tttGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACAt                    >  3:514/1‑149 (MQ=255)
                                                                                                                                   tagtaACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACATCGGa                <  2:84586/149‑1 (MQ=255)
                                                                                                                                     gtaACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACATCGGACt              >  4:230408/1‑149 (MQ=255)
                                                                                                                                       aaCTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACATCGGACtt             <  4:297197/148‑1 (MQ=255)
                                                                                                                                             cGAAACGGTTCGGCGTCTGCAAAGTCAGCAGTCAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCTCAGGCCGCGACAATACTCGGCACATCGGACTTCAGCgg       <  4:300202/148‑1 (MQ=255)
                                                                                                                                                aaCGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACATCGGACTTCAGCGGGtt    >  3:365382/1‑148 (MQ=255)
                                                                                                                                                 aCGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACATCGGACTTCAGCGGGTTGa  <  8:260378/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACTACACATTTATCGAGGGAACAGATCCAACGTGCGCTCCACCTCGTCAATATCGATCCGGAAACCATCGCCTTCCGGTAGGCCTACC‑ACAAAGCCGAAATGCTGGTGATCACGATCCGTCAGGTATTTGTAGTAACTGACGAAACGGTTCGGCGGCTGCAAAGTCAGCAGACAGCCTCCCGCCTGCAGCACATTCACGCCATGCACCAACTGGCTACCCTCCACCCCGATGACCACTTTAGCCCCGGCGCAGGCCGCGACAATACTCGGCACATCGGACTTCAGCGGGTTGA  >  NC_002947/3896783‑3897075

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