Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,463,925 T→A 8.8% L130* (TTG→TAG)  PP_3956 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,463,9250TA 8.8% 152.1 / 2.8 53L130* (TTG→TAG) PP_3956hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/28);  new base A (3/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCCTTGCACGCCGGGCGAACGCCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCG  >  NC_002947/4463780‑4464071
                                                                                                                                                  |                                                                                                                                                        
cAGCCCTTGCACGCCGGGCGAACGCCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGa                                                                                                                                                       <  4:293078/149‑1 (MQ=255)
      ttGCACGCCGGGCGAACGCCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGAcagcag                                                                                                                                                 >  3:109520/1‑149 (MQ=255)
        gCACGCCGGGCGAACGCCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAgg                                                                                                                                                <  1:150087/148‑1 (MQ=255)
         cACGCCGGGCGAACGCCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGc                                                                                                                                               <  4:317029/148‑1 (MQ=255)
                ggCGAACGCCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCgc‑‑c                                                                                                                                     <  8:298180/149‑1 (MQ=255)
                   gAACGCCCGAGGCTGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCgc‑‑cctgc                                                                                                                                >  2:246304/1‑149 (MQ=255)
                               cTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGa                                                                                                                     <  8:82811/149‑1 (MQ=255)
                                   tGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAgg                                                                                                                <  6:14381/149‑1 (MQ=255)
                                   tGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAg                                                                                                                 <  1:246304/148‑1 (MQ=255)
                                    gACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAggg                                                                                                               >  2:275709/1‑149 (MQ=255)
                                      cTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAggg                                                                                                               >  2:191079/1‑147 (MQ=255)
                                       tACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAg                                                                                                           >  8:48777/1‑149 (MQ=255)
                                        aCCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAg                                                                                                           >  3:254592/1‑148 (MQ=255)
                                         ccAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGTCGCCGTCAGGGCAgg                                                                                                          <  1:62407/148‑1 (MQ=255)
                                         ccAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGa                                                                                                         >  8:275762/1‑149 (MQ=255)
                                          cAGGTCTGCACGTCTGACGGGTGCCAGCCACTGTCCTAAGGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGAc                                                                                                        <  1:284882/149‑1 (MQ=255)
                                             gTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCg                                                                                                      >  3:108245/1‑148 (MQ=255)
                                                     gACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACtt                                                                                            <  8:46479/149‑1 (MQ=255)
                                                        aGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCAc                                                                                              >  1:70256/1‑144 (MQ=255)
                                                        aGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCAc                                                                                              <  2:70256/144‑1 (MQ=255)
                                                                  cagccaCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGt                                                                               <  5:170788/149‑1 (MQ=255)
                                                                      caCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCaa                                                                           <  4:109520/149‑1 (MQ=255)
                                                                       aCTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAAATCCACTTCAGCCCACCACGTGCaa                                                                           <  7:48777/149‑1 (MQ=255)
                                                                                ttGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTAGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAgg                                                                 >  2:258741/1‑149 (MQ=255)
                                                                                ttGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTAGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAgg                                                                 >  2:262988/1‑149 (MQ=255)
                                                                                  ggTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGt                                                                <  5:176304/148‑1 (MQ=255)
                                                                                   gTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTAAGGAAAGCAGGCGAACCCC‑‑CCTGCCCTGACGGCCACAGGGCAGCACCGAATCCACTTCAGCCCACCACGTGCAAGTAATACAGGTg                                                               >  1:271350/1‑149 (MQ=255)
                                                                                   gTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTgg                                                              <  2:77153/149‑1 (MQ=255)
                                                                                   gTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTgg                                                              >  2:290018/1‑149 (MQ=255)
                                                                                        gACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTc                                                         <  8:184204/149‑1 (MQ=255)
                                                                                         aCCGCCGAGCTAACGAATTGCTGCAGTAGACTTCTGGGCTTGTTAAAATCCGCTTTTGTCCAGCAGGCGAACCCCTACACTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGt                                                           <  1:275709/149‑1 (MQ=255)
                                                                                           cgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGc                                                                    <  6:308699/135‑1 (MQ=255)
                                                                                           cgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTgg                                                              <  3:63799/141‑1 (MQ=255)
                                                                                           cgccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTgg                                                              >  4:63799/1‑141 (MQ=255)
                                                                                             ccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTAGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACATGGCCTGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                     <  1:262988/149‑1 (MQ=255)
                                                                                             ccgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTAGGACAGCAGGCGAAACGC‑‑CCTGCCCTGACGGCGACAGGGCTTGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                     <  1:258741/149‑1 (MQ=255)
                                                                                              cgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCcca                                                                                     <  7:306576/115‑1 (MQ=255)
                                                                                              cgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCcca                                                                                     >  8:306576/1‑115 (MQ=255)
                                                                                              cgAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                     <  1:290018/147‑1 (MQ=255)
                                                                                                 gCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATc                                                >  5:144475/1‑149 (MQ=255)
                                                                                                   gAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                     <  4:254540/142‑1 (MQ=255)
                                                                                                   gAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTAc                                                     >  3:254540/1‑142 (MQ=255)
                                                                                                      cGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGa                                           <  5:317388/149‑1 (MQ=255)
                                                                                                        cATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGCATg                                         >  4:120926/1‑149 (MQ=255)
                                                                                                         aTTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGt                                        >  7:10712/1‑149 (MQ=255)
                                                                                                                cAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAAt                                 >  6:345284/1‑149 (MQ=255)
                                                                                                                cAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAAt                                 <  7:305766/149‑1 (MQ=255)
                                                                                                                        ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGt                         >  6:184137/1‑149 (MQ=255)
                                                                                                                        ttCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGt                         >  7:216943/1‑149 (MQ=255)
                                                                                                                                       aaTCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCa          <  2:271350/149‑1 (MQ=255)
                                                                                                                                        aTCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCat         <  8:216943/149‑1 (MQ=255)
                                                                                                                                         tCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCatg        <  8:10712/149‑1 (MQ=255)
                                                                                                                                           cGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCatgat      >  2:74546/1‑149 (MQ=255)
                                                                                                                                                ttGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCAACCCCATGATGTCg  >  6:9878/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                        
CAGCCCTTGCACGCCGGGCGAACGCCCGAGGCTGGTGACTACCAGGTCTGCACGACAGACGGGTGCCAGCCACTGTCCTATTGGTTAGGACCGCCGAGCGAACGCATTGCTGCAGTAGACTTCTGGGCCTGTGAAAATCCGCCTTTGGACAGCAGGCGAACCGC‑‑CCTGCCCTGACGGCGACAGGGCAGGACCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCG  >  NC_002947/4463780‑4464071

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: