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.A93.1% 86.4 / 0.4 29intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (10/17);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGC  >  NC_002947/278608‑278864
                                                                                                                                       |                                                                                                                           
tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                <  6:93395/148‑1 (MQ=255)
tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                <  7:100191/148‑1 (MQ=255)
   tCGCGACAGACGTGTGATGGTCAGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                             >  4:82985/1‑148 (MQ=255)
      cgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTAcacaca                                                                                                          <  6:23719/148‑1 (MQ=255)
       gACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACAt                                                                                                         <  2:48210/148‑1 (MQ=255)
          aGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGcc                                                                                                      >  5:72026/1‑148 (MQ=255)
            aCGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACACGCTCCTACACACAAGCCt                                                                                                     >  2:127544/1‑147 (MQ=255)
                  gCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                               <  1:147670/147‑1 (MQ=255)
                            tggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGt                                                                                    <  4:119583/148‑1 (MQ=255)
                             ggtccggtATATATACCCTAATGATTTATAAAATTTTTAAGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATATGTGCACTGTa                                                                                   <  8:136416/148‑1 (MQ=255)
                                  ggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTATGAGc                                                                              <  7:11080/148‑1 (MQ=255)
                                    tatatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCgg                                                                            <  3:82985/148‑1 (MQ=255)
                                         ataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGttt                                                                        <  7:112108/147‑1 (MQ=255)
                                          taCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCGCTCCTACACACATGCCTCAAGATAGGTGCACTGGAGGAGCGGGTTTAc                                                                      >  1:144221/1‑148 (MQ=255)
                                           aCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAcc                                                                     <  1:39923/148‑1 (MQ=255)
                                              cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgc                                                                  <  7:106901/148‑1 (MQ=255)
                                                  tGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTGCACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGaaa                                                              >  5:148681/1‑148 (MQ=255)
                                                       tataAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                         >  7:137406/1‑148 (MQ=255)
                                                         taAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                        >  6:145150/1‑147 (MQ=255)
                                                               ttttATGAAATAACTTTTTGTAATATTTTTATACGGGGTATTTGTTACCTAGGTGCCGGGCTCTTCGCGGGTAAACGCGCTCCTACACACCTCCCTCAAGATCTGTGCACGGTAGCTGCGGGTTTACCCGCGAAAGGGCAGGAACAgg                                                 <  6:91819/148‑1 (MQ=255)
                                                                 ttattaATTAATTTTTTGTATTTTTTTTTTATGGGGTTTTTTTTACTGTGGGGCGGGGCTCTTCGGGGTTAAACCCGCTCATACACACTTGCCTCAATTTGTGTCCACTTTTGGTGGGGTTTTACCCGCGAAAGGGCCGGAACAGGt                                                <  2:170586/147‑1 (MQ=255)
                                                                 ttATGAACTAACTTTTTGTAATATTTATATACGGGGTATTTGTTACTGTGGTGCCGGGCTGTTCGGGGTTAAACCCGCTCCTACACACTTGCCTCAAGATCTGTGCACTGTTGGCGGGGTTTTACCCGCGAAAGGGCCGGAACAGGt                                                <  8:71584/147‑1 (MQ=255)
                                                                 ttATGAACTAACTTTTTGTAATATTTATAGACGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCCGCTCCTACACACATCCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGt                                                <  8:120238/147‑1 (MQ=39)
                                                                     taactaattTTTTGTAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTAc                                            <  6:56451/146‑1 (MQ=255)
                                                                     gAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCACCTAAACAAATGCCTCAAGATATGTGCACTTTAGGAGCGGGTTTACCCGCGATAGGGCCGGAACAGGTTTTCt                                           >  8:71010/1‑148 (MQ=255)
                                                                             tttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAAACCCCTCCGACACACATGCCTCAAGATATGTGGACTGTAGGAGCGGGGTTTCCCGCGAAAGGGCCGGCACAGGGTTACTCGGAGGcc                                   >  1:161720/1‑148 (MQ=255)
                                                                                             aGACGGGGTTTTTTTTACAGCGGGGCCGGGCTCTTCGGGGGAAACCCCGCGCCTACACACATGCCTCAAGAGCTGTGCATTTTCGGCGGGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAAc                    <  6:22435/147‑1 (MQ=255)
                                                                                                             aCAGCGGTGCGGGGCTCTTCGGGGGTAAACCCGCTCCTACACACTTGCCTCAAGATCTGTGCACTGTAGGAGGGGTTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGcc    <  2:161720/147‑1 (MQ=255)
                                                                                                              cAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCACCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgc  >  2:57105/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                           
TGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGC  >  NC_002947/278608‑278864

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