Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 960,082 C→T 8.7% I117I (ATC→ATT)  PP_0821 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947960,0820CT 8.7% 128.4 / 2.7 46I117I (ATC→ATT) PP_0821hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base T (2/2);  total (24/22)
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.02e-01

GCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGGTGTCATGCTCCGCGTGGCC  >  NC_002947/959935‑960229
                                                                                                                                                   |                                                                                                                                                   
gcagcTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGcatca                                                                                                                                                    >  5:327166/1‑149 (MQ=255)
  agcTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATg                                                                                                                                                  <  1:289839/149‑1 (MQ=255)
                   gATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGAcc                                                                                                                                  <  6:327166/148‑1 (MQ=255)
                   gATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACcg                                                                                                                                 >  5:127624/1‑149 (MQ=255)
                          ggCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCg                                                                                                                          >  5:264274/1‑149 (MQ=255)
                                   agagCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAg                                                                                                                 >  7:23696/1‑149 (MQ=255)
                                    gagCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGc                                                                                                                >  5:104440/1‑149 (MQ=255)
                                      gCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCAGCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCgg                                                                                                              <  6:259255/149‑1 (MQ=255)
                                        cTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGt                                                                                                            <  8:23696/149‑1 (MQ=255)
                                           cGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCt                                                                                                         <  6:264274/149‑1 (MQ=255)
                                            gggATGTTCAGCGCAAGTGGCTGGCATACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGTCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCtt                                                                                                        <  2:270307/149‑1 (MQ=255)
                                                  ttCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCt                                                                                                  <  7:171057/149‑1 (MQ=255)
                                                   tCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCcgc                                                                                                                                      <  2:27018/112‑1 (MQ=255)
                                                   tCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCcgc                                                                                                                                      >  1:27018/1‑112 (MQ=255)
                                                     aGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTggg                                                                                               <  6:127624/149‑1 (MQ=255)
                                                      gcgcAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATTATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTggg                                                                                               >  4:291319/1‑148 (MQ=255)
                                                      gcgcAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATTATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTggg                                                                                               >  4:296116/1‑148 (MQ=255)
                                                        gcAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAg                                                                                            >  6:321897/1‑149 (MQ=255)
                                                               gctggcCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTa                                                                                                       >  1:328482/1‑131 (MQ=255)
                                                               gctggcCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTa                                                                                                       <  2:328482/131‑1 (MQ=255)
                                                                ctggcCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATTATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCt                                                                                    <  3:291319/149‑1 (MQ=255)
                                                                ctggcCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATTATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCt                                                                                    <  3:296116/149‑1 (MQ=255)
                                                                   gcCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTgaa                                                                                 <  6:282847/149‑3 (MQ=255)
                                                                   gcCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTg                                                                                   <  1:65430/147‑1 (MQ=255)
                                                                   gcCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGtt                                                                                 <  6:104440/149‑1 (MQ=255)
                                                                      tACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAg                                                                              >  5:336768/1‑149 (MQ=255)
                                                                          gTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCAt                                                                          <  5:198270/149‑1 (MQ=255)
                                                                               ggCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGCTTg                                                                     >  8:26848/1‑149 (MQ=255)
                                                                                  aaTTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGcc                                                                  >  2:67467/1‑149 (MQ=255)
                                                                                      gCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCcaca                                                              >  4:305470/1‑149 (MQ=255)
                                                                                         aGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACAttt                                                           >  1:119631/1‑149 (MQ=255)
                                                                                         aGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACAttt                                                           <  1:67467/149‑1 (MQ=255)
                                                                                                   aCAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCgg                                                 >  6:280569/1‑149 (MQ=255)
                                                                                                    cAGGCCAGTGGCTGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCggg                                                <  6:269363/149‑1 (MQ=255)
                                                                                                     aGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACAt                                                             >  2:216208/1‑135 (MQ=255)
                                                                                                     aGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACAt                                                             <  1:216208/135‑1 (MQ=255)
                                                                                                         cAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAg                                           >  1:299487/1‑149 (MQ=255)
                                                                                                              gCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCcgc                                       >  6:51225/1‑148 (MQ=255)
                                                                                                               cGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCcgcg                                      >  8:176324/1‑148 (MQ=255)
                                                                                                                  gACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAaca                                  >  2:155295/1‑149 (MQ=255)
                                                                                                                      cTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCgg                              <  2:119631/149‑1 (MQ=255)
                                                                                                                       tGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGc                             >  2:253495/1‑149 (MQ=255)
                                                                                                                         gcgcAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGc                             <  7:176324/147‑1 (MQ=255)
                                                                                                                           gcAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTg                         >  1:315373/1‑149 (MQ=255)
                                                                                                                            cAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGc                        >  1:236013/1‑149 (MQ=255)
                                                                                                                                     ggTGGCAGCTGAATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGGTGTCAt               <  5:51225/149‑1 (MQ=255)
                                                                                                                                         gCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGGTGTCATGCTc           >  8:122085/1‑149 (MQ=255)
                                                                                                                                                 atcatGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGgtg                   <  2:100011/133‑1 (MQ=255)
                                                                                                                                                 atcatGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGgtg                   >  1:100011/1‑133 (MQ=255)
                                                                                                                                                  tcatGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGGTGTCATGCTCCGCGTGGcc  <  5:321897/149‑1 (MQ=255)
                                                                                                                                                  tcatGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGGTGTCATGCTCCGCGTGGcc  >  3:129839/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAGCTGATGGGCAAGCTGGATGCCGGGCAGCAGAAGAGCCTGCGGGATGTTCAGCGCAAGTGGCTGGCCTACCGTGATGGCAATTGCCAGTTCCATGTACAGGCCAGTGGCGGGACCCTGGCGCAGCTGGAAGGTGGCAGCTGCATCATGGACATGACCCGCGACCGCGCTGCGGAACTGGAGCGGGTGCTTAGCCCTGGGCAGTAGCGCCTGTTCAGGCATGATTGGCCCACATTTGCAGGCGCGGGTTCAGCCGCGAACACCGGCTTTGCCGGTGTCATGCTCCGCGTGGCC  >  NC_002947/959935‑960229

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