Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,187,362 A→G 8.9% C83R (TGT→CGT)  PP_1936 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,187,3620AG 8.9% 120.4 / 2.9 45C83R (TGT→CGT) PP_1936membrane protein
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

CATAGAATAGTCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGAATTATAGGTGCGTAATTTGATCCTAAAT  >  NC_002947/2187226‑2187505
                                                                                                                                        |                                                                                                                                               
cATAGAATAGTCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGAt                                                                                                                                        >  6:218227/1‑146 (MQ=255)
cATAGAATAGTCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGAt                                                                                                                                        <  5:218227/146‑1 (MQ=255)
   aGAATAGTCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCg                                                                                                                                  >  4:68141/1‑149 (MQ=255)
      aTAGTCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGa                                                                                                                                <  4:52044/148‑1 (MQ=255)
       tAGTCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGa                                                                                                                                >  4:170451/1‑147 (MQ=255)
          tCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAgtg                                                                                                                           >  4:108942/1‑149 (MQ=255)
            aaGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTc                                                                                                                         >  5:285192/1‑149 (MQ=255)
            aaGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTc                                                                                                                         >  4:98308/1‑149 (MQ=255)
                  gCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCgg                                                                                                                    <  3:68141/148‑1 (MQ=255)
                   cAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGt                                                                                                                   <  1:138856/148‑1 (MQ=255)
                         cccAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTcc                                                                                                             >  6:72181/1‑148 (MQ=255)
                         cccAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTcc                                                                                                             <  5:72181/148‑1 (MQ=255)
                          ccAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAgt                                                                                                                            <  1:128713/132‑1 (MQ=255)
                             gTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGCTGCTAGCAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAACTTTTCATTATGCGCACAGTGTCTACGGTTGCTCCCctct                                                                                                        >  4:238645/1‑149 (MQ=255)
                                cGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTga                                                                                                     <  3:170451/149‑1 (MQ=255)
                                  gATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGAAAGTGTCTTCGGTTGCTCCCCGCTTga                                                                                                     >  7:243565/1‑147 (MQ=255)
                                      cGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGaa                                                                                                <  3:98308/148‑1 (MQ=255)
                                               aaaCACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATg                                                                                       >  8:67437/1‑148 (MQ=255)
                                               aaaCACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCACAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAAtaattaaag                                                                                       >  6:329007/1‑146 (MQ=255)
                                                        aTGGAGCGATTTTAGTCTCGGTGGTAGGAGCGAATTTGCCCAGAGACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCGCGCTTGAGAGTTTAATTAATGCATaaaagc                                                                              >  4:222668/1‑146 (MQ=255)
                                                           gAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAAt                                                                           <  3:238645/148‑1 (MQ=255)
                                                            aGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATgg                                                                         >  4:71647/1‑149 (MQ=255)
                                                             gCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGAt                                                                                                                                        >  1:262250/1‑85 (MQ=255)
                                                             gCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGAt                                                                                                                                        <  2:262250/85‑1 (MQ=255)
                                                                 ttttAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAAtaattaat                                                                                        <  5:149529/129‑1 (MQ=255)
                                                                 ttttAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAAtaattaat                                                                                        >  6:149529/1‑129 (MQ=255)
                                                                      gTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAg                                                                <  3:71647/148‑1 (MQ=255)
                                                                         tcGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACGGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTga                                                                                                     <  8:139213/108‑1 (MQ=255)
                                                                         tcGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACGGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTga                                                                                                     >  7:139213/1‑108 (MQ=255)
                                                                         tcGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACGGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTga                                                                                                     <  2:287488/108‑1 (MQ=255)
                                                                         tcGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACGGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTga                                                                                                     >  1:287488/1‑108 (MQ=255)
                                                                          cGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGt                                                             >  2:140057/1‑147 (MQ=255)
                                                                           ggTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGTTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAg                                                          <  1:12893/149‑1 (MQ=255)
                                                                             tGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAggg                                                        <  6:285192/149‑1 (MQ=255)
                                                                              gCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGc                                                       >  8:294394/1‑149 (MQ=255)
                                                                                  ggAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTg                                                    >  4:7005/1‑148 (MQ=255)
                                                                                       gAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCa                                              <  7:67437/149‑1 (MQ=255)
                                                                                            tgaccAGATACCTAAAATTTATGCATCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACACTGTCTTCGGTTGCTCCCCGCTTTAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATg                                          <  5:329007/145‑1 (MQ=255)
                                                                                                        gAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGaa                             <  1:140057/149‑1 (MQ=255)
                                                                                                            aTTAATGAAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGAATta                          >  3:144198/1‑148 (MQ=255)
                                                                                                               aaTGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGAATTATAgg                      >  1:302834/1‑149 (MQ=255)
                                                                                                                  gCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGAATTATAGGTGc                   >  6:305499/1‑149 (MQ=255)
                                                                                                                              tCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAggg                                                        <  8:22787/100‑1 (MQ=255)
                                                                                                                              tCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGAATTATAGGTGCGTAATTTGATc        <  3:222668/148‑1 (MQ=255)
                                                                                                                                    cTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGAATTATAGGTGCGTAATTTGATCCTAAAt  <  3:108942/148‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                               
CATAGAATAGTCAAGATTGCAGGGGCCCAGTACGGATCCGCTGCCTCAAACACCTCATGGAGCGATTTTAGTCTCGGTGCTAGGAGCGAATTTGCCCAGATACCGAAAATTAATGCAGCGAGAAGATCTATTCTACAGCTTTTGATTATGCGGACAGTGTCTTCGGTTGCTCCCCGCTTGAGAGAATAATTAATGCATAAAATCAATGGCATTGAAAGGGTTAGGGCTTGAACCCAGATGTTTTCTTTTAGAATTATAGGTGCGTAATTTGATCCTAAAT  >  NC_002947/2187226‑2187505

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