Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 278,742:1 +A 100% intergenic (‑237/+77) PP_0224 ← / ← sctC DszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947278,7421.A97.6% 131.3 / ‑2.9 41intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (21/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278608‑278885
                                                                                                                                       |                                                                                                                                                
tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                     >  6:47189/1‑148 (MQ=255)
       gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATg                                                                                                                             >  6:319857/1‑149 (MQ=255)
       gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATg                                                                                                                             >  6:170718/1‑149 (MQ=255)
          aGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGTTTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACACGCTCCTAAACACATGcc                                                                                                                           >  5:265360/1‑148 (MQ=255)
             cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                                       <  5:319857/149‑1 (MQ=255)
              gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                              >  2:129829/1‑141 (MQ=255)
              gtgtGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                              <  1:129829/141‑1 (MQ=255)
                gtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                              <  6:112520/123‑3 (MQ=255)
                gtgCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                              >  5:112520/1‑121 (MQ=255)
                        ccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCAc                                                                                                            >  5:140887/1‑149 (MQ=255)
                         cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacacac                                                                                                                                <  1:275019/128‑1 (MQ=255)
                         cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacacac                                                                                                                                >  2:275019/1‑128 (MQ=255)
                            tggtccggtATATATACCCTAATGATCTATAAAATTATTCTCAACTAACTTTTTGGAATATTTATAGAAGGGGGATGTGTTACAGCGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTCTGCACTGTa                                                                                                        <  6:140887/149‑1 (MQ=255)
                                 cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGc                                                                                                   >  1:232985/1‑149 (MQ=255)
                                   gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                                                 <  5:47189/149‑1 (MQ=255)
                                      tatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGt                                                                                               <  5:102615/148‑1 (MQ=255)
                                            cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcg                                                                                        <  5:170718/149‑1 (MQ=255)
                                             ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                                  <  6:305611/122‑1 (MQ=255)
                                             ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                                  >  5:305611/1‑122 (MQ=255)
                                             ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                       >  8:323906/1‑149 (MQ=255)
                                             ccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACACGCTCCTAAAAACATGCCTAAAGATATGTGCACTGTAGGAGAGGGTTTaacagc                                                                                       >  7:143109/1‑144 (MQ=255)
                                               tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACt                                                                                                           >  5:286668/1‑127 (MQ=255)
                                               tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACt                                                                                                           <  6:286668/127‑1 (MQ=255)
                                                    aTTAATAAAATTTTTATTAACTAACTTTTTTTAATATTTATAGAAGGGGGATTTGTTACACAGGTGCCGGGCTCTTCGCGGTTAAACCCGCTCCTACACACATGCCTCAATATCTTTGCACTGTAGGAGGGGTTTTACCCGCGAAAgg                                                                                 <  8:143109/148‑1 (MQ=255)
                                                          aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGa                                                                           >  7:255006/1‑148 (MQ=255)
                                                                  tATGATTTAATTTTTTGAAATTTTTTTAGTTGGTGGTTTTTTTAGAGTGTGGCGGGGCTGTTCCGGTTTAACCCCCCTCCTACAAACTTCCATCATGTTCTGTGCACTTTAGGAGGGGTTTTACCCGCGAAAGGGCCGGAACAGGttt                                                                   <  2:196875/148‑1 (MQ=255)
                                                                    tGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCACTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACt                                                                <  6:9214/149‑1 (MQ=255)
                                                                          aactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATTCCTCAAGATCTGTGGACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                              >  6:278412/1‑129 (MQ=255)
                                                                          aactTTTTGGAATATTAATACAAGCGGGATTTGTAACAGAGGTCCCGGGCGCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGCAGCGGGTTTACCCGCGAAAGGGcc                                                                              <  5:278412/129‑1 (MQ=255)
                                                                          aactTTTTGCAATATTAATAGAAGGGGGATTTGTTACAGAGGTCCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATCCCTCAACATCTGTGCACTGTAGCAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAgg                                                          <  4:55505/149‑1 (MQ=255)
                                                                                      tatTAATAGAAGGGGTAGTGGTTCCAGAGGTCCCGGCCTCATCGCGGCTAAACCCGCTCCTCCCCACATGCCTCAAGATCTGTGCACTGTACTAGGGGGTTTACCCGCGAAAGGGCCGTAACAGGTTTACTCGTAGGCCGCAGTCCCGa                                              <  6:265360/149‑1 (MQ=255)
                                                                                                       tttGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                                    >  7:238819/1‑142 (MQ=255)
                                                                                                       tttGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                                    <  8:238819/142‑1 (MQ=255)
                                                                                                               agagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgcg                      >  8:240810/1‑148 (MQ=255)
                                                                                                                  ggTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACg                  >  4:237116/1‑149 (MQ=255)
                                                                                                                      ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCATCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgc                                     <  4:254720/126‑1 (MQ=255)
                                                                                                                      ccGGGCTCTTCGCGGGTAAACCCGCGCCTACGCACCTGCCTCAAGGTCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACCGGTTTACTCGGAGGCCGCAGTCCCGAGGAACATgc                                     >  3:254720/1‑126 (MQ=38)
                                                                                                                            tctTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttccttc        >  6:77880/1‑149 (MQ=255)
                                                                                                                                cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCgg      <  2:232985/147‑1 (MQ=255)
                                                                                                                                 gcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGgcgc   <  5:77880/149‑1 (MQ=255)
                                                                                                                                  cgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc                         >  6:253398/1‑126 (MQ=255)
                                                                                                                                  cgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc                         <  5:253398/126‑1 (MQ=255)
                                                                                                                                  cgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCt  >  6:129423/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                
TGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278608‑278885

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