Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947278,7421.A70.0% 68.0 / 28.8 30intergenic (‑237/+77)PP_0224/sctCDszC family monooxygenase/sulfur compound ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (3/6);  new base A (11/10);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑T‑A‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTC  >  NC_002947/278602‑278875
                                                                                                                                              |                                                                                                                                       
tACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑T‑A‑A‑CCCGCt                                                                                                                               >  7:306874/1‑149 (MQ=255)
  cAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCta                                                                                                                              >  3:170025/1‑148 (MQ=255)
    gCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTcc                                                                                                                             >  5:237439/1‑148 (MQ=255)
    gCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑T‑A‑A‑CCCGCTCCTa                                                                                                                           <  8:306874/149‑1 (MQ=255)
     cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑T‑A‑A‑acc                                                                                                                                  >  3:270492/1‑141 (MQ=255)
     cTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑T‑A‑A‑acc                                                                                                                                  <  4:270492/141‑1 (MQ=255)
      tGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCGTac                                                                                                                          >  4:239963/1‑149 (MQ=255)
                 gACGTGTGCTGGTCCGGTGGTCCTGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTc                                                                                                               <  7:766/149‑1 (MQ=255)
                   cGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGCTTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑T‑A‑ACCCCGCTCCTACACACATGCCTCaa                                                                                                             <  6:327340/149‑1 (MQ=255)
                         cTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑ACCGCTGCTTCACACATGCCTCAAGGTCtg                                                                                                       >  3:41925/1‑149 (MQ=255)
                          tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCtgt                                                                                                      <  2:221387/149‑1 (MQ=255)
                                ggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTg                                                                                                <  7:128584/149‑1 (MQ=255)
                                       cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGc                                                                                         >  2:95570/1‑149 (MQ=255)
                                             atataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGtt                                                                                    >  6:186020/1‑148 (MQ=255)
                                                taCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAc                                                                                 <  4:138183/148‑1 (MQ=255)
                                                    cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGtt                                                                                    >  1:201118/1‑141 (MQ=255)
                                                    cTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGtt                                                                                    <  2:201118/141‑1 (MQ=255)
                                                     tAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTT‑A‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCcgcg                                                                            >  1:141474/1‑148 (MQ=255)
                                                          aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAggg                                                                      <  5:186020/149‑1 (MQ=255)
                                                               taAAATTTTTATGAACTAACTTTTTGCAATATTTATAGACGGGGTATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTCCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCgg                                                                  <  6:237439/148‑1 (MQ=255)
                                                                           taactaactTTTTGCAATTTTTATACACGGGGTATTTGTTACCGCGGGGCCGGGCTCTTCGCGGG‑T‑A‑ACCCCGCTCCAACACACTTGCCTCAAGATCTGGGCACTGTCGGCGGGGTTTAACCCGCGAAAGGGCCGGAACAGGTTTACt                                                      <  2:141474/147‑1 (MQ=255)
                                                                                 cttttttGTAATTTTTATCGAAGGGGCTTTTGTCTCTGGGGGGCGGGGCCCTTCGCGGT‑T‑ACC‑CCCGCTCCACCACACATCCCCAAAGAGCTCTGCTGGGTCGGCGCGGGTTAACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAgg                                                <  4:41925/146‑1 (MQ=255)
                                                                                              tttATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAg                                   <  1:95570/148‑1 (MQ=255)
                                                                                                 ataGAAGGGGTATTTTTTACCGCGCGCCCGGGCTCTTCGCGGG‑T‑A‑ACCCCGCTCCAACACACAGCCCTCAAGATCTGCGCACTGTAGCCGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGaa                                <  6:25193/148‑1 (MQ=25)
                                                                                                  tatatgggggATTTTTTCCTGAGGTGGCGGTCTTTTGGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCGTCATGATTTGTGAAGTTTCGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAAc                               <  8:163303/143‑1 (MQ=255)
                                                                                                  taGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCGACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACACCGAGGCCGCAGTCCCGAGGa                                 >  7:163303/1‑146 (MQ=39)
                                                                                                       gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACCCATGCCTCAAGATCTGTGGACTGTATGAGCGGGTTTACCCGCGCAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTgcg                          >  7:49847/1‑148 (MQ=255)
                                                                                                       gggggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGGAGTGCCGAGGAACTTtcg                          >  7:56081/1‑148 (MQ=255)
                                                                                                                        ggTCCCGGGCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACg        <  3:71907/149‑1 (MQ=255)
                                                                                                                              gggCTCTTCGCGGG‑TAA‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACAGTAGGAGTGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttc  >  1:308259/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                       
TACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGG‑T‑A‑A‑CCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTC  >  NC_002947/278602‑278875

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: ‑