Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,892 T→G 12.1% Q631P (CAA→CCA)  PP_1033 ← sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,8920TG12.1% 79.5 / 3.5 37Q631P (CAA→CCA) PP_1033sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (12/20);  new base G (3/2);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAA  >  NC_002947/1178751‑1179022
                                                                                                                                             |                                                                                                                                  
cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGc                                                                                                                            <  2:83010/150‑1 (MQ=255)
   gCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCttt                                                                                                                         >  1:528971/1‑150 (MQ=255)
   gCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgacga                                                                                                                         <  1:569471/150‑5 (MQ=255)
   gCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTga                                                                                                                            <  2:321454/147‑2 (MQ=255)
   gCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTga                                                                                                                            >  1:321454/1‑146 (MQ=255)
      tGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTcc                                                                                                                      <  1:393684/150‑1 (MQ=255)
           gCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                 <  2:120660/150‑1 (MQ=255)
                   tCTGGTTCTTGTTCGTTCAGATTGCTTATTCCTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTc                                                                                                                       >  1:186969/1‑136 (MQ=255)
                   tCTGGTTCTTGTTCGTTCAGATTGCTTATTCCTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTc                                                                                                                       <  2:186969/136‑1 (MQ=255)
                              ttcgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTAt                                                                                              <  1:661621/150‑1 (MQ=255)
                                       aTTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                     <  2:528971/150‑1 (MQ=255)
                                           cTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg                                                                                                                                 <  1:210778/102‑1 (MQ=255)
                                           cTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg                                                                                                                                 >  2:210778/1‑102 (MQ=255)
                                            ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTc                                                                                <  1:51869/150‑1 (MQ=255)
                                            ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTc                                                                                >  2:599987/1‑150 (MQ=255)
                                            ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCGGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTc                                                                                <  1:105998/150‑1 (MQ=255)
                                                  ttctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACt                                                                          >  2:580591/1‑150 (MQ=255)
                                                     tGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg                                                                                                                                >  2:294548/1‑93 (MQ=255)
                                                     tGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg                                                                                                                                <  1:294548/93‑1 (MQ=255)
                                                             tGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATc                                                                                                           <  2:250325/106‑1 (MQ=255)
                                                             tGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATc                                                                                                           >  1:250325/1‑106 (MQ=255)
                                                               gcgaaTAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATg                                                             >  1:595363/4‑150 (MQ=255)
                                                                          tttCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCTATGAACTCGTCACCCTCCACGCGTGTGCGGTGCTTTTCTTGCTGGTTAACCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATATATTTGTACgg                                                  >  1:617166/1‑150 (MQ=255)
                                                                          tttCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACgg                                                  <  1:346135/150‑1 (MQ=255)
                                                                                  aGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTcca                                          >  2:247192/1‑150 (MQ=255)
                                                                                       ttATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAg                                     <  1:580591/150‑1 (MQ=255)
                                                                                        tataGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCTTTGGTGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTTGTCAGGGTGTTCTAACTGTATTCGTATATGTATTTGTACGGGTGCTCCACCAAgg                                    <  2:521659/150‑1 (MQ=255)
                                                                                             gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTc                                                                                        >  1:37524/1‑93 (MQ=255)
                                                                                             gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTc                                                                                        <  2:37524/93‑1 (MQ=255)
                                                                                                            ctcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAg                <  2:617166/150‑1 (MQ=255)
                                                                                                               aCTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAg                                     <  1:421390/126‑1 (MQ=255)
                                                                                                               aCTGCGATGAACTCGTCACCCTCCACGCGTGGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAg                                     >  2:421390/1‑126 (MQ=255)
                                                                                                                  gCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGCCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTg          <  2:651463/150‑1 (MQ=255)
                                                                                                                          cTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATgaa  >  2:642319/1‑150 (MQ=255)
                                                                                                                          cTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGACAGTCGGTCAAAGAACGTGGACATgaa  <  1:606266/150‑1 (MQ=255)
                                                                                                                             gTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTgg                      <  1:334074/127‑1 (MQ=255)
                                                                                                                             gTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTgg                      >  2:334074/1‑127 (MQ=255)
                                                                                                                                             |                                                                                                                                  
CCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAA  >  NC_002947/1178751‑1179022

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