Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,9470CT27.4% 102.4 / 16.8 63E276K (GAA→AAA) cas1multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (13/32);  new base T (16/1);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  NC_000913/2878703‑2879195
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
ccAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTATTTCTcccgc                                                                                                                                                                                                                                                        >  1:98706/1‑248 (MQ=255)
    gccAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATCCCCCTGCACGTACCTCCAACAACCAGATGGCTAAACTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCAGTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCa                                                                                                                                                                                                                                                    <  2:687032/250‑1 (MQ=255)
           aTTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGt                                                                                                                                                                                                                                             <  1:241744/251‑1 (MQ=255)
            tttGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCCAGCACGTc                                                                                                                                                                                                                                            >  1:289431/1‑251 (MQ=255)
               gTTCCCAGACCATTTCACGAAATTTTGCGGATACATCACCTACATATACCCCTGCACGTCCCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCACCATCCTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGTCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAATCACGTCctc                                                                                                                                                                                                                                         <  2:193215/251‑1 (MQ=255)
                   ccAGATCATTTCACGAAGTATTGCGGATAGATCACCTACATATACCCCTGCACGTACCTGCAACAACCAGAGGGCTAAGCTGCCTCGGAAGCGCGGAGGTACAGTTGCAGTGACCACGACCAACATACTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGTCCGGCGGTTGTATTTCTCCAGCGTCAAGCGCGTACTCtatc                                                                                                                                                                                                                                     <  2:37881/251‑2 (MQ=255)
                     aGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGTCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAg                                                                                                                                                                                                                                   <  1:345096/251‑1 (MQ=255)
                     aGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAg                                                                                                                                                                                                                                   <  1:206796/251‑1 (MQ=255)
                                  aaTTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTg                                                                                                                                                                                                                      <  2:351821/251‑1 (MQ=255)
                                            gATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAAGCTGCCTCTTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAATCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCTGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAAGTTGGCTAATGttt                                                                                                                                                                                                           <  2:426131/251‑1 (MQ=255)
                                              tACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCCGCGGCAAGCACGTCCTCTATAAGCGGAATCACATTGTCTAATGttttt                                                                                                                                                                                                         >  1:568417/1‑250 (MQ=255)
                                                   cACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCACCGGCAAGGACGCCCCCTATACACGGGATCCATTTGGGTAATGTTTTtctact                                                                                                                                                                                                    >  2:444842/1‑251 (MQ=255)
                                                          atatACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGAACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTATTTCTCCATCGGCAAGCACGACCTCTCTAAGCGGGATCCATTTGGCTAATGTTTttctaactcgaaa                                                                                                                                                                                              >  2:22783/1‑242 (MQ=255)
                                                             taCCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCCGCGGCCAACACGTCCTCTTTTAACGGGATCAATTTTGCTAATGTTTTTCTACTGCGAAAAat                                                                                                                                                                                           >  1:606731/1‑250 (MQ=255)
                                                              cccccAGCACGAACCTCCAACAACAAGATAGAGAATCTGCCTCGTAAGCGCGGGAGTACAATTTCAGTGACCACGACCAACAAACTCAATACCGCGACTCCGATGGCCTGCAGCTCCCACTGAAACAGGAAGCGGAATGGCAACAGACTGTGCAACTTCAGGTGGGGACGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAAGCAATTTGGCTAATGTTCTACTACTGCGAAAAatat                                                                                                                                                                                        <  2:613115/250‑1 (MQ=255)
                                                                    gCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCt                                                                                                                                                                                    <  2:633100/249‑1 (MQ=255)
                                                                    gCACGTACCTCCAACAACCAGATGGCTAAGCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCTCGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGTCAATCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGc                                                                                                                                                                                  <  2:283404/251‑1 (MQ=255)
                                                                     cACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGAATCTTCAGGTGGGGCCGGCGGTTGGTATTCTCCCCAGGCAAGAACGTCCTCTTTAAGGGGGAAAAATTTGGCTTATGTTTTATTAATGCGAAAAAAATACCTGcc                                                                                                                                                                                 >  2:533776/1‑250 (MQ=255)
                                                                      aCGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCAc                                                                                                                                                                                <  1:437069/251‑1 (MQ=255)
                                                                             tCCAACAACCAAATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCAACGGCAAGCAAGGCCTCTTTAAACGGAATCAATTTTGGTAATGTTTTACTACTGCGAAAAAAATACCTGCACGCCAAAc                                                                                                                                                                         >  1:83969/1‑251 (MQ=255)
                                                                              ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACg                                                                                                                                                                        <  1:154131/251‑1 (MQ=255)
                                                                              ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACg                                                                                                                                                                        <  1:21197/251‑1 (MQ=255)
                                                                                  caacCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGACt                                                                                                                                                                    <  2:495001/251‑1 (MQ=255)
                                                                                  caacCAGATGGCTAATCTGCCGCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACt                                                                                                                                                                    <  2:107717/251‑1 (MQ=255)
                                                                                                    gCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTc                                                                                                                                                     <  1:68862/247‑1 (MQ=255)
                                                                                                       tCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAACGGCAAGCACGCCCCCTATAAACGGAATCCAATTGGCTAATGTTTTATTACTGCGAAAAAAATTCCATCACGCCAAACAGGCTTcctgggtcgggcccaccggg                                                                                                                                               >  2:379930/1‑233 (MQ=255)
                                                                                                          tAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAACGGCAAGCACGTCCTCTTTAAGCGGAATCAATTTTGCTAATGGTTTACTACTGCGAAAAATATCCCCGCACGCCAAACGGACTTCCCCGTCCGGCTCACCAAGGtt                                                                                                                                            >  1:556078/1‑250 (MQ=255)
                                                                                                           aaGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGAATCTTCAGGAGGGGCCGGCGGTTGTATTTCTCCAGCGGGAAGCCCCTCCTCCATAAACGGAATCAATTTGGCTAATGTTTTACTACTTCGCAAAATATCCCTGCCCGCCCAACGGACTTCCCCGTCCCGCTCACCCGGGGTac                                                                                                                                          >  1:61355/1‑251 (MQ=255)
                                                                                                                  gAGGTACATTTTCAGTGTCCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGAATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCCTCTTCAGGAGGGGCCGGCGGTTGTGTTTCTCCCGCGGCAACGAAGCCCTCTCTCAGCGGGATCTATTTTGCTAATGTTTTACTCCTGCGCAAAATATACCTGCAACCCAAACCGACTTTCCCGTCCGGcacccccggggtacaccgcgc                                                                                                                                   >  2:297974/1‑231 (MQ=255)
                                                                                                                     gTACATGTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTAt                                                                                                                                <  1:613997/251‑1 (MQ=255)
                                                                                                                             tCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTATATTTTTTCCGCGGCAAGCACGGACTCTTTTAGGGGAAACAATTTTGGTAAAATTTTACTACTGCGAAAAATATTCCTGCACGACAAACGGAATTTCCGGTCCGGCCCAACAAGGGTACGAAGCGCtcactccaaa                                                                                                                         >  1:682225/1‑241 (MQ=255)
                                                                                                                                    ccACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCgg                                                                                                                <  1:49858/251‑1 (MQ=255)
                                                                                                                                        gACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCAAGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGGTAAGCAGCGCTATCTCCAAAGCTTTCGGTac                                                                                                             >  1:454229/1‑250 (MQ=255)
                                                                                                                                         aCCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTGTTTTTCCAGCGGGAAGCACGGCCTCTATAAGCGGAATCAATTTTGCTAATGTTTTAATACTGGGGAAAATATCCCTTCACCCCAAAAGGACTTCCCCGGCCCGCTCACCAGGGGTACGACGCGCCATCTCAAAAGCTTTCGGTaca                                                                                                            >  2:613958/1‑250 (MQ=255)
                                                                                                                                         aCCAACATACTCATTTCAGCTACTCCGATGGCCCGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGATTTTTTTCAGAGGCAAGCACGTCCTCTATAAAAGGAATCAATTTAGCTAATGTTTTACTACTTCGAAAAATATCCCTGCCAGCCCAACTGACTTTCCAGTCCTGGCCACCACGGTTACAGCCCCCTATCTCAAAAGCTTTCGGTaacc                                                                                                           >  1:584978/1‑248 (MQ=255)
                                                                                                                                               aTACTCATTTCAGCTACGTCGATGGCCTGCATCTCCCAGGGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGTGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCGAAACGGACTTCCCGGTCCGGCGCACCAGGGTTACGACGCGCTATCTCAAAAGCGTTCGGTAAAACAgtgt                                                                                                     <  2:79273/251‑1 (MQ=255)
                                                                                                                                                       ttCAGCGACTCCGCTGCCCAGAATCTCCCAGTGAAGAAGGAGTCGGAATGGTACCAGGCTGTCCATCTTCAGGTGGTCTCGGCATTTGTATTTCTCCAGCTTCAAGCTCGTCTTCTATAAGCGGAACCAATTTGGCTAATGTTTTAATATTGCGAAAAATATCACTGCTCGCCAAACGTTCTTCCCGGCCCGGCTAACCAGGCTTACGACGCGCTATCTCGAAAGCTTTCGGTACAACAGTGTCAAAttt                                                                                              <  1:484881/250‑1 (MQ=255)
                                                                                                                                                         cAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTGTTTTTCCAGCGGCAAGCACGTCCTCTATAAACGGGATCAATTTTGCTAATGTTTTTATACTGCGAAAAATATTCCTGCCCCCCAAACGGACTTTCCGGGCCGGCCCACCCGGGTTACGACCCGCTATTTCCAAACCTTTCGGTACAACAGAGGCAAATTTAat                                                                                           >  1:64108/1‑251 (MQ=255)
                                                                                                                                                                        ccTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATc                                                                            <  1:143998/251‑1 (MQ=255)
                                                                                                                                                                        ccTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGGGGGGCCGGCGGTTTGATTTTTCCAGCGGCCAGCCCGGCCTCTTTAAACGGGATCAATTTTGGTAATTGTTTTCTACTGCGGAAAAAAACCCTGCACGCCCAAAGGACTTCCCGGGCCGGCCCCCCAGGGGTAACACGCGCCATCTCCAAAACTTTTGGGTCAACAGGGTCCAATTTAAAGATGTCTGCAATATc                                                                            >  2:209372/1‑251 (MQ=255)
                                                                                                                                                                             aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAATTTGGGTAATGTTTTTCTCCTGCGCAAAATATCCCTGCACGCCCCACGGACTTCCCGGTCCGGCTCACCAGGGTTAACACCCGCTATCTCAAAAGCTTTTGGTGCAACAGTGTCAAATTTAATTATGTCTGCCATATCGTaa                                                                        >  1:350482/1‑250 (MQ=255)
                                                                                                                                                                             aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGGGGGGCCGGCGGTTGTATTTCTCCAAAGGCAAACCCGTGCTCTTTAAACGGAATCAATTTGGCTAATGTTTTACTACTGTGGAAAATATCCCTGCACGCCAAACGGGATTCCCGGTCCGGCTAACCAGGGGTAACACGCGGCATCTCAAAAAGCTTCTGGACAACCGTGGCAAATTTTATAATTTCTGCGATATCGTaa                                                                        >  1:128574/1‑250 (MQ=255)
                                                                                                                                                                                 cccAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCGTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATAGCCCTGCACGCCAAACTGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACa                                                                     <  2:98706/249‑1 (MQ=255)
                                                                                                                                                                                 cccAGTGAAACAGGAAGCGGAATGACAACAGGCTGTGCATCTGCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGTCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTCTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACaaa                                                                   <  1:377166/251‑1 (MQ=255)
                                                                                                                                                                                   cAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATGTGGCTAATGTTGTACTACTGCGAAAAATATCCCTGCACGGCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCaaaa                                                                                                                         <  1:683502/196‑1 (MQ=255)
                                                                                                                                                                                    aGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCACGGCAAGCACGTCCTCTATAAACGGGATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGGCCGGCTCACCAGGGTTACGACGCGCTATCTCCAAAACTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACaaa                                                                   >  1:129081/1‑248 (MQ=255)
                                                                                                                                                                                       gAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGAGGGAAGCCCCTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGGCCGGCTCACCAGGGGTACGGCGCGCCATCTCAAAAGCTTTCGGTCCAACCGTGTCAAATTTAATGATGTCTGCCATATCGTAAACAAAGGaaa                                                              >  1:465823/1‑250 (MQ=255)
                                                                                                                                                                                        aaaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGGAAGCACCGCCTCTATAAGCGGAATCAATTTTGCTAATGTTTTACTCCTGCGAAAATAATCCCTGCACGCCCAACGGACTTTCCCGTCCCGCTCACCCGGGTTACGACGCGCTATCTCAAAAGGTTTCGGTACAACAGTGTCCAATTTAATGATGTCTGCCatagacaaaacaaaagaa                                                               >  2:185305/1‑233 (MQ=255)
                                                                                                                                                                                            aGGAAGCGGAATGGCAACAGGCTGTTCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATAGCGTAAACAAAGGAAAGAGGCt                                                        <  1:599002/251‑1 (MQ=255)
                                                                                                                                                                                            aGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCCAGCACGTCCTCTATAAACGGGATCAATTTGGCTAATGTTTTTCTTCTGCGAAAAAAAACCCTGGACGCCAAACGGGCTTCCCGGGCCGGCCCCCCAGGGTTTCGACGCGCTATCTCCCAAGCTTTCGGTACCACAGGGTCCAATTTAATTGTGTCTTCAATATTCTAAACAAAgggaaaagggt                                                        >  2:310855/1‑242 (MQ=255)
                                                                                                                                                                                            aGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAACACGTCCTCTCTAAGCGGGATCAATTTGGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCCAACGGACTTCCCGGTCCGGCTCACCCGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAAAGATGTCTGCCATATCCTAAAAAAAGGAAAGAgggt                                                        >  1:614076/1‑249 (MQ=255)
                                                                                                                                                                                                aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTcc                                                    <  1:325829/251‑1 (MQ=255)
                                                                                                                                                                                                 gCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTCTTTTTCCAGCAGGAAGCCCGGTCTCTATAAGGGGGAACAAATTGGCTAATGTTTTACTACTGCGGAAAATAATCCCTCCCCCCAAAACGACTTCCCCGTCCGGGTCCCCCGGGGTACGCAGCCCCTTCTCCACAACTTTTGGTACACAAAAGTCAAATTTAATGATGTATGCACTAATGTAAACAAAGGaaaaggggttttct                                                   >  2:636526/1‑240 (MQ=255)
                                                                                                                                                                                                    gAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCaa                                                                                                                           <  1:491025/177‑1 (MQ=255)
                                                                                                                                                                                                        ggCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAACGGCAAGCACGTCCCCTATAAGCGGGATCAATTTGGCTAAAGTTTTTCTACTTCGAAAAATATCCCTGCCCGCCCAACGGACTTCCCGGTCCCGCTCACCAGGGGTACGACGCCCCATCTCCAAAGCTTTCGGTACAACAATGGCAAATTTAATGATGTCTGCCATATCGTAAACAAAGGAAAAAAGCTTTCCttttttcc                                            >  1:273424/1‑244 (MQ=255)
                                                                                                                                                                                                         gCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGAATTTTTCCACCGACAAACACATCCCCTTTAAGCGGGATCAATTTTGCCAATGTTTTAATACTTCGAAAAAAATCCCTTCACCCCAAAACGACTTCCCCGTCCGGCTCACCCAGGTTACGACCCCCTTTCTCAAAAACTTTCGGTTCAACAGTGTCAAATTTTATGATTTCTTCAATATCTTTAACAAAGGAAAAATGCTTTCCTTTATGccc                                           >  1:67742/1‑249 (MQ=255)
                                                                                                                                                                                                                   gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCaa                                   <  1:636526/249‑1 (MQ=255)
                                                                                                                                                                                                                   gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAAt                                  <  1:379930/250‑1 (MQ=255)
                                                                                                                                                                                                                           tCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTg                             <  1:310855/247‑1 (MQ=255)
                                                                                                                                                                                                                                tGGGGCCGGCGGTTGTTTTTTTCCCGCGGCAAGAACATCCTCTTTAAACGGAATCAATTTTGTTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGGCTTCCCGGTCCGGCGCCCCAGGGGTACGACGCGCTATCTCACAAGATTTTGGGACAAAAATGGCAAAATTAATGATTTCTGCAATATCGTAAAAAAAGGAAAGAGGCATTCCTGTATGCACAAACCCAAAAGCGGGTGCAAAAc                    >  2:285346/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccGGCGGTTGTATTGCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCTCGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACgcgc                                                                                                                                   <  1:337645/136‑1 (MQ=255)
                                                                                                                                                                                                                                        gaggTTGTATTTATCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGATTTACTACTGCGAAACATATCCCTGCACGCCAAACGGACTTCCCGGTCAGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGACAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCGTTCCTATATGCACAAACCAAATAGCTGGAGCATAACCAGCTGc             >  2:578263/3‑250 (MQ=255)
                                                                                                                                                                                                                                           gTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGt         <  2:350482/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    tctcCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgccg  <  1:41706/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
CCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  NC_000913/2878703‑2879195

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