Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 278,744 +A intergenic (‑239/+75) 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.A82.6% 57.5 / 8.1 23intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base A (9/10);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTG  >  NC_002947/278599‑278860
                                                                                                                                                 |                                                                                                                        
cGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                      <  6:98199/148‑3 (MQ=255)
                     aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTc                                                                                                 <  3:143639/148‑1 (MQ=255)
                             tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACCCATGCCTCAAGATCtg                                                                                         >  6:28116/1‑148 (MQ=255)
                             tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtg                                                                                         >  5:86349/1‑148 (MQ=255)
                                tccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                        <  5:104298/146‑1 (MQ=255)
                                tccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGCGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                        <  2:13443/146‑1 (MQ=255)
                                   ggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACt                                                                                   >  3:138159/1‑148 (MQ=255)
                                     tggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACt                                                                                   >  5:84076/1‑146 (MQ=255)
                                     tggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGt                                                                                 <  7:6496/148‑1 (MQ=255)
                                        tccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                >  8:122377/1‑146 (MQ=255)
                                          cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCGCAAGATCTGTGCACTGTAGGAg                                                                            >  1:146028/1‑148 (MQ=255)
                                              atatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCggg                                                                       <  2:118360/148‑1 (MQ=255)
                                                    aCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTTAAAACGCTACTACACACATGCCTCAAGATCTGTTCACTGTTGGGGCGGGTTTAcc                                                                 >  3:4368/1‑148 (MQ=255)
                                                    aCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTAACCACATGCCTCAAGATCTGGGCCCTGTAGGAGCGGGTTTAcc                                                                 >  8:41947/1‑148 (MQ=255)
                                                       cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                              >  3:39331/1‑148 (MQ=255)
                                                         aaTGCTTTATAAAATTTTTATGATGTAACTTTTTTGAATTTTTTTATAAGGGGGATTTTTTACTGAGGTGCCGGGGTCTTCGCTGGTAACCCCGCTCCTAAACACATGCCTCAAGATCTGTGCAGTGTAGGAGCGGGTTTACCCGCGa                                                            <  4:99320/148‑1 (MQ=255)
                                                              tttATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGATGGGTGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAgg                                                        <  8:45042/147‑1 (MQ=255)
                                                                                actaactTTTTGTAATATTTATATACGGGGGATTTGTTACAGAGGGGCGGGGCTCTTCGCGGTTAAACCCGCTCCTACACACATGCCTCAAGATCTGGGCACTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTc                                      <  6:170015/147‑1 (MQ=255)
                                                                                                atTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGa                     <  6:86349/148‑1 (MQ=255)
                                                                                                atTTATAGAAGGGGGATTTGTTACAGAGGGGCCGGGCGCTTCGCGGGTAACCCCGCTCCTACACACATGCCTCACCAGCGGTCCACTGTCGGCGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGa                     <  3:34537/148‑1 (MQ=37)
                                                                                                  ttATAGACGGGGGATTTGTTCCAGCGGGGCCGGGCTCTGCGCGGGTAAACCCGCTCCTACACACAGGCCTCAAGATCTATGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAg                    <  7:93118/147‑1 (MQ=38)
                                                                                                          gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCGCCGCCCCCCCGGCCTCAAGCTCGGGGCACTGTAGGCGCGGGTGTACCCGCGAAAGGGCCGGGACAGGGGTAATCGGGGGCCGCCCTCCCGAGGaccgtgg            >  6:11651/1‑141 (MQ=255)
                                                                                                                 ttGTTACTGCGGGGCGGGGCTCTTCGGGGTTAAACCCGCTCCTACAAACATGCCTAAATTTCTGTGCACTGTTGGTGGGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGaa     <  4:4368/147‑1 (MQ=255)
                                                                                                                   gTTACAGAGGTGCCGGGCTCTTCGCGGGTTAAACCCCTCCAACACACAAGCCGCAAGATATGGGGCCAGTAGGCGCGGGGTTTACCCCGAAAGGGCCCGCACACGGTTAAAAGGCGGCCGCAGTGCCCCGGCACATGGGCaacaacag  >  6:101120/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                        
CGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTG  >  NC_002947/278599‑278860

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: