Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947181,8120AG28.2% 51.5 / 23.2 40intergenic (+26/‑80)PP_r01/PP_5SB5S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (12/16);  new base G (5/6);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAA  >  NC_002947/181672‑181950
                                                                                                                                            |                                                                                                                                           
gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGGCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAgg                                                                                                                                       <  1:35112/147‑1 (MQ=17)
  cgacCATAGAGTATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAg                                                                                                                                    <  3:202033/148‑1 (MQ=17)
   gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                  <  5:22421/149‑1 (MQ=255)
           agCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCAACGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTc                                                                                                                         <  1:309328/149‑1 (MQ=17)
            gCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAGGGGT‑TTTgcga                                                                                                                        <  8:111393/149‑4 (MQ=11)
              aTTGGAACCCCCAGATCCCATCCCCAAATCAGTAGTGAAACGATGCATAGCCGAAGGGAGTGTGGGGTTTGCCCAAGTGAGAGTAGGTGATAGTCAAGATTAATTTAGCAAAAACCCAAAAAGCACATGCAGGTAGGGGt‑tttttgtt                                                                                                                       >  4:80085/1‑143 (MQ=2)
              aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATAGGGGT‑TTTgcgaaa                                                                                                                      <  8:270706/149‑6 (MQ=11)
              aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTACCTGCATGTGCAGATAGGGGT‑TTTgcgaaa                                                                                                                      <  8:232375/149‑6 (MQ=11)
                 ggAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTaa                                                                                                                    <  3:80085/148‑1 (MQ=32)
                    accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAg                                                                                                                <  6:180641/149‑1 (MQ=17)
                    accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                  >  6:246354/1‑132 (MQ=11)
                    accaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCACAAAACCCCTACCTGCATGTGCAGATAggg                                                                                                                                  <  5:246354/132‑1 (MQ=11)
                    accaccTGATCCCATCACGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAg                                                                                                                <  7:1435/149‑1 (MQ=17)
                      caccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAgg                                                                                                               >  4:267077/1‑148 (MQ=32)
                                                    aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAgg                                                                                 >  7:134045/1‑148 (MQ=255)
                                                       gATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACg                                                                             >  6:147134/1‑149 (MQ=255)
                                                            aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAggg                                                                                                                                  <  1:205403/92‑1 (MQ=33)
                                                            aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCa                                                                        <  2:218761/149‑1 (MQ=255)
                                                                 cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcac                                                                    >  3:186540/1‑148 (MQ=255)
                                                                       tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATACAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAt                                                             >  7:237499/1‑149 (MQ=255)
                                                                          tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTc                                                          <  3:267077/149‑1 (MQ=255)
                                                                            tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                        <  1:236497/149‑1 (MQ=255)
                                                                            tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCtt                                                        >  1:308208/1‑149 (MQ=255)
                                                                                    ccccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                >  6:274192/1‑149 (MQ=255)
                                                                                     cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa                                                >  5:24949/1‑148 (MQ=255)
                                                                                        aTGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATaga                                            <  5:147134/149‑1 (MQ=255)
                                                                                                gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATACAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGa                                    <  8:237499/149‑1 (MQ=255)
                                                                                                   ggTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTACACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTACTAGGGAAGGGCCAGCACACAGAATTTATTGACGACCATAGAGCATTGGAac                                  >  1:55653/1‑148 (MQ=255)
                                                                                                             aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatccc                       <  5:274192/149‑1 (MQ=255)
                                                                                                              aGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacac                                                                  >  2:101000/1‑105 (MQ=255)
                                                                                                              aGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGcacac                                                                  <  1:101000/105‑1 (MQ=255)
                                                                                                                     tttCGCAAAACCCCTATCTGCACTTGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa               <  6:24949/149‑1 (MQ=255)
                                                                                                                       agacAAAACCCCTACCTGCATGTGCAGATAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACt             >  4:29368/4‑149 (MQ=37)
                                                                                                                          cAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag          >  3:201824/1‑149 (MQ=255)
                                                                                                                           aaaaCCCCTACCTGCATGTGCAGATAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                          <  6:305080/132‑1 (MQ=255)
                                                                                                                           aaaaCCCCTACCTGCATGTGCAGATAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat                          >  5:305080/1‑132 (MQ=255)
                                                                                                                              aCCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGt      <  4:186540/149‑1 (MQ=255)
                                                                                                                                cccTACCTGCATGTGCAGATAGGGGTATTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGt      >  1:3804/1‑148 (MQ=37)
                                                                                                                                cccTACCTGCATGTGCAGATAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTg     >  3:114867/1‑148 (MQ=255)
                                                                                                                                  cTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGaaa  <  8:134045/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTAGGGGT‑TTTGTCTTTAAGTAGGAACTACAGAGATTCGCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAA  >  NC_002947/181672‑181950

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