Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,9470CT26.9% 98.2 / 10.1 52E276K (GAA→AAA) cas1multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (13/25);  new base T (14/0);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.39e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCC  >  NC_000913/2878699‑2879194
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                              
tcCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTtc                                                                                                                                                                                                                                                                <  1:115655/249‑1 (MQ=255)
         ccAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTT‑cccgctgcc                                                                                                                                                                                                                                                     >  2:4954/1‑243 (MQ=255)
         ccAGCTATTTGTTCCCAGATCAGTTCACGAATTGTGGCGGATACATCCCCTCCAGATACCCCTGCACGTACCTCCAACAACCAGATGTCTAATCTGCCTCGTAAGCGCGGAGTTACATTTTCAGTGACCACGTCCAACATACTCATTTCAGCTACTACGATGTCCTGCATCTCGCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGTGGTTGTATTTCT‑CCAGCGG‑Caa                                                                                                                                                                                                                                                   <  2:455501/251‑1 (MQ=255)
                tttGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTT‑CCAGCGG‑CAAACACGTc                                                                                                                                                                                                                                            >  1:104317/1‑251 (MQ=255)
                     tCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCctct                                                                                                                                                                                                                                        <  1:204074/250‑1 (MQ=255)
                        cAGATCCTTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATCCTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTT‑CCAACGG‑CAAGcccctactctctc                                                                                                                                                                                                                                     >  2:157802/1‑238 (MQ=255)
                        cAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATaa                                                                                                                                                                                                                                    <  1:252424/251‑1 (MQ=255)
                             cATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTT‑CCAGCGG‑GCAGAACGTCCTCTATAAGCggg                                                                                                                                                                                                                               >  1:51074/1‑250 (MQ=255)
                                      aaTTTTTGCGGATACATCACATACCTATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCCGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACCGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTATTTTT‑ACAGCGG‑CACCAAAGTCCACTATAAGCGG‑aaaaacattt                                                                                                                                                                                                                     >  2:415514/1‑243 (MQ=255)
                                            tGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGTACGTCCTCTATAAGCGG‑AATCAATTTGGCTaa                                                                                                                                                                                                                <  1:602682/250‑1 (MQ=255)
                                               ggATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCCCTTTAAGCGG‑AATCCAttttggtattggt                                                                                                                                                                                                            >  2:210056/1‑241 (MQ=255)
                                                     aTCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTT‑CCACCGG‑CAAGCACGACCTCTCTAAGCGG‑AATCAATTTTGCCCATGtttttctt                                                                                                                                                                                                      >  2:346832/1‑247 (MQ=255)
                                                       cgccTACATATGCCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGCCCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGTCCGGCGGTTGTATTTCC‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTtactac                                                                                                                                                                                                     <  1:553719/248‑1 (MQ=255)
                                                       cACCTACATAAACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGT‑AATCAATTTGGCTAATGTTTTactact                                                                                                                                                                                                    <  2:136666/251‑1 (MQ=255)
                                                            aCATATACCCCTGCACGTACCTCCAACATCCAGATTGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCATTGACCACGACCACCATACGCATTTCAGCTACTCCGATGGCCTGCTTCTCCCAGTGAACCAGGAATCGGAATGGCAACAGTCTTTGCTTCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCTG‑CAAGCACGTCCTCTTGAAGCGG‑AACCTATTTGGCTATTTGTTTACTACTGct                                                                                                                                                                                                 <  2:539072/249‑2 (MQ=255)
                                                                 taCCCATGCACGTACCTCCAAAAACCAGATGGCTAATCTGCCTCGTACGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTT‑CCACCGGGCAA‑CACGTCCTCTATTAGCGG‑GAGCAATTTGGCTTATGTTTTATCACTGCGGAAAat                                                                                                                                                                                           >  2:203378/1‑250 (MQ=255)
                                                                   ccccTGCACGTTCCTCCATCAACCAGATGGCTAATCTGACGCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTGGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCTAAAAATATc                                                                                                                                                                                        <  2:497463/251‑1 (MQ=255)
                                                                        gCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCt                                                                                                                                                                                     <  2:149525/249‑1 (MQ=255)
                                                                            gTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTCT‑CCACCGG‑CCAGCCCCCTCTCTATAAACGG‑AATCAATTTTGCTAATGTTTTTCTACTGCGaaaaaaaaccccgccccc                                                                                                                                                                               >  2:556439/1‑238 (MQ=255)
                                                                                  ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCaaa                                                                                                                                                                           <  1:481540/249‑1 (MQ=255)
                                                                                      caacCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCCA‑CGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAAATTGGGTCATTTTTTATTTCTGCGGACAATATTCCTGCAACCCAAA‑CGGACt                                                                                                                                                                    >  1:200657/1‑251 (MQ=255)
                                                                                                              tAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTT‑TCCGCGG‑GAAACACGTCCTCTATAAGCGG‑AATCAATTTTGCTAAAGTTTTTAAACTGCGAAAAATATCCCTGCACGCCCAA‑CGGACTTTCCCGTCCCGCTCCCCAggggt                                                                                                                                             >  1:14848/1‑248 (MQ=255)
                                                                                                                   gcgGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTacg‑acg                                                                                                                                      <  1:395111/251‑1 (MQ=255)
                                                                                                                        ggTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTCT‑CCAGCGG‑CAAGCACCTCCTCTATAAGCGG‑AATCAATTTGGCTAAAGTTTTTCTTCTGCGAAAAAAATCCCTGCACCCCAAA‑ACGACTTCCCGGGCCGGCCCACCAgggataagccccgct                                                                                                                                   >  1:590071/1‑238 (MQ=255)
                                                                                                                         gTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTg                                                                                                                                                                                    <  1:187332/201‑1 (MQ=255)
                                                                                                                              ttttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCACCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAAATTGGGTAATGGTTTACTACTACGAAAAAAATCCCTGCACGCCAAA‑CGGACTTCCCCGTCCCGCTCACCAGGGTTacg‑aaccgccatctcc                                                                                                                            >  1:466709/1‑238 (MQ=255)
                                                                                                                                  cAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTT‑TCAGCGG‑CAAGCACGGCCTCTATAAACGG‑GAACAATTTGGCTAAAGTTTTACTACTGCGGAAAATAACCCCGCACGCCACA‑CGGACTTCCCGGTCCGGCTCACCAGGGGTACG‑ACCCGCTAtccccaaaac                                                                                                                       >  2:610935/1‑243 (MQ=255)
                                                                                                                                  cAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGc                                                                                                                      <  1:451813/251‑1 (MQ=255)
                                                                                                                                      gACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCttt                                                                                                                   <  1:210056/250‑1 (MQ=255)
                                                                                                                                      gACCACGACCAACATACTCATGTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCTGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTc                                                                                                                  <  2:522363/251‑1 (MQ=255)
                                                                                                                                             aCCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTT‑CCCGCGG‑CAAGCACCGCCTCTATAAGCGG‑GATCAATTTTGCTAATGGTTTACTACTGCGGAAAAAATCCCTGCCCCCCAAA‑CGGACTTCCCGGTCCGGCCCACCAGGGGTACG‑CCCCGCCCTCTC‑CAAAGCTTTCggggca                                                                                                            >  2:284812/1‑246 (MQ=255)
                                                                                                                                              ccAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTacaa‑c                                                                                                         <  1:37040/251‑1 (MQ=255)
                                                                                                                                                   attcTCATTTCAGCTACTCCGAAGGCCTGCATCTCCCATTGAAACAGGACGCGTAATGGCAACAGGCTGTGCAGCTTCAGGTGGGGCCGGCTTTTGTATGTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATTTTTTACTACTGCGAAAAATTTCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACTAGGGTTACT‑ACGCGCTATCTC‑AAAAGCTTTCGTTACAA‑CAgtgt                                                                                                    <  2:429730/248‑1 (MQ=255)
                                                                                                                                                   aTACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAgtgt                                                                                                    <  1:314737/251‑1 (MQ=255)
                                                                                                                                                        cATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTATTTCT‑CCCGCGG‑CAAGCACATCCTCTCTAAGCGG‑AATCAATTTTGCTTATGTTTTTCTACTGCGACAAATATCCCTGCCCCCCACA‑AGGGCTTCCCCGTCCGGCTCCCCAAGGGTACG‑ACGCCCTATCCC‑AAAAGCCTTctgtgccacagagccaaa                                                                                                 >  2:300277/1‑233 (MQ=255)
                                                                                                                                                         aTTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAAtt                                                                                              <  1:81578/251‑1 (MQ=255)
                                                                                                                                                             cAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CCAGCACGTCCTCTATAAGCGG‑AATCAATTTTGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑AGGACTTCCCGGTCCGGCCCACCCGGGGTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAat                                                                                          >  1:196839/1‑251 (MQ=255)
                                                                                                                                                                    tCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTT‑CCAGCGG‑CCAGCACGTCCCCTATAAGCGGGAAT‑AATTTGGGTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCCGTCCGGCTCACCCGGGGTACG‑ACGCGCTATCTC‑AAAAACTTTCGGGACCA‑CCGGGTCAAATTTTATGATGTCt                                                                                   >  2:467630/1‑251 (MQ=255)
                                                                                                                                                                     ccGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTc                                                                                    <  1:400741/249‑1 (MQ=255)
                                                                                                                                                                          ggCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAat                                                                              <  2:104317/250‑1 (MQ=255)
                                                                                                                                                                            ccTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTTTAAGCGG‑AATCAATTTTGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CCGACTTCCCGGTCCGGCCCACCAGGGTTACG‑ACACGCTATCTC‑CAAAGCTTTCGGTACAA‑CAGTGTCCAATTTTATGATGTCTGCAATATc                                                                           >  1:390307/1‑251 (MQ=255)
                                                                                                                                                                             cTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTATTTTT‑CCAGCGG‑GCAGCACGGCCTCTATAAACGG‑AAACAATTTGGCTAATGTTTTTCTACCGGGAAAAATATCCCTGCACGCCAAAAGGGCCTTCCCGG‑CCCGCTCCCCAGGGGTACG‑GAGCGCTATCTCCAAAACCTTT‑GGGACAA‑CCGGGGCAAAATTTATGATGTTTGCAAAATCg                                                                          >  1:277722/1‑251 (MQ=255)
                                                                                                                                                                                           gAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTATTTCT‑CCAGCGG‑AAAGAACGTCCCCTATAAACGG‑AATCAAATTGGCTAATGTTTTACTACTTCGAAAAAAATCCCTTCACGCCAAA‑CGGACTTCCCGGTTCGGGTTACCAGGGGTACGCAC‑CGCTATCTC‑AAAAAGTTTCGGTACAAACGGGGTCAA‑TCTTATGGTGTCTGCAATATAGGAAACAAAGGaa                                                              >  2:407250/1‑249 (MQ=255)
                                                                                                                                                                                                   aaGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTATTTTT‑CCAGCGG‑CAAGCACGTCCTCTTTAAGAGG‑AATCAATTTGGCTAATGTTTTTTTACTGCGGAAAAAATCCCTGCACGCCCAA‑CGGACTTCCCGGTCCGGCTCACCCGGGGTACG‑ACGCGCTATCTC‑AAAAGCTTTCTGTACAA‑CAGTGTCAAATTTAATTATGTCTGCAATATCGTAAACAAAAGAAAAAGGcgttt                                                    >  2:122285/1‑247 (MQ=255)
                                                                                                                                                                                                    aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTCT‑TCAGCGG‑CCAACACGTCCTCTTTAAGAGG‑GATCAAATTGGCTAATGTTTTAATACTTCGAAAAAAATACCTGCCCGCCAAA‑CGGACTTTCCGGGCCGGCTCACCAGGGGTACG‑AACCGCTATCTC‑AAAAACCTTCGGGACAA‑CAATGGCAAAATTAAAGATGTCTGCAATATCGTAAACAAAAGAAAGAGGCTTTcc                                                   >  1:554596/1‑251 (MQ=255)
                                                                                                                                                                                                            ggCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGc                                            <  1:317692/250‑1 (MQ=255)
                                                                                                                                                                                                                  aGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCt                                                       <  2:639538/233‑1 (MQ=255)
                                                                                                                                                                                                                  aGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAACGG‑CCAGCAAGGCCGCCTTAAGCGG‑AATCAATTTTGCTAATGTTTTAATAATGCGAAAAATATTCCTCCACGCCAAA‑CCGACTTCCCGGTCCGGGTCAACCGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGAGGCAAAATTAAAGATGGCTGCAAAATCGTAAACAAAAGAAAGAgggtagatcggaagagcaaacc                                     >  1:639538/1‑231 (MQ=255)
                                                                                                                                                                                                                                       ggCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAAcc                  >  1:146544/1‑249 (MQ=255)
                                                                                                                                                                                                                                          cggcggTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTg             <  1:203378/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cggTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCa           >  2:552211/1‑250 (MQ=255)
                                                                                                                                                                                                                                               gTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTAATGCGAAAAATATCCCTGCACGCCAAA‑CGGCCTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGAGCATAACCAGCTGCAAgg        >  2:498026/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATc     <  1:4954/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      tcct‑cCAGCGG‑CTAGCACGTCCTCTATAAACGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgct  <  2:394128/248‑2 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                              
TCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCT‑CCAGCGG‑CAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAA‑CGGACTTCCCGGTCCGGCTCACCAGGGTTACG‑ACGCGCTATCTC‑AAAAGCTTTCGGTACAA‑CAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCC  >  NC_000913/2878699‑2879194

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