Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,148,799 G→T 18.0% T151K (ACG→AAG)  PP_1007 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,148,7990GT18.0% 123.9 / 16.6 50T151K (ACG→AAG) PP_1007transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base T (5/4);  total (25/25)
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.24e-01

TCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTG  >  NC_002947/1148655‑1148942
                                                                                                                                                |                                                                                                                                               
tCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAgg                                                                                                                                            >  1:301224/1‑150 (MQ=255)
  ggggACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAggtg                                                                                                                                          <  2:719905/150‑1 (MQ=255)
   gtgACCTGGAAGTAAGGCTCGCCCTGTAGCAGGCGGGCGACGTGGCGTCCGTGCTGGCGTGCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCTAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAggtgg                                                                                                                                         <  2:166612/148‑1 (MQ=255)
        cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCTTACGGTGGTCGGc                                                                                                                                    <  2:397682/150‑1 (MQ=255)
        cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAggtg                                                                                                                                          <  1:218265/144‑1 (MQ=255)
        cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAggtg                                                                                                                                          >  2:218265/1‑144 (MQ=255)
        cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGc                                                                                                                                    >  1:266054/1‑150 (MQ=255)
          ggAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCt                                                                                                                                  >  2:140521/1‑150 (MQ=255)
                   ccTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCa                                                                                                                         >  2:556741/1‑150 (MQ=255)
                          cTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                  <  1:724271/150‑1 (MQ=255)
                           tGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCa                                                                                                                               >  2:579425/1‑136 (MQ=255)
                           tGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCa                                                                                                                               <  1:579425/136‑1 (MQ=255)
                                     gggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACgg                                                                                                       <  2:648495/150‑1 (MQ=255)
                                                    cGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACca                                                                                        >  1:427598/1‑150 (MQ=255)
                                                     gTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcag                                                                                       <  2:184748/150‑1 (MQ=255)
                                                                            ccGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTTCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                >  2:559544/1‑150 (MQ=255)
                                                                            ccGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                <  2:427598/150‑1 (MQ=255)
                                                                                 ttGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACTACGGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTc                                                           >  2:707040/1‑150 (MQ=255)
                                                                                     tgtTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGt                                                       <  2:301224/150‑1 (MQ=255)
                                                                                           gcagcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTgc                                                 >  1:224632/1‑150 (MQ=255)
                                                                                             agcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTg                                                                       >  2:438418/1‑126 (MQ=255)
                                                                                             agcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTg                                                                       <  1:438418/126‑1 (MQ=255)
                                                                                               caCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAg                                             <  2:592043/150‑1 (MQ=255)
                                                                                                 cTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgc                                           <  1:140521/150‑1 (MQ=255)
                                                                                                       cgcCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCcgc                                     >  1:155832/1‑150 (MQ=255)
                                                                                                        gcCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                  >  2:80393/1‑72 (MQ=255)
                                                                                                        gcCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                  <  1:80393/72‑1 (MQ=255)
                                                                                                         ccGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCgcgc                                   <  1:556741/150‑1 (MQ=255)
                                                                                                          cGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGcc                                  >  2:233303/1‑150 (MQ=255)
                                                                                                           gTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCa                                 >  1:576579/1‑150 (MQ=255)
                                                                                                            tCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTa                                                                                                              >  1:566016/1‑72 (MQ=255)
                                                                                                            tCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTa                                                                                                              <  2:566016/72‑1 (MQ=255)
                                                                                                                ccAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGAcc                            >  1:523955/1‑150 (MQ=255)
                                                                                                                   gTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCAGGACCTGc                         >  2:278617/1‑150 (MQ=255)
                                                                                                                       cAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtg                     <  2:135859/150‑1 (MQ=255)
                                                                                                                                 cgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGc           <  2:576579/150‑1 (MQ=255)
                                                                                                                                  gctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGAAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCa          >  1:537647/1‑150 (MQ=255)
                                                                                                                                      gCCGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCAtgtg      >  2:456756/1‑150 (MQ=255)
                                                                                                                                       ccGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGc     >  2:558301/1‑150 (MQ=255)
                                                                                                                                       ccGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGc     <  2:155832/150‑1 (MQ=255)
                                                                                                                                        cGACCACCTTACGGTGGTCGGCGAGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCg    >  2:255735/1‑150 (MQ=255)
                                                                                                                                        cGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGa                                                     <  1:80397/101‑1 (MQ=255)
                                                                                                                                        cGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGa                                                     >  2:80397/1‑101 (MQ=255)
                                                                                                                                         gACCACCTTACGGTGGTCGGCGAGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGt   <  1:255735/150‑1 (MQ=255)
                                                                                                                                         gACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGt             >  1:622357/1‑140 (MQ=255)
                                                                                                                                         gACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGt             <  2:622357/140‑1 (MQ=255)
                                                                                                                                         gACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGt   <  1:278617/150‑1 (MQ=255)
                                                                                                                                          accaccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTg  <  1:233303/150‑1 (MQ=255)
                                                                                                                                            caccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACgg                                                                                                       <  1:751842/47‑1 (MQ=255)
                                                                                                                                            caccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACgg                                                                                                       >  2:751842/1‑47 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTG  >  NC_002947/1148655‑1148942

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