Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,606,450 A→T 17.6% intergenic (‑158/+53) PP_1407 ← / ← PP_t64 sulfate transporter/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,606,4500AT17.6% 78.9 / 11.2 34intergenic (‑158/+53)PP_1407/PP_t64sulfate transporter/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base T (4/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACT  >  NC_002947/1606301‑1606590
                                                                                                                                                     |                                                                                                                                            
cTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCaa                                                                                                                                              >  1:22366/1‑150 (MQ=255)
 tttCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAACAAAGAAACAAGGAGCCAAg                                                                                                                                             >  2:169879/1‑150 (MQ=255)
           tttCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTcc                                                                                                                                        >  2:738582/1‑145 (MQ=255)
           tttCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTcc                                                                                                                                        <  1:738582/145‑1 (MQ=255)
           tttCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTg                                                                                                                                   >  2:364449/1‑150 (MQ=255)
                   aTCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGAt                                                                                                                           <  1:735826/150‑1 (MQ=255)
                             aatatCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAATGCAGAATAAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTg                                                                                                                 <  2:192871/149‑1 (MQ=255)
                                 tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAgg                                                                                                                                            <  1:240618/119‑1 (MQ=255)
                                 tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAgg                                                                                                                                            >  2:240618/1‑119 (MQ=255)
                                 tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGAc                                                                                                                        >  1:197882/1‑139 (MQ=255)
                                 tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGAc                                                                                                                        <  2:197882/139‑1 (MQ=255)
                                   cGCCACAAGCCGTGTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGGAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGTACGc                                                                                                           <  2:189928/150‑1 (MQ=255)
                                                    cTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATTGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAtg                                                                                          <  1:435094/150‑1 (MQ=255)
                                                             gcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                                 >  2:504203/1‑150 (MQ=255)
                                                                 aGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTaga                                                                             <  1:169879/150‑1 (MQ=255)
                                                                  ggCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTagag                                                                            <  2:22366/150‑1 (MQ=255)
                                                                                     cGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAgg                                                                                                                                            <  1:514835/67‑1 (MQ=255)
                                                                                     cGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAgg                                                                                                                                            >  2:514835/1‑67 (MQ=255)
                                                                                            gTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTag                                                                              <  1:221412/122‑1 (MQ=255)
                                                                                            gTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTag                                                                              >  2:221412/1‑122 (MQ=255)
                                                                                            gTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTg                                                  <  2:770124/150‑1 (MQ=255)
                                                                                               ggTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTg                                                                                                                               >  1:621339/1‑70 (MQ=255)
                                                                                               ggTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTg                                                                                                                               <  2:621339/70‑1 (MQ=255)
                                                                                                                   gTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTaa                           >  2:427764/1‑150 (MQ=255)
                                                                                                                                  aaaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCCGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCt            <  1:571873/150‑1 (MQ=255)
                                                                                                                                    aaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAgg                                                                           <  2:313288/85‑1 (MQ=255)
                                                                                                                                    aaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAgg                                                                           >  1:313288/1‑85 (MQ=255)
                                                                                                                                        tcaaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCg                                                                                   >  2:214077/3‑73 (MQ=38)
                                                                                                                                        tcaaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCg                                                                                   <  1:214077/71‑1 (MQ=38)
                                                                                                                                          aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGtttt                                   >  2:443241/1‑119 (MQ=255)
                                                                                                                                          aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGtttt                                   <  1:443241/119‑1 (MQ=255)
                                                                                                                                           aaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGAc   >  1:617401/1‑150 (MQ=255)
                                                                                                                                           aaaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGAc   <  2:758872/150‑1 (MQ=255)
                                                                                                                                            aaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACt  >  2:269880/1‑150 (MQ=255)
                                                                                                                                            aaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGGAGAGGGAATAGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACt  >  1:181247/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
CTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACT  >  NC_002947/1606301‑1606590

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