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.

TCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCC  >  NC_000913/2878699‑2879194
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                              
tcCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTtc                                                                                                                                                                                                                                                                <  1:115655/249‑1 (MQ=255)
         ccAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTcccgctgcc                                                                                                                                                                                                                                                     >  2:4954/1‑243 (MQ=255)
         ccAGCTATTTGTTCCCAGATCAGTTCACGAATTGTGGCGGATACATCCCCTCCAGATACCCCTGCACGTACCTCCAACAACCAGATGTCTAATCTGCCTCGTAAGCGCGGAGTTACATTTTCAGTGACCACGTCCAACATACTCATTTCAGCTACTACGATGTCCTGCATCTCGCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGTGGTTGTATTTCTCCAGCGGCaa                                                                                                                                                                                                                                                   <  2:455501/251‑1 (MQ=255)
                tttGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCAGCGGCAAACACGTc                                                                                                                                                                                                                                            >  1:104317/1‑251 (MQ=255)
                     tCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCctct                                                                                                                                                                                                                                        <  1:204074/250‑1 (MQ=255)
                        cAGATCCTTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATCCTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTTCCAACGGCAAGcccctactctctc                                                                                                                                                                                                                                     >  2:157802/1‑238 (MQ=255)
                        cAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATaa                                                                                                                                                                                                                                    <  1:252424/251‑1 (MQ=255)
                             cATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCAGCGGGCAGAACGTCCTCTATAAGCggg                                                                                                                                                                                                                               >  1:51074/1‑250 (MQ=255)
                                      aaTTTTTGCGGATACATCACATACCTATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCCGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACCGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTATTTTTACAGCGGCACCAAAGTCCACTATAAGCGGaaaaacattt                                                                                                                                                                                                                     >  2:415514/1‑243 (MQ=255)
                                            tGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGTACGTCCTCTATAAGCGGAATCAATTTGGCTaa                                                                                                                                                                                                                <  1:602682/250‑1 (MQ=255)
                                               ggATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCCCTTTAAGCGGAATCCAttttggtattggt                                                                                                                                                                                                            >  2:210056/1‑241 (MQ=255)
                                                     aTCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTTCCACCGGCAAGCACGACCTCTCTAAGCGGAATCAATTTTGCCCATGtttttctt                                                                                                                                                                                                      >  2:346832/1‑247 (MQ=255)
                                                       cgccTACATATGCCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGCCCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGTCCGGCGGTTGTATTTCCCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtactac                                                                                                                                                                                                     <  1:553719/248‑1 (MQ=255)
                                                       cACCTACATAAACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGTAATCAATTTGGCTAATGTTTTactact                                                                                                                                                                                                    <  2:136666/251‑1 (MQ=255)
                                                            aCATATACCCCTGCACGTACCTCCAACATCCAGATTGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCATTGACCACGACCACCATACGCATTTCAGCTACTCCGATGGCCTGCTTCTCCCAGTGAACCAGGAATCGGAATGGCAACAGTCTTTGCTTCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCTGCAAGCACGTCCTCTTGAAGCGGAACCTATTTGGCTATTTGTTTACTACTGct                                                                                                                                                                                                 <  2:539072/249‑2 (MQ=255)
                                                                 taCCCATGCACGTACCTCCAAAAACCAGATGGCTAATCTGCCTCGTACGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCACCGGGCAACACGTCCTCTATTAGCGGGAGCAATTTGGCTTATGTTTTATCACTGCGGAAAat                                                                                                                                                                                           >  2:203378/1‑250 (MQ=255)
                                                                   ccccTGCACGTTCCTCCATCAACCAGATGGCTAATCTGACGCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTGGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCTAAAAATATc                                                                                                                                                                                        <  2:497463/251‑1 (MQ=255)
                                                                        gCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCt                                                                                                                                                                                     <  2:149525/249‑1 (MQ=255)
                                                                            gTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTCTCCACCGGCCAGCCCCCTCTCTATAAACGGAATCAATTTTGCTAATGTTTTTCTACTGCGaaaaaaaaccccgccccc                                                                                                                                                                               >  2:556439/1‑238 (MQ=255)
                                                                                  ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCaaa                                                                                                                                                                           <  1:481540/249‑1 (MQ=255)
                                                                                      caacCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCCACGGCAAGCACGTCCTCTATAAGCGGAATCAAATTGGGTCATTTTTTATTTCTGCGGACAATATTCCTGCAACCCAAACGGACt                                                                                                                                                                    >  1:200657/1‑251 (MQ=255)
                                                                                                              tAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTTCCGCGGGAAACACGTCCTCTATAAGCGGAATCAATTTTGCTAAAGTTTTTAAACTGCGAAAAATATCCCTGCACGCCCAACGGACTTTCCCGTCCCGCTCCCCAggggt                                                                                                                                             >  1:14848/1‑248 (MQ=255)
                                                                                                                   gcgGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTacgacg                                                                                                                                      <  1:395111/251‑1 (MQ=255)
                                                                                                                        ggTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTCTCCAGCGGCAAGCACCTCCTCTATAAGCGGAATCAATTTGGCTAAAGTTTTTCTTCTGCGAAAAAAATCCCTGCACCCCAAAACGACTTCCCGGGCCGGCCCACCAgggataagccccgct                                                                                                                                   >  1:590071/1‑238 (MQ=255)
                                                                                                                         gTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTg                                                                                                                                                                                    <  1:187332/201‑1 (MQ=255)
                                                                                                                              ttttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCACCGGCAAGCACGTCCTCTATAAGCGGAATCAAATTGGGTAATGGTTTACTACTACGAAAAAAATCCCTGCACGCCAAACGGACTTCCCCGTCCCGCTCACCAGGGTTacgaaccgccatctcc                                                                                                                            >  1:466709/1‑238 (MQ=255)
                                                                                                                                  cAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCAGCGGCAAGCACGGCCTCTATAAACGGGAACAATTTGGCTAAAGTTTTACTACTGCGGAAAATAACCCCGCACGCCACACGGACTTCCCGGTCCGGCTCACCAGGGGTACGACCCGCTAtccccaaaac                                                                                                                       >  2:610935/1‑243 (MQ=255)
                                                                                                                                  cAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGc                                                                                                                      <  1:451813/251‑1 (MQ=255)
                                                                                                                                      gACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCttt                                                                                                                   <  1:210056/250‑1 (MQ=255)
                                                                                                                                      gACCACGACCAACATACTCATGTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCTGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTc                                                                                                                  <  2:522363/251‑1 (MQ=255)
                                                                                                                                             aCCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCCGCGGCAAGCACCGCCTCTATAAGCGGGATCAATTTTGCTAATGGTTTACTACTGCGGAAAAAATCCCTGCCCCCCAAACGGACTTCCCGGTCCGGCCCACCAGGGGTACGCCCCGCCCTCTCCAAAGCTTTCggggca                                                                                                            >  2:284812/1‑246 (MQ=255)
                                                                                                                                              ccAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTacaac                                                                                                         <  1:37040/251‑1 (MQ=255)
                                                                                                                                                   attcTCATTTCAGCTACTCCGAAGGCCTGCATCTCCCATTGAAACAGGACGCGTAATGGCAACAGGCTGTGCAGCTTCAGGTGGGGCCGGCTTTTGTATGTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATTTTTTACTACTGCGAAAAATTTCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACTAGGGTTACTACGCGCTATCTCAAAAGCTTTCGTTACAACAgtgt                                                                                                    <  2:429730/248‑1 (MQ=255)
                                                                                                                                                   aTACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAgtgt                                                                                                    <  1:314737/251‑1 (MQ=255)
                                                                                                                                                        cATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTATTTCTCCCGCGGCAAGCACATCCTCTCTAAGCGGAATCAATTTTGCTTATGTTTTTCTACTGCGACAAATATCCCTGCCCCCCACAAGGGCTTCCCCGTCCGGCTCCCCAAGGGTACGACGCCCTATCCCAAAAGCCTTctgtgccacagagccaaa                                                                                                 >  2:300277/1‑233 (MQ=255)
                                                                                                                                                         aTTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAAtt                                                                                              <  1:81578/251‑1 (MQ=255)
                                                                                                                                                             cAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCCAGCACGTCCTCTATAAGCGGAATCAATTTTGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAAAGGACTTCCCGGTCCGGCCCACCCGGGGTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAat                                                                                          >  1:196839/1‑251 (MQ=255)
                                                                                                                                                                    tCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCCAGCACGTCCCCTATAAGCGGGAATAATTTGGGTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCCGTCCGGCTCACCCGGGGTACGACGCGCTATCTCAAAAACTTTCGGGACCACCGGGTCAAATTTTATGATGTCt                                                                                   >  2:467630/1‑251 (MQ=255)
                                                                                                                                                                     ccGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTc                                                                                    <  1:400741/249‑1 (MQ=255)
                                                                                                                                                                          ggCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAat                                                                              <  2:104317/250‑1 (MQ=255)
                                                                                                                                                                            ccTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTTTAAGCGGAATCAATTTTGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACCGACTTCCCGGTCCGGCCCACCAGGGTTACGACACGCTATCTCCAAAGCTTTCGGTACAACAGTGTCCAATTTTATGATGTCTGCAATATc                                                                           >  1:390307/1‑251 (MQ=255)
                                                                                                                                                                             cTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTATTTTTCCAGCGGGCAGCACGGCCTCTATAAACGGAAACAATTTGGCTAATGTTTTTCTACCGGGAAAAATATCCCTGCACGCCAAAAGGGCCTTCCCGGCCCGCTCCCCAGGGGTACGGAGCGCTATCTCCAAAACCTTTGGGACAACCGGGGCAAAATTTATGATGTTTGCAAAATCg                                                                          >  1:277722/1‑251 (MQ=255)
                                                                                                                                                                                           gAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTATTTCTCCAGCGGAAAGAACGTCCCCTATAAACGGAATCAAATTGGCTAATGTTTTACTACTTCGAAAAAAATCCCTTCACGCCAAACGGACTTCCCGGTTCGGGTTACCAGGGGTACGCACCGCTATCTCAAAAAGTTTCGGTACAAACGGGGTCAATCTTATGGTGTCTGCAATATAGGAAACAAAGGaa                                                              >  2:407250/1‑249 (MQ=255)
                                                                                                                                                                                                   aaGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTATTTTTCCAGCGGCAAGCACGTCCTCTTTAAGAGGAATCAATTTGGCTAATGTTTTTTTACTGCGGAAAAAATCCCTGCACGCCCAACGGACTTCCCGGTCCGGCTCACCCGGGGTACGACGCGCTATCTCAAAAGCTTTCTGTACAACAGTGTCAAATTTAATTATGTCTGCAATATCGTAAACAAAAGAAAAAGGcgttt                                                    >  2:122285/1‑247 (MQ=255)
                                                                                                                                                                                                    aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTCTTCAGCGGCCAACACGTCCTCTTTAAGAGGGATCAAATTGGCTAATGTTTTAATACTTCGAAAAAAATACCTGCCCGCCAAACGGACTTTCCGGGCCGGCTCACCAGGGGTACGAACCGCTATCTCAAAAACCTTCGGGACAACAATGGCAAAATTAAAGATGTCTGCAATATCGTAAACAAAAGAAAGAGGCTTTcc                                                   >  1:554596/1‑251 (MQ=255)
                                                                                                                                                                                                            ggCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGc                                            <  1:317692/250‑1 (MQ=255)
                                                                                                                                                                                                                  aGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCt                                                       <  2:639538/233‑1 (MQ=255)
                                                                                                                                                                                                                  aGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAACGGCCAGCAAGGCCGCCTTAAGCGGAATCAATTTTGCTAATGTTTTAATAATGCGAAAAATATTCCTCCACGCCAAACCGACTTCCCGGTCCGGGTCAACCGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGAGGCAAAATTAAAGATGGCTGCAAAATCGTAAACAAAAGAAAGAgggtagatcggaagagcaaacc                                     >  1:639538/1‑231 (MQ=255)
                                                                                                                                                                                                                                       ggCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAAcc                  >  1:146544/1‑249 (MQ=255)
                                                                                                                                                                                                                                          cggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTg             <  1:203378/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCa           >  2:552211/1‑250 (MQ=255)
                                                                                                                                                                                                                                               gTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTAATGCGAAAAATATCCCTGCACGCCAAACGGCCTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGAGCATAACCAGCTGCAAgg        >  2:498026/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATc     <  1:4954/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      tcctcCAGCGGCTAGCACGTCCTCTATAAACGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgct  <  2:394128/248‑2 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                              
TCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCC  >  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: