Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,507,790 C→A 15.6% F32L (TTC→TTA) ‡ PP_3105 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,507,7900CA15.6% 106.1 / 11.4 45F32L (TTC→TTA) ‡PP_3105membrane protein
Reads supporting (aligned to +/- strand):  ref base C (17/21);  new base A (4/3);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

TTTCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCGCTGCTGACAC  >  NC_002947/3507645‑3507938
                                                                                                                                                 |                                                                                                                                                    
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:1834964/150‑1 (MQ=255)
                                                   tGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                               <  1:1806909/150‑1 (MQ=255)
                                                   tGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                               <  1:1429608/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)
                                                                                                                                                gcATTGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCGCTGCTGacac  >  1:366954/2‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTTCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCGCTGCTGACAC  >  NC_002947/3507645‑3507938

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: