Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947699,7220AG41.5% 32.0 / 57.7 53intergenic (+115/‑129)PP_t06/PP_23SDtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base G (11/11);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTT  >  NC_002947/699576‑699803
                                                                                                                                                   |                                                                                 
ggTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGAc                                                                                  <  7:47407/148‑1 (MQ=37)
   cAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGa                                                                               >  4:134932/1‑148 (MQ=37)
         ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa                                                                         <  3:59229/148‑1 (MQ=37)
            gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGt                                                                       <  8:27432/147‑1 (MQ=37)
             aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGttt                                                                     <  4:149301/148‑1 (MQ=37)
             aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGtt                                                                      >  1:153569/1‑147 (MQ=37)
               cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa                                                                   <  4:69361/148‑1 (MQ=37)
                 cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAct                                                                 <  4:97857/148‑1 (MQ=37)
                       gCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt                                                           >  1:133964/1‑148 (MQ=37)
                         tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtgg                                                         >  4:37145/1‑148 (MQ=37)
                         tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg                                                         >  2:107779/1‑148 (MQ=37)
                          ccaccaCTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGCTTTCGGGCTTTTGTCAGATCGTTCTGTAAAAAGTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCCCTGCTGGTGGa                                                        >  8:61781/1‑148 (MQ=11)
                             ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa                                                     <  8:164217/148‑1 (MQ=37)
                                 tttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGc                                                  >  5:96518/1‑147 (MQ=37)
                                   tCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                               >  2:3027/1‑148 (MQ=37)
                                     gCCGTACGCAGTTACTTGTCCGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGGCAGATCGTTCTTTAAAAATTCGGATATGTGATCGATATAGGCTGAACCCCAGTTTCACTGCTGGTGGATCAGGCTAAg                                              <  4:118935/147‑1 (MQ=11)
                                     gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                              <  2:82049/147‑1 (MQ=37)
                                     gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                              <  3:100113/147‑1 (MQ=37)
                                           cGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCTAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATTTAGGCTGAACACCAGTTTCACTGCTGATGGATCAGGCTAAGGTAAAAt                                       >  8:70963/1‑148 (MQ=255)
                                            gCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGTACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAtt                                      >  1:161898/1‑148 (MQ=255)
                                            gCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAtt                                      >  1:120202/1‑148 (MQ=37)
                                             cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAttt                                     <  2:179003/148‑1 (MQ=37)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCTGATCGTTCATTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  4:83948/148‑1 (MQ=255)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  1:73395/148‑1 (MQ=37)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  3:2888/148‑1 (MQ=37)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  5:115215/148‑1 (MQ=37)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  1:107779/148‑1 (MQ=37)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  2:109975/148‑1 (MQ=37)
                                                ttACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtg                                    >  8:63408/1‑146 (MQ=35)
                                                   cTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAg                                <  3:119760/147‑1 (MQ=34)
                                                    ttGTGAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGtt                              >  6:110571/1‑148 (MQ=21)
                                                    ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGtt                              <  4:86974/148‑1 (MQ=21)
                                                    ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGtt                              >  2:22647/1‑148 (MQ=34)
                                                    ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGtt                              >  1:152070/1‑148 (MQ=34)
                                                      gTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                            >  2:95911/1‑148 (MQ=21)
                                                       tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTc                             <  5:36069/146‑1 (MQ=34)
                                                              ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                    <  3:163569/148‑1 (MQ=21)
                                                                 gAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                  <  2:33895/147‑1 (MQ=34)
                                                                  aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag                 <  5:20715/147‑1 (MQ=21)
                                                                  aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag                 >  1:167571/1‑147 (MQ=34)
                                                                    aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa              >  4:16500/1‑148 (MQ=21)
                                                                    aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa              <  1:10708/148‑1 (MQ=21)
                                                                    aTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa              <  8:28957/148‑1 (MQ=34)
                                                                     tcaaCATTCGAGCT‑GAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc             >  6:30769/3‑148 (MQ=2)
                                                                      gAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                                                                >  6:167146/1‑80 (MQ=39)
                                                                      gAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc             >  7:79292/1‑147 (MQ=34)
                                                                      gAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc             >  3:92021/1‑147 (MQ=34)
                                                                       aaCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGa                                                                               <  5:167146/80‑1 (MQ=39)
                                                                         cATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTAc         <  6:56018/148‑1 (MQ=21)
                                                                         cATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTAc         <  7:19851/148‑1 (MQ=34)
                                                                               aGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtt    >  1:165074/1‑148 (MQ=32)
                                                                               aGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAt     <  8:15676/147‑1 (MQ=32)
                                                                                 cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtttt  >  6:161200/1‑148 (MQ=11)
                                                                                                                                                   |                                                                                 
GGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTT  >  NC_002947/699576‑699803

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