Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,280,831 T→C 8.0% I77T (ATC→ACC)  PP_5742 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,280,8310TC 8.0% 139.6 / 2.9 50I77T (ATC→ACC) PP_5742membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base C (2/2);  total (24/26)
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.46e-01

GAAGAAGTGGCCAGACAGTGTTGCTGATGTTATTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATTGTGGCTGTTGAGGAGCGTGGCGATTATAT  >  NC_002947/1280691‑1280978
                                                                                                                                            |                                                                                                                                                    
gaagaaGTGGCCAGACAGTGTTGCTGATGTTATTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAAttt                                                                                                                                              >  8:11479/1‑149 (MQ=255)
          ccAGACAGTGTTGCTGATGTTATTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTc                                                                                                                                    <  5:254870/149‑1 (MQ=255)
               cAGTGTTGCTGATGTTATTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGaaaa                                                                                                                               <  3:137173/149‑1 (MQ=255)
                      gCTGATGTTATTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGaa                                                                                                                           <  4:102192/146‑1 (MQ=255)
                      gCTGATGTTATTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGaa                                                                                                                           >  3:102192/1‑146 (MQ=255)
                               aTTGGCACCCTGTGGTCATGGTTAACTGCGCTTTCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGaa                                                                                                                           <  3:163535/137‑1 (MQ=255)
                               aTTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATATTTTTGTATCGTTCGAAAATGAAGTTCTCATCCCa                                                                                                               >  1:88901/1‑149 (MQ=255)
                               aTTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGaa                                                                                                                           >  4:163535/1‑137 (MQ=255)
                                   gCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGaa                                                                                                                           <  6:216761/133‑1 (MQ=255)
                                   gCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGaa                                                                                                                           >  5:216761/1‑133 (MQ=255)
                                              tCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGctct                                                                                                >  7:146798/1‑149 (MQ=255)
                                                   gTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTcgcg                                                                                            >  5:227598/1‑148 (MQ=255)
                                                       aCTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAAt                                                                                        <  6:227598/148‑1 (MQ=255)
                                                           cgcTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGa                                                                                     <  8:142558/147‑1 (MQ=255)
                                                            gcTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATc                                                                                  <  5:141385/149‑1 (MQ=255)
                                                                  tGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATg                                                                                                                             <  1:97564/100‑1 (MQ=255)
                                                                  tGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATg                                                                                                                             >  2:97564/1‑100 (MQ=255)
                                                                  tGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCg                                                                                >  1:149791/1‑145 (MQ=255)
                                                                   gAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATc                                                                                                                  <  4:45771/110‑1 (MQ=255)
                                                                   gAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATc                                                                                                                  >  3:45771/1‑110 (MQ=255)
                                                                        gTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCt                                                                       <  2:149791/148‑1 (MQ=255)
                                                                              tCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGc                                                                   >  3:155206/1‑146 (MQ=255)
                                                                                aaaCTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGTCTTCTGTTAGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGa                                                               <  5:81015/148‑1 (MQ=255)
                                                                                aaaCTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGAt                                                              >  2:153349/1‑149 (MQ=255)
                                                                                aaaCTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTACCATAATTTTGTATCGtt                                                                                                                                     <  2:129800/78‑1 (MQ=255)
                                                                                aaaCTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTACCATAATTTTGTATCGtt                                                                                                                                     >  1:125605/1‑78 (MQ=255)
                                                                                aaaCTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTACCATAATTTTGTATCGtt                                                                                                                                     <  2:125605/78‑1 (MQ=255)
                                                                                aaaCTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTACCATAATTTTGTATCGtt                                                                                                                                     >  1:129800/1‑78 (MQ=255)
                                                                                      tCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGt                                                                                                                         <  4:95026/84‑1 (MQ=255)
                                                                                      tCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGt                                                                                                                         >  3:95026/1‑84 (MQ=255)
                                                                                      tCATTTTACCGATATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATgaga                                                          <  1:237521/147‑1 (MQ=255)
                                                                                        aTTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTca                                                      <  2:88901/149‑1 (MQ=255)
                                                                                             aCCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTcagcag                                                  >  7:159414/1‑148 (MQ=255)
                                                                                                 cTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAg                                              >  2:28188/1‑148 (MQ=255)
                                                                                                   aTGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTAc                                           >  2:97082/1‑149 (MQ=255)
                                                                                                      gATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGaa                                                                                          <  2:42244/99‑1 (MQ=255)
                                                                                                      gATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGaa                                                                                          >  1:42244/1‑99 (MQ=255)
                                                                                                      gATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTg                                         <  1:28188/148‑1 (MQ=255)
                                                                                                             cTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTa                                  <  7:11479/148‑1 (MQ=255)
                                                                                                              tatcgtgATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTAtt                               <  7:128493/143‑1 (MQ=255)
                                                                                                              tAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTAtt                               <  1:97082/149‑1 (MQ=255)
                                                                                                                gCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTAta                               >  8:128493/1‑146 (MQ=255)
                                                                                                                gCGGGATTTGGGCTTCTGATCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTAt                                >  5:39638/1‑146 (MQ=255)
                                                                                                                gCGGGATTTGGGCTTCTGATCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTAt                                <  6:39638/146‑1 (MQ=255)
                                                                                                                 cGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATtgtg                            >  6:41873/1‑149 (MQ=255)
                                                                                                                  gggATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATTGTgg                           >  4:47214/1‑149 (MQ=255)
                                                                                                                          ggCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATTGTGGCTGTTga                    <  8:146798/148‑1 (MQ=255)
                                                                                                                                  ttCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATTGTGGCTGTTGAGGAGCGt             <  5:41873/147‑1 (MQ=255)
                                                                                                                                      ctgctATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATTGTGGCTGTTGAGGAGCGTGGCGAt       >  3:15126/1‑149 (MQ=255)
                                                                                                                                         ctATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTTCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTACTTGTGGCTGTTGAGGAGCGTGGCGATta     <  1:153349/149‑1 (MQ=255)
                                                                                                                                            tCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATTGTGGCTGTTGAGGAGCGTGGCGATTatat  <  4:15126/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GAAGAAGTGGCCAGACAGTGTTGCTGATGTTATTGGCACCCTGTGGTCATGGTTAACTGCGCTTGCTGAATGGTTTACTCAAACTATCATTTTACCGCTATGGATTTTACTAGCGGGATTTGGGCTTCTGTTCGCTGCTATCATAATTTTGTATCGTTCGAAAATGAAGTTCTCATCCCAGTTGGAAAGCGCTCTCGCGAAATTGAATCCGAAGCTTCCTGTGCTAGATGAGAGTCAGCAGAAAGTACTGGCAGTTATTGTGGCTGTTGAGGAGCGTGGCGATTATAT  >  NC_002947/1280691‑1280978

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: