Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,9470CT42.9% 22.6 / 11.8 28E276K (GAA→AAA) cas1multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (3/13);  new base T (12/0);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.54e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCT  >  NC_000913/2878700‑2879182
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
ccGCCAGTCCATCTATTTGTTCCCAGAACATTTCACGAATTTTTGCGGATACATCACCTACATATTCCCCTGCCCGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGATGTACATTGTCAGTGACCACGACCAACATACTCATTTCAGCGACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTcc                                                                                                                                                                                                                                                <  2:323695/251‑1 (MQ=255)
          cgCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTGCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGGATCTTCAGGTGGGTCCGGCGGTTGTATTTCTCCAGCGGCAAGc                                                                                                                                                                                                                                      <  1:428279/250‑1 (MQ=255)
                    tCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCctct                                                                                                                                                                                                                             <  1:454458/250‑1 (MQ=255)
                                tCACGAATTTTTGCGGATACATCACCTACATAGACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGCGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTTCTCTATA‑AG‑CGGAATc                                                                                                                                                                                                               <  2:41459/250‑1 (MQ=255)
                                      aTTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTATTTTTTCAGCGGCAAGCACGGCCTCTTTA‑AG‑CGGGAACAAttttg                                                                                                                                                                                                        >  1:279413/1‑249 (MQ=255)
                                            gCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑At                                                                                                                                                                                                  <  2:400912/250‑1 (MQ=255)
                                                    aTCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACAAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTGTTTTTCCAGAGGCACGCACGTCCCCTATA‑AA‑CGGAATCAATTTGGCTC‑ATGtttttat                                                                                                                                                                                          >  1:188239/1‑247 (MQ=255)
                                                          tACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTGTTTTCCCACGGCAACCACGGCCTCTATA‑AG‑CGGGATCAATTTTGCTA‑ATGTTTTTATACTGCg                                                                                                                                                                                    >  1:377697/1‑250 (MQ=255)
                                                             atatACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTAGA‑AG‑CGGAATCAATTTGGCTA‑AGGTTTTACTACTGCGaaaa                                                                                                                                                                                <  2:334271/251‑1 (MQ=255)
                                                                                   aacaacCAGAAGGCTAATCTGCCTCGTAAACGCAGAGGTACATGTCCAGTGACCACGACCAACATACTCAATGCAGCTACTCCGATGGCCTGCATCTCCCAGGGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCTGCGGTTGTATTGCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CCGAATCAATTTGGCTA‑ATGTGTTACTACTGCGAAAAATATCCCTGCACTCCAAA‑CGGa                                                                                                                                                         <  2:279413/251‑1 (MQ=255)
                                                                                                            gTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCCAGCCCGTCCTCTATA‑AG‑CGGGATCAATTTGGCTC‑CTGTTTTACTTCTGCGGAAAATATCCCTGCAAGCCCAA‑CGGAC‑TT‑CCCGGGCCGGCTCACCAGGGtt                                                                                                                              >  1:83661/1‑251 (MQ=255)
                                                                                                                   cgGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGGAAGCACCGCCTCTATT‑AGACGGA‑TCAATTTTGCTA‑ATGTTTTACTTCTGCGAAAAATATCCCTTCACGCCAAA‑CGGAC‑CT‑TCCGGTCCGGCTCACCAGGGTTACGAcgc                                                                                                                       >  1:77293/1‑251 (MQ=255)
                                                                                                                   cgGAGGAACATTTTAAGTGACCCCGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCACGGCAAGCACGGCCTCTATACAA‑CG‑AAAAAATTTGGCTA‑ATGTTTTTTTAATTCGAAAAATATACCTCCACGCCCAA‑CAGAC‑TT‑CCCCGTcccggctaccccggctttgcccc                                                                                                                       >  2:78565/1‑230 (MQ=255)
                                                                                                                    ggAGGTACATTTTCAGTGACCCCGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGTTTTTTCCAGCGGGCAGCACGCCCCCTTTA‑AG‑CGGGAGCAATTTTGCTT‑ATGTTTTACTTCTGCGGAAAAAAACCCTTCACCCCAAA‑CGGAC‑TT‑CCCGGTCCGGcccccccaggtagagacgg                                                                                                                       >  2:374179/1‑232 (MQ=255)
                                                                                                                    ggAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGGCTCACCAGGGt                                                                                                                               <  1:374179/242‑1 (MQ=255)
                                                                                                                                                            cAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGGACGGGCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACCCCCAA‑CGGAC‑CT‑CCCGGTCCGGGTCCCCAGGGTTACGACGCGCTATCTCAAAAACTTTCggggcacaagtgtcacattttat                                                                              >  1:97496/1‑230 (MQ=255)
                                                                                                                                                                    ccGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAACGGCAAGCACGTCCTCTATA‑AG‑CGGGATCCATTTGGTTA‑ATGTTTTTCTACTGCGAAAAATATCCCTTCACGCCAAC‑CGGGCATT‑CCCG‑TCTGGCTCACCCGGGTTACCACGCGCTCGCTCCAAAGCTTTCGGTACAACCGGGTCAAATTTAATGATTTCTg                                                                      >  2:89892/1‑251 (MQ=255)
                                                                                                                                                                     cGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCCCGCACGTCCTCTCTT‑AG‑CGGGATCCATTTTGCTA‑ATTTTTTACTACAGCGGAAAATAACCCTTCACGCCACA‑CGGAC‑CT‑CCCGGTCCGGCCCACCAGGGGTTCGAAGCGCTATCTCAAAAACTTTCGGGACAACAGAGTAAAAATTAATGTTGTCTGc                                                                     >  2:104475/1‑251 (MQ=255)
                                                                                                                                                                     cGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCCGACGCCAGCACGTCCCCTATA‑AG‑CGGGATAAAATTTGCTT‑ATGTTTTTCTACTGCGGCAAAAAACCCCGCACGCCCAA‑CGGCC‑TTTCCCGG‑CCCGCTCCCCACGGGTACGACGCGCCATCTCAAAAACTTTCCGGCCAACCGTGGCCAATTTAATGATGTCTGc                                                                     >  1:442355/1‑251 (MQ=255)
                                                                                                                                                                                aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAAATTGGCTA‑ATGTTTTTCTACTGCGCAAAATAACCCTGCACGCCAAA‑AGGAC‑TT‑TCCCGTCCGGCTCCCCAGGGTTAAGACGCGCTATCTCAAAAACTTTCGGTACAACAGTGTCAAATTTAATTATGTCTGCAATATTGTaa                                                           >  1:118702/1‑250 (MQ=255)
                                                                                                                                                                                        gTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGGGAGGACGTCCTCTTTA‑AG‑AGGAAACAATTTGGGTA‑ATGTTTTACTACTTGGAAAAATATACCTGCCCGGCAAA‑CGGAC‑CT‑CCCCGTCCGGCCCACCAGGGTTTCGACACGCTTTCTCAACAGCTTTCGGTGCAACAGAGTCCAAATTAATGAATTCTTCCATATCGTAAACAAAGGa                                                   >  2:147413/1‑250 (MQ=255)
                                                                                                                                                                                         tGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCCGCGGCAAGCACGGCCTCTTTA‑AA‑CGGAATCAATTTTGCTACAGGTTTT‑CTTCTGCGGAAAAAATACCTCCACACCAAAACGG‑C‑CT‑CCCCGTCCTGCTCCCCAGGGGTTCGCCGCGCCAACTCAAAAGATTTCGGTTCAACAGGGGTAAAATTAATTATGTCTgacaaaatgcaaaaaaaagaa                                                  >  2:339191/1‑230 (MQ=255)
                                                                                                                                                                                                     cGGAATGGCAACAGGCTGTGCAACTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAAGTTAATGATGTCTGGAATATCGTAAACAAAGGAAAGAGGTTTTCCTg                                     <  2:54321/251‑1 (MQ=255)
                                                                                                                                                                                                       gAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCCGCGGGAAGCCCGTCCCCTATA‑AG‑CGGGATCAATTTGGCTA‑ATGGTTTACTACTGCGAAAAATATTCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGGCCCACCAGGGGTACGACGCGCTATTTCAAAAGCTTTTGGTACAACAGTGTGCAATTTAATGATGTCTGCAATAACGTAAACCAAGGAAAAAAGCTTTCCTGt                                    >  1:375966/1‑250 (MQ=255)
                                                                                                                                                                                                                      gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGGCTCACCAGGGTTacgacg                                                                                                                        <  2:311144/151‑1 (MQ=255)
                                                                                                                                                                                                                                  gTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATaa        <  2:377697/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggggCCTGCGTTTGTATTGCGCCAGCGGCATGCACGTCCACTATA‑AG‑CTGAATCAATTTGACTA‑ATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCTGCTCACCAGGGTTACGACGCGCTATCTCAAAATCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTCCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAAGAGCTGGTGCATaccc      <  2:249569/251‑4 (MQ=255)
                                                                                                                                                                                                                                        cctgcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGTCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCt  <  2:362213/248‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATA‑AG‑CGGAATCAATTTGGCTA‑ATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGAC‑TT‑CCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCT  >  NC_000913/2878700‑2879182

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