Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,695,470 T→C 8.2% I272V (ATC→GTC)  PP_4155 ← amine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,695,4700TC 8.2% 131.1 / 2.3 49I272V (ATC→GTC) PP_4155amine oxidase
Reads supporting (aligned to +/- strand):  ref base T (17/27);  new base C (2/2);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCATCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGCACAGGCCGACAGTTGCTCGACGG  >  NC_002947/4695322‑4695610
                                                                                                                                                    |                                                                                                                                            
cccATCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGAt                                                                                                                                              <  8:96333/149‑1 (MQ=255)
 ccATCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGAtg                                                                                                                                             >  3:251644/1‑149 (MQ=255)
  cATCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGAtgt                                                                                                                                            <  6:308356/149‑1 (MQ=255)
  cATCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGAtgt                                                                                                                                            <  8:25047/149‑1 (MQ=255)
    tCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGTAACGTTTTTATAACCGCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTgg                                                                                                                                          <  1:104833/149‑1 (MQ=255)
    tCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTgg                                                                                                                                          <  6:4546/149‑1 (MQ=255)
     ccccAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTggg                                                                                                                                         >  2:249261/1‑149 (MQ=255)
      cccAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGt                                                                                                                                        <  2:155175/149‑1 (MQ=255)
                ccAACGCTCACTGCTGTATTTCAGGTTACGCAGGGGTAATGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTg                                                                                                                              <  8:188288/149‑1 (MQ=255)
                         ccTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGc                                                                                                                       <  6:105767/147‑1 (MQ=255)
                         ccTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCg                                                                                                                      <  5:318454/148‑1 (MQ=255)
                         ccTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCg                                                                                                                      <  8:79581/148‑1 (MQ=255)
                          ctgctgTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCAACCGGGAATGTCCAAATGCACAAGTCGCACACGACTACATAAACCATGACGCGCGCCCTCGACAAGACGATTTGGGTGACAGGACTGTTCAGCCGGa                                                                                                                    >  1:161607/1‑149 (MQ=255)
                           tgctgTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGCCCAGTCGAGGTGGGTAACAGGACTGTTCAGGCGGa                                                                                                                    >  4:36915/1‑148 (MQ=255)
                           tgctgTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGACGTGGGTGAc                                                                                                                                     >  5:39286/1‑131 (MQ=255)
                           tgctgTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGACGTGGGTGAc                                                                                                                                     >  5:34549/1‑131 (MQ=255)
                           tgctgTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGACGTGGGTGAc                                                                                                                                     <  6:39286/131‑1 (MQ=255)
                           tgctgTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGACGTGGGTGAc                                                                                                                                     <  6:34549/131‑1 (MQ=255)
                             ctgTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGata                                                                                                                  <  7:121677/148‑1 (MQ=255)
                                tATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATc                                                                                                                >  4:101161/1‑147 (MQ=255)
                                  tttCAGTTAACGCATGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCATGCGGATATCCAg                                                                                                             <  7:282323/148‑1 (MQ=255)
                                  tttCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGtt                                                                                                                            <  2:154317/133‑1 (MQ=255)
                                  tttCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGtt                                                                                                                            >  1:154317/1‑133 (MQ=255)
                                  tttCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAg                                                                                                             <  6:241962/148‑1 (MQ=255)
                                    tCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGt                                                                                                           <  8:216190/148‑1 (MQ=255)
                                    tCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGt                                                                                                           <  6:50830/148‑1 (MQ=255)
                                      aGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGc                                                                                                         <  4:251644/148‑1 (MQ=255)
                                       ggTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACATGACTGCTCAGGCGGATATCCAGGTGCtt                                                                                                       <  7:203037/149‑1 (MQ=255)
                                        gTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTc                                                                                                      <  8:188942/149‑1 (MQ=255)
                                            cGCAGTGGTAAAGCGCAACGTTTTCCGAACCCCCACCGGGACTGTGCACATGCAGAAGTCGCACACGGCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCCGGCGGATATCCAGGTGCTTCGACa                                                                                                  <  8:209827/149‑1 (MQ=255)
                                               aGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGc                                                                                                >  4:91528/1‑148 (MQ=255)
                                                 gggTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCTCCATGACGCGCGCGCTCGACCAGTCGATGTGGGGGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAg                                                                                             <  5:257009/149‑1 (MQ=255)
                                                       aGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGAc                                                                                        >  4:16738/1‑148 (MQ=255)
                                                               gTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCt                                                                               <  3:16738/149‑1 (MQ=255)
                                                                          ccccACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTg                                                                    >  4:189760/1‑149 (MQ=255)
                                                                            ccACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTgg                                                                   >  1:137370/1‑148 (MQ=255)
                                                                               ccGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGgcag                                                                >  8:257388/1‑148 (MQ=255)
                                                                                cGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGgcagca                                                              <  3:91528/149‑1 (MQ=255)
                                                                                      tgtgCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGAt                                                        >  2:138948/1‑149 (MQ=255)
                                                                                        tgCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACTATc                                                       <  2:161607/148‑1 (MQ=255)
                                                                                          cAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCAcc                                                    <  3:101161/149‑1 (MQ=255)
                                                                                                         cacaCGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCctc                                     <  1:138948/149‑1 (MQ=255)
                                                                                                            acGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTAGCAGCATACGATCACCATGCCCGACCCCCTCtt                                   >  2:197850/1‑148 (MQ=255)
                                                                                                                cTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATc                               >  5:43178/1‑148 (MQ=255)
                                                                                                                    attcaCCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAg                           <  1:249261/148‑1 (MQ=255)
                                                                                                                         ccATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATGACCATGCCCGCCCCCCTCTGCATCCAAGGCACAg                     >  4:26678/1‑149 (MQ=255)
                                                                                                                          cATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGca                        >  6:97319/1‑145 (MQ=255)
                                                                                                                          cATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGca                        <  5:97319/145‑1 (MQ=255)
                                                                                                                                  gcgcgcTCGACCAGTCGATGAGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGAACAGGCCGACAGt            >  7:62791/1‑149 (MQ=255)
                                                                                                                                        tCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGCACAGGCCGACAGTTGCTc       >  3:148183/1‑148 (MQ=255)
                                                                                                                                           aCCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGCACAGGCCGACAGTTGCTCGACg   >  3:173499/1‑149 (MQ=255)
                                                                                                                                           aCCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGCACAGGCCGACAGTTGCTCGACg   <  2:137370/149‑1 (MQ=255)
                                                                                                                                            ccAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGCACAGGCCGACAGTTGCTCGACgg  <  6:43178/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCCATCCCCAGATGCGCCAACGCTCCCTGCTGTATTTCAGGTAACGCAGGGGTAAAGCGCAACGTTTTCAGAACCCCCACCGGGACTGTGCAAATGCAGAAGTCGCACACGTCTTCATTCACCATGACGCGCGCGCTCGACCAGTCGATGTGGGTGACAGGACTGTTCAGGCGGATATCCAGGTGCTTCGACAGCAAGTCGACAAGCTCGCTGTAGCCGCCTGGCAGCATACGATCACCATGCCCGACCCCCTCTTCATCCAAGGCACAGGCCGACAGTTGCTCGACGG  >  NC_002947/4695322‑4695610

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