Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,507,789 T→G 15.9% F32C (TTC→TGC) ‡ PP_3105 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,507,7890TG15.9% 90.3 / 13.8 44F32C (TTC→TGC) ‡PP_3105membrane protein
Reads supporting (aligned to +/- strand):  ref base T (17/20);  new base G (4/3);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

TTTCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCGCTGCTGACA  >  NC_002947/3507645‑3507937
                                                                                                                                                |                                                                                                                                                    
tttCGTGGCTGTGTACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTTGCGCTGCCATTCATgg                                                                                                                                                 <  2:929479/150‑1 (MQ=255)
tttCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATgg                                                                                                                                                 <  1:1854932/150‑1 (MQ=255)
                tACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCACCGCCATGAATGGCggc                                                                                                                                             >  2:57588/1‑137 (MQ=255)
                tACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCACCGCCATGAATGGCggc                                                                                                                                             <  1:57588/138‑2 (MQ=255)
                  cccTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGa                                                                                                                               <  2:30235/150‑1 (MQ=255)
                  cccTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGa                                                                                                                               <  2:1862609/150‑1 (MQ=255)
                        gTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTa                                                                                                                         <  1:272121/150‑1 (MQ=255)
                                 tGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGa                                                                                                                <  1:1363496/150‑1 (MQ=255)
                                  gACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGAc                                                                                                               >  1:1359406/1‑150 (MQ=255)
                                     ggggAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCt                                                                                                            <  1:973403/150‑1 (MQ=255)
                                      gggAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCtg                                                                                                           <  1:871136/150‑1 (MQ=255)
                                       ggAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTg                                                                                                                                            <  2:1446966/116‑1 (MQ=255)
                                       ggAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTg                                                                                                                                            >  1:1446966/1‑116 (MQ=255)
                                        gAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCtgtg                                                                                                         >  2:810361/1‑150 (MQ=255)
                                           ggATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTgtggtg                                                                                                      >  2:1806909/1‑150 (MQ=255)
                                              ttAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATGAATGGCggc                                                                                                                                             >  1:1866570/1‑107 (MQ=255)
                                              ttAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATGAATGGCggc                                                                                                                                             <  2:1866570/108‑2 (MQ=255)
                                                 aaTGTTTACCACAGAACGTTTTTTTTTTTTTTTTTGTTCGGTGTGGCTGCTGGGCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGTCGGTGCCAGCGATCCTGATCATTTTTTTTATTACTCTGTGGTGCATCaa                                                                                                >  2:456651/1‑150 (MQ=255)
                                                 aaTGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGa                                                                                                                <  1:341054/134‑1 (MQ=255)
                                                 aaTGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGa                                                                                                                >  2:341054/1‑134 (MQ=255)
                                                 aaTGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCaa                                                                                                >  1:1200270/1‑150 (MQ=255)
                                                   tGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                              <  1:1429608/150‑1 (MQ=255)
                                                   tGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                              <  1:1806909/150‑1 (MQ=255)
                                                   tGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                              <  1:1834964/150‑1 (MQ=255)
                                                           acaGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTAc                                                                                      >  2:1876880/1‑150 (MQ=255)
                                                               aaCGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCa                                                                                  >  1:431289/1‑150 (MQ=255)
                                                                  gTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTg                                                                               >  1:61921/1‑150 (MQ=255)
                                                                       tttAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGtt                                                                                                                    <  1:132214/108‑1 (MQ=255)
                                                                       tttAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGtt                                                                                                                    >  2:132214/1‑108 (MQ=255)
                                                                         tAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCNCtt                                                                        <  2:1200270/150‑1 (MQ=255)
                                                                                aTCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGt                                                                 <  2:431289/150‑1 (MQ=255)
                                                                                        ggTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGGGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGaa                                                         <  2:11399/150‑1 (MQ=255)
                                                                                           gtgGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAAttt                                                      <  1:1855764/150‑1 (MQ=255)
                                                                                                tgctgGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCa                                                 >  2:509429/1‑150 (MQ=255)
                                                                                                           atatTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAat                                      >  1:570294/1‑150 (MQ=255)
                                                                                                                 gCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCg                                >  2:1581292/1‑150 (MQ=255)
                                                                                                                   ccTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGt                              >  1:1824050/1‑150 (MQ=255)
                                                                                                                     tgatgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTg                            <  2:1824050/150‑1 (MQ=255)
                                                                                                                              ccGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGt                   >  2:1063967/1‑150 (MQ=255)
                                                                                                                                   gTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCa              >  2:1649093/1‑150 (MQ=255)
                                                                                                                                    cggcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCAt             >  2:1845966/6‑150 (MQ=255)
                                                                                                                                     ggcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                                     <  2:1059095/89‑1 (MQ=38)
                                                                                                                                     ggCACCGCCATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                                     >  1:1059095/1‑93 (MQ=38)
                                                                                                                                           gccATTCATGGCGGTGCCAGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCgctgct      <  1:456651/150‑1 (MQ=255)
                                                                                                                                               ttCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCCATATGTCGGTTGCACCAATGTTCGCATCGCTGCTGaca  <  1:520886/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTTCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCGCTGCTGACA  >  NC_002947/3507645‑3507937

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