Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,101,394 A→G 8.1% I65T (ATT→ACT)  PP_3609 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,101,3940AG 8.1% 135.4 / 2.5 49I65T (ATT→ACT) PP_3609membrane protein
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base G (2/2);  total (26/23)
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.86e-01

ACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCGCACCGGCAAGCAACGCCATGGCCACAG  >  NC_002947/4101247‑4101541
                                                                                                                                                   |                                                                                                                                                   
accaccaGTACCAGGACCCCTGCTGCGCCAAGTTAAGGGGTGAGTGCGGCGGCGCCGGCGACATACATTGCCCATCACCCCCCACCAATCCACAACCACCAAGGACCCTGCAGGGCCTTGCCCATCTCATGAGCACGTACCCGGAGaa                                                                                                                                                     >  4:219015/1‑148 (MQ=255)
accaccaGTACCAGGAACCCTGCTGCGCCAAGTTTAGGCGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAAt                                                                                                                                                    <  7:265592/149‑1 (MQ=255)
    ccaGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCAACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATcag                                                                                                                                                 >  1:159375/1‑148 (MQ=255)
        tACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAggg                                                                                                                                            <  3:55729/149‑1 (MQ=255)
            aGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGc                                                                                                                                        <  5:8707/149‑1 (MQ=255)
            aGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTAAATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGCCCCATGCAGGGCCTTCCCCATGTCGAGCTCCGGAACCCCGATAATCAGCAGGGAAg                                                                                                                                         >  2:341110/1‑148 (MQ=255)
               aaCCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGAt                                                                                                                                     <  3:273314/149‑1 (MQ=255)
                aCCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGAt                                                                                                                                     >  5:85172/1‑148 (MQ=255)
                 cccTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATca                                                                                                                                   >  8:12350/1‑149 (MQ=255)
                 cccTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATca                                                                                                                                   <  6:85172/149‑1 (MQ=255)
                    tgctgcGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATcag                                                                                                                                                 >  5:90889/1‑132 (MQ=255)
                    tgctgcGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATcag                                                                                                                                                 <  6:90889/132‑1 (MQ=255)
                       tgcGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATcaccac                                                                                                                               >  6:130620/1‑147 (MQ=255)
                         cgcCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCAccacca                                                                                                                             <  2:82589/147‑1 (MQ=255)
                            cAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCGCCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATg                                                                                                                        >  4:167635/1‑149 (MQ=255)
                            cAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGACAATCAGCAGGGCAGCGATCACCACCAGCATg                                                                                                                        >  1:151772/1‑149 (MQ=255)
                                   aGGTGTGAGTGCGTCGTCGACTGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGACCATGTCTGGCGCCGGTACCCGTATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCAcc                                                                                                                  <  3:219015/148‑1 (MQ=255)
                                         gAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATc                                                                                                                                    <  6:331191/124‑1 (MQ=255)
                                         gAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATc                                                                                                                                    >  5:331191/1‑124 (MQ=255)
                                           gTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGc                                                                                                                       >  7:109478/1‑135 (MQ=255)
                                           gTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGc                                                                                                                       <  8:109478/135‑1 (MQ=255)
                                              cggcggcgCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGAtgtg                                                                                                      <  5:190291/149‑1 (MQ=255)
                                              cggcggcgCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGAtgt                                                                                                       >  7:341781/1‑148 (MQ=255)
                                                   gcgcCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGc                                                                                                 >  2:283556/1‑149 (MQ=255)
                                                    cgcCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAATGcc                                                                                                <  8:341781/149‑1 (MQ=255)
                                                              atacatTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGc                                                                                       >  2:225049/1‑148 (MQ=255)
                                                               tacatTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCa                                                                                     <  5:130620/149‑1 (MQ=255)
                                                                acatTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAg                                                                                                  <  1:148661/135‑1 (MQ=255)
                                                                     gCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGcc                                                                                <  1:341110/148‑1 (MQ=255)
                                                                                  caccaATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACCGCCAGTGCCCAAGTGCCCGACCt                                                                  >  4:312538/1‑149 (MQ=255)
                                                                                    ccaATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAGTCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGc                                                                >  5:115816/1‑149 (MQ=255)
                                                                                    ccaATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAGTCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGc                                                                >  7:136234/1‑149 (MQ=255)
                                                                                          cacaACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCgcag                                                           <  7:12350/148‑1 (MQ=255)
                                                                                           acaACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATcacc                                                                                                                                 <  6:349697/77‑1 (MQ=255)
                                                                                           acaACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATcacc                                                                                                                                 >  5:349697/1‑77 (MQ=255)
                                                                                               ccaccaAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTAc                                                     >  4:52302/1‑149 (MQ=255)
                                                                                                   caAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGccc                                                                                               >  2:41248/1‑103 (MQ=255)
                                                                                                   caAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGccc                                                                                               <  1:41248/103‑1 (MQ=255)
                                                                                                    aaGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACCGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCg                                                >  7:305255/1‑149 (MQ=255)
                                                                                                        ccccTGCAGGGCCTTGCCCATGTCTGGCGCCGGTTCCCGGATAGTCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCt                                            <  6:115816/149‑1 (MQ=255)
                                                                                                        ccccTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAGTCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCt                                            <  8:136234/149‑1 (MQ=255)
                                                                                                          ccTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGaa                                                                                                    >  6:23037/1‑91 (MQ=255)
                                                                                                          ccTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGaa                                                                                                    <  5:23037/91‑1 (MQ=255)
                                                                                                                     ttGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAacac                               >  2:189310/1‑149 (MQ=255)
                                                                                                                       gCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGTCAGCTATCACCACCATCATGCTCACCGTGAACGATGTGAAAGCCCCCCACAGCGAGTGCACAAGTGACCGACCTGCCGCAGCATTACCTGCGTCCTGAAACGGCAacac                               <  2:159375/147‑1 (MQ=255)
                                                                                                                         ccATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTg                                                  <  2:48294/126‑1 (MQ=255)
                                                                                                                         ccATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAAcaccgca                           >  4:222566/1‑149 (MQ=255)
                                                                                                                                tGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCGCACCGgc                      >  5:165381/1‑147 (MQ=255)
                                                                                                                                    gcCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCGCACCGgcaagcaa                >  4:309659/1‑149 (MQ=255)
                                                                                                                                            ccGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCGCACCGGCAAGCAACGCCATgg        >  3:171681/1‑149 (MQ=255)
                                                                                                                                                   aTCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCGCACCGGCAAGCAACGCCATGGCCACAg  <  3:167635/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCACCAGTACCAGGAACCCTGCTGCGCCAAGTTTAGGGGTGAGTGCGGCGGCGCCGGCTACATACATTGCCCCGAGCACGCCACCAATCCACAACCACCAAGGCCCCTGCAGGGCCTTGCCCATGTCTGGCGCCGGTACCCGGATAATCAGCAGGGCAGCGATCACCACCAGCATGCTCACCGTCAACGATGTGAAAGCCCCCCACAGCCAGTGCCCAAGTGCCCGACCTGCCGCAGCATTACCTGCGGCCTGAAACGGCAACACCGCACCGGCAAGCAACGCCATGGCCACAG  >  NC_002947/4101247‑4101541

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