Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,9470CT45.2% 42.5 / 19.9 42E276K (GAA→AAA) ygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base C (4/19);  new base T (19/0);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 2.08e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTA  >  NC_000913/2878802‑2879087
                                                                                                                                                 |                                                                                                                                            
tGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTtct                                                                                                                                             >  2:801273/1‑147 (MQ=255)
          gcgcGGAGGTACATTTTCAGTGACCACGACCAACAGACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCa                                                                                                                                <  1:615342/150‑1 (MQ=255)
           cgcgGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCCAg                                                                                                                                  >  1:804676/1‑147 (MQ=255)
                         ttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCACCGGCAAGGCCGTCCCCTATAAGCg                                                                                                                 >  1:721972/1‑150 (MQ=255)
                         ttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCgg                                                                                                                >  2:229508/1‑151 (MQ=255)
                         ttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCg                                                                                                                 <  1:862889/150‑1 (MQ=255)
                           cAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGa                                                                                                               <  2:141509/150‑1 (MQ=255)
                           cAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATCATTTTCTCCAAAGGCacagacgtc                                                                                                                            >  1:705497/1‑129 (MQ=255)
                                aCCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCCAACACGTCCTCTATAAGCGGGATCaa                                                                                                          >  2:938790/1‑150 (MQ=255)
                                  cACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAAttt                                                                                                       <  1:117301/151‑1 (MQ=255)
                                     gcccAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTTTTTTTTTTCAGCGGCGAACATGTGCTCTCTAAACGGaaaaaa                                                                                                          >  1:932769/3‑141 (MQ=255)
                                         aaCATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCCAAGAAGTCCTCTATAAGCGGAATCAATTTTGGTaa                                                                                                 >  1:486101/1‑150 (MQ=255)
                                                    ttCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTACAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtactac                                                                                      <  1:628813/150‑1 (MQ=255)
                                                     tCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTactact                                                                                     <  2:551328/150‑1 (MQ=255)
                                                        gCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGa                                                                                 <  1:881924/151‑1 (MQ=255)
                                                           aCTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGaa                                                                                <  1:932231/149‑1 (MQ=255)
                                                           aCTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTATTTCTCCAGCGGGAAACACGTCCTTTTTAAGCGGAATCCATTTTGGTAATGTTTTACTACTGcccaa                                                                               >  2:487772/1‑146 (MQ=255)
                                                             tCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCAGCGGCAAACAAGTCCTCTATAAACGGAATCAATTTGGCTAATGTTTTACTACTGCGaaaaa                                                                             >  2:231885/1‑150 (MQ=255)
                                                              ccGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTTCCCACGGGCAGCACGGCCTCCATAAGCGGGATCAATTTTGCTAATGTTTTTCTTCTTCGgaaaa                                                                             >  2:435478/1‑149 (MQ=255)
                                                               cGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAatat                                                                          <  1:1018632/151‑1 (MQ=255)
                                                                   ggCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCt                                                                      <  1:631471/151‑1 (MQ=255)
                                                                     ccTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCAAGCACGTCCCCCATTAACGGAATAAATTTTGCTAATGTTTTACTAACGCGAAAAAAATtcc                                                                       >  1:310629/1‑148 (MQ=255)
                                                                      cTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGGACGTCCTCTATAAACGGAAACAATTTGGCTAAAGTTTTTCTACTGCGGAAAATAAccc                                                                       >  1:153789/1‑147 (MQ=255)
                                                                        gCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCa                                                                   <  2:887813/149‑1 (MQ=255)
                                                                        gCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCAc                                                                  <  1:202486/150‑1 (MQ=255)
                                                                          aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAACGGCAAGCACGTCATCTATAAGCGGAATCAATTTTGCTAATGTTTTACTACTGCGAAAAATATCCCTTCACGcc                                                               >  1:740321/1‑151 (MQ=255)
                                                                             tcCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTTCCAGCGGCAAGCACCTCCTCTATAAGCGGAATAAATTTTGCTAATGTTTTAATACTGCGAAAAATATCCCTGCCCGCCaa                                                             >  2:523689/1‑150 (MQ=255)
                                                                                   tGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGGAAGCACGTTCTCTATAAACGGGAACAATTTTGGTAAAGTTTTACTACTGCGGAAAAAATTCCTGCACCCCAAACGGAc                                                       >  2:930030/1‑150 (MQ=255)
                                                                                          ggAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCCAGCACCTCCCCTATAAACGGGATCAAATTGGCTAAAGTTTTACTACCGCGGAAAATATCCCTCCACCCCAAACGGACCTccccg                                                >  2:612904/1‑148 (MQ=255)
                                                                                          ggAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCCGCGGCCAGCACGTCCTCTCTAAGCGGGATCCATTTTGCTAATGTTTTATTACTGCGAAAAATATCCCTGCCCGCCAAACCGGCTTccccg                                                >  1:154022/1‑148 (MQ=255)
                                                                                            aaGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAACGGAAACAAATTTGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTTCCGGTcc                                             >  2:458470/1‑151 (MQ=255)
                                                                                            aaGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTcc                                             <  1:879823/151‑1 (MQ=255)
                                                                                             aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCAGCGGGAAGGACGCGCCCTTTAAAGGGGAACAAATTTGCTAATGTTTTACTACTGCGAAAAAAATTCCCGCACGCCCAACGGGCTTCCCCGTcc                                             >  1:629116/1‑150 (MQ=255)
                                                                                                ggAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGGCCTCTTTAAGCGGAATCCATTTGGGTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCCGTCCgg                                           >  2:733560/1‑149 (MQ=255)
                                                                                                     ggCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCAcc                                     <  1:243879/150‑1 (MQ=255)
                                                                                                      gCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAg                                   <  2:686806/151‑1 (MQ=255)
                                                                                                           aGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAggggt                               >  2:958434/1‑148 (MQ=255)
                                                                                                                     tCTTCAGGTGGGGCCGGCGGTTTTATTTTTCCAGCGGCGAGCACCTCCTCTATAAACGGGATCAATTTGGCTAATGTTTTACTACTGCGGAAAATATCCCTGCACGCCAAAAGGACATCCCCGTCCGGGTCACCAGGGTTACGAcgcg                       >  2:272119/1‑148 (MQ=255)
                                                                                                                        tCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGGAAGCACGTCCTCTATAAGCGGAATCAATTTTGCTTATGTTTTTCTACTGCGAAAAAAAATCCCGCACGCCAAAAGGACTTCCCGGTCCGGCCCACCAGGGGTTAGACGCGCTAt                   >  1:402696/1‑149 (MQ=255)
                                                                                                                        tCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTAt                   <  1:877731/149‑1 (MQ=255)
                                                                                                                            gTGGGGTAGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCACGGTCCGGCTCACCAGGGTTACGACGCGCTATCGCaaa             <  1:366465/151‑1 (MQ=255)
                                                                                                                              ggggCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCaaaa            <  1:694504/150‑1 (MQ=255)
                                                                                                                                        gTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGGCTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTa  <  2:271782/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTA  >  NC_000913/2878802‑2879087

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑