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.A94.4% 109.3 / 0.3 36intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (14/20);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278596‑278885
                                                                                                                                                       |                                                                                                                                                 
tGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTa                                                                                                                                                  >  2:117888/1‑148 (MQ=255)
        atcttGATCGCGACAGAAGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTc                                                                                                                                         <  7:349928/145‑1 (MQ=255)
           cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTa                                                                                                                                      <  2:259044/149‑1 (MQ=255)
             gATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTc                                                                                                                                         >  8:349928/1‑144 (MQ=255)
              aTCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATATAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTacac                                                                                                                                   <  2:21202/149‑1 (MQ=255)
                 gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTAcacaca                                                                                                                                >  4:196905/1‑149 (MQ=255)
                  cgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                              <  5:159416/149‑1 (MQ=255)
                  cgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                                              <  2:131185/149‑1 (MQ=255)
                        aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGcc                                                                                                                           >  4:58591/1‑146 (MQ=255)
                        aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGcc                                                                                                                           <  3:58591/146‑1 (MQ=255)
                         cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                                       <  6:168715/149‑1 (MQ=255)
                                  gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                                <  8:229311/147‑1 (MQ=255)
                                      ggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                          <  4:273363/149‑1 (MQ=255)
                                       gtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTaca                                                                                                                                    <  5:267901/123‑1 (MQ=255)
                                       gtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTaca                                                                                                                                    >  6:267901/1‑123 (MQ=255)
                                        tggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                                        <  6:222496/149‑1 (MQ=255)
                                            ccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAAACCGCTACTACAAACATGCCTCAAGATATGTGCACTGTAGGAg                                                                                                    >  1:373776/1‑149 (MQ=255)
                                               gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCt                                                                                                                                       >  8:288961/1‑112 (MQ=255)
                                               gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCt                                                                                                                                       <  7:288961/112‑1 (MQ=255)
                                                    tataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGttt                                                                                             <  6:23948/148‑1 (MQ=255)
                                                        cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcg                                                                                        <  8:294380/149‑1 (MQ=255)
                                                           tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                                      >  2:103528/1‑148 (MQ=255)
                                                             aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCCTCCCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                                       <  1:117888/148‑1 (MQ=255)
                                                              tGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAAACCGCTCCTACACACATGCCGAAAGATATGTGCACTGTAGGAGCGGGTTTACCCGCGAAAg                                                                                  >  2:324833/1‑149 (MQ=255)
                                                               gATTTATAAAATTTTTATGAACTAACTTTTTGGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGAT‑‑CTTCGAGGGTAAACCCGCTCCTACACACATGCATCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAg                                                                                  >  1:40047/1‑149 (MQ=255)
                                                                    ataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAAAACGCTACTTAACAAATTCGGCAATATATGTGGACTGTATGAGGGGGTTTTCCCGCGAAAAGGCCgg                                                                            >  5:262656/1‑149 (MQ=255)
                                                                            tttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATTTCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAgg                                                                      <  1:103528/148‑1 (MQ=39)
                                                                                     taactTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACAAACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAg                                                           >  5:151859/1‑149 (MQ=255)
                                                                                           tttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                            <  2:41481/142‑1 (MQ=255)
                                                                                           tttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGa                                                            >  1:41481/1‑142 (MQ=255)
                                                                                           tttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGCGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCg                                                       <  1:102076/147‑1 (MQ=255)
                                                                                                           gAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAATCCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTg                                      >  5:89983/1‑149 (MQ=255)
                                                                                                            aaGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgc                                     <  6:151859/149‑1 (MQ=255)
                                                                                                                                gTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGc                          <  7:284905/140‑1 (MQ=255)
                                                                                                                                gTGCCGGGCT‑‑CTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGc                          >  8:284905/1‑140 (MQ=255)
                                                                                                                                                cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGCTCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGgcg    >  3:354988/1‑149 (MQ=255)
                                                                                                                                                  cgGGTAATCCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCt  <  6:89983/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
TGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCT‑‑CTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCT  >  NC_002947/278596‑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: