Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,649,986 T→A 10.0% K310* (AAG→TAG) ‡ PP_4113 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,649,9860TA10.0% 133.8 / 6.1 50K310* (AAG→TAG) ‡PP_4113hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base A (2/3);  total (22/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 = 6.98e-01

ATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGG  >  NC_002947/4649840‑4650131
                                                                                                                                                  |                                                                                                                                                 
aTCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGc                                                                                                                                                 >  3:190207/1‑149 (MQ=255)
        cgcCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGcccc                                                                                                                                         <  1:153030/149‑1 (MQ=255)
         gcCGACCCGGTGAATGCCCCAGCTCCTGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCGGGGCTGGGCAAGCACAGccctg                                                                                                                                        >  3:257517/1‑147 (MQ=255)
           cGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCCCTGGc                                                                                                                                      <  4:257517/149‑1 (MQ=255)
                 ggTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCCCTGGcttcgag                                                                                                                                <  4:280989/149‑7 (MQ=255)
                 ggTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAg                                                                                                                                >  8:83972/1‑149 (MQ=255)
                 ggTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAg                                                                                                                                <  2:124381/149‑1 (MQ=255)
                 ggTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAg                                                                                                                                <  7:274890/149‑1 (MQ=255)
                 ggTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAg                                                                                                                                <  7:279435/149‑1 (MQ=255)
                   tGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGcc                                                                                                                              <  4:190207/149‑1 (MQ=255)
                                tCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGcc                                                                                                                                                <  2:22473/118‑1 (MQ=255)
                                tCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGcc                                                                                                                                                >  1:22473/1‑118 (MQ=255)
                                tCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGccc                                                                                                                             <  6:104310/137‑1 (MQ=255)
                                tCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGccc                                                                                                                             >  5:104310/1‑137 (MQ=255)
                                 ccAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAg                                                                                                                >  8:129389/1‑149 (MQ=255)
                                  cAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGc                                                                                                               >  3:138657/1‑149 (MQ=255)
                                   aGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGc                                                                                                               >  7:263719/1‑148 (MQ=255)
                                           aCCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGcc                                                                                                      >  2:216249/1‑149 (MQ=255)
                                             ccGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                    >  1:221741/1‑149 (MQ=255)
                                              cGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGc                                                                                                   <  5:61880/149‑1 (MQ=255)
                                                cACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGa                                                                                                 <  8:263719/149‑1 (MQ=255)
                                                         ttCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGc                                                                                        >  1:294464/1‑149 (MQ=255)
                                                           cGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCt                                                                                          <  5:277590/145‑1 (MQ=255)
                                                            gACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGt                                                                                     <  7:83972/149‑1 (MQ=255)
                                                                  gTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTg                                                                               >  4:12396/1‑149 (MQ=255)
                                                                     cACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCg                                                                            >  1:326444/1‑149 (MQ=255)
                                                                       ccGAGGCAGCGATGTTGTCGGTGATACCGGCATCCGAGTTGCGCGCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCggg                                                                          <  7:144155/149‑1 (MQ=255)
                                                                                        gCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGc                                                         <  1:210173/149‑1 (MQ=255)
                                                                                           ttGATACCGTCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACGCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTa                                                      <  5:270733/148‑1 (MQ=255)
                                                                                                aCCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAgg                                                  >  5:78989/1‑148 (MQ=255)
                                                                                                  cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                               >  2:90540/1‑149 (MQ=255)
                                                                                                  cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                               <  1:216249/149‑1 (MQ=255)
                                                                                                      aTCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCt                                           >  2:37508/1‑149 (MQ=255)
                                                                                                      aTCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCt                                           >  5:299133/1‑149 (MQ=255)
                                                                                                         cGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTgg                                         <  2:294464/148‑1 (MQ=255)
                                                                                                          gAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGcc                                       <  8:59099/149‑1 (MQ=255)
                                                                                                          gAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGcc                                       >  2:190153/1‑149 (MQ=255)
                                                                                                            gTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCtt                                     <  4:138657/149‑1 (MQ=255)
                                                                                                                     tACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                              <  3:12396/147‑1 (MQ=255)
                                                                                                                        cATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCTGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTg                                          <  2:96183/132‑1 (MQ=255)
                                                                                                                        cATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCAGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTg                                          >  1:96183/1‑132 (MQ=255)
                                                                                                                                  aaGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATg               <  6:78989/149‑1 (MQ=255)
                                                                                                                                  aaGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATg               <  1:37508/149‑1 (MQ=255)
                                                                                                                                    gCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCGGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGc             <  7:129389/149‑1 (MQ=255)
                                                                                                                                     ccggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCtt                                     <  6:117539/121‑1 (MQ=255)
                                                                                                                                     ccggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCtt                                     >  5:117539/4‑124 (MQ=255)
                                                                                                                                         ggCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                      >  1:8329/1‑135 (MQ=255)
                                                                                                                                         ggCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                      <  2:8329/135‑1 (MQ=255)
                                                                                                                                          gCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAAcc       >  2:13449/1‑149 (MQ=255)
                                                                                                                                                cTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCTGCCCCCTCGATAATGGCTGAACCCCggg  <  1:190153/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATCGCGCGCGCCGACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGG  >  NC_002947/4649840‑4650131

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