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.

CCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  NC_000913/2878703‑2879195
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
ccAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTT‑TATTTCTcccgc                                                                                                                                                                                                                                                        >  1:98706/1‑248 (MQ=255)
    gccAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATCCCCCTGCACGTACCTCCAACAACCAGATGGCTAAACTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCAGTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCa                                                                                                                                                                                                                                                    <  2:687032/250‑1 (MQ=255)
           aTTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGt                                                                                                                                                                                                                                             <  1:241744/251‑1 (MQ=255)
            tttGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCCAGCACGTc                                                                                                                                                                                                                                            >  1:289431/1‑251 (MQ=255)
               gTTCCCAGACCATTTCACGAAATTTTGCGGATACATCACCTACATATACCCCTGCACGTCCCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCACCATCCTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGTCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCAGCGGCAATCACGTCctc                                                                                                                                                                                                                                         <  2:193215/251‑1 (MQ=255)
                   ccAGATCATTTCACGAAGTATTGCGGATAGATCACCTACATATACCCCTGCACGTACCTGCAACAACCAGAGGGCTAAGCTGCCTCGGAAGCGCGGAGGTACAGTTGCAGTGACCACGACCAACATACTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGTCCGGCGGTTG‑TATTTCTCCAGCGTCAAGCGCGTACTCtatc                                                                                                                                                                                                                                     <  2:37881/251‑2 (MQ=255)
                     aGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGTCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAg                                                                                                                                                                                                                                   <  1:345096/251‑1 (MQ=255)
                     aGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAg                                                                                                                                                                                                                                   <  1:206796/251‑1 (MQ=255)
                                  aaTTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTg                                                                                                                                                                                                                      <  2:351821/251‑1 (MQ=255)
                                            gATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAAGCTGCCTCTTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAATCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCTGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAAGTTGGCTAAT‑Gttt                                                                                                                                                                                                           <  2:426131/251‑1 (MQ=255)
                                              tACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCCGCGGCAAGCACGTCCTCTATAAGCGGAATCACATTGTCTAAT‑Gttttt                                                                                                                                                                                                         >  1:568417/1‑250 (MQ=255)
                                                   cACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCACCGGCAAGGACGCCCCCTATACACGGGATCCATTTGGGTAAT‑GTTTTtctact                                                                                                                                                                                                    >  2:444842/1‑251 (MQ=255)
                                                          atatACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGAACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTG‑TATTTCTCCATCGGCAAGCACGACCTCTCTAAGCGGGATCCATTTGGCTAAT‑GTTTttctaactcgaaa                                                                                                                                                                                              >  2:22783/1‑242 (MQ=255)
                                                             taCCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCCGCGGCCAACACGTCCTCTTTTAACGGGATCAATTTTGCTAAT‑GTTTTTCTACTGCGAAAAat                                                                                                                                                                                           >  1:606731/1‑250 (MQ=255)
                                                              cccccAGCACGAACCTCCAACAACAAGATAGAGAATCTGCCTCGTAAGCGCGGGAGTACAATTTCAGTGACCACGACCAACAAACTCAATACCGCGACTCCGATGGCCTGCAGCTCCCACTGAAACAGGAAGCGGAATGGCAACAGACTGTGCAACTTCAGGTGGGGACGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAAGCAATTTGGCTAAT‑GTTCTACTACTGCGAAAAata‑t                                                                                                                                                                                        <  2:613115/250‑1 (MQ=255)
                                                                    gCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCt                                                                                                                                                                                    <  2:633100/249‑1 (MQ=255)
                                                                    gCACGTACCTCCAACAACCAGATGGCTAAGCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCTCGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGTCAATCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGc                                                                                                                                                                                  <  2:283404/251‑1 (MQ=255)
                                                                     cACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGAATCTTCAGGTGGGGCCGGCGGTTGGTATT‑CTCCCCAGGCAAGAACGTCCTCTTTAAGGGGGAAAAATTTGGCTTAT‑GTTTTATTAATGCGAAAAAAA‑TACCTGcc                                                                                                                                                                                 >  2:533776/1‑250 (MQ=255)
                                                                      aCGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCAc                                                                                                                                                                                <  1:437069/251‑1 (MQ=255)
                                                                             tCCAACAACCAAATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTTCAACGGCAAGCAAGGCCTCTTTAAACGGAATCAATTTTGGTAAT‑GTTTTACTACTGCGAAAAAAA‑TACCTGCACGCCAAAc                                                                                                                                                                         >  1:83969/1‑251 (MQ=255)
                                                                              ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACg                                                                                                                                                                        <  1:154131/251‑1 (MQ=255)
                                                                              ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACg                                                                                                                                                                        <  1:21197/251‑1 (MQ=255)
                                                                                  caacCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTTCTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACt                                                                                                                                                                    <  2:495001/251‑1 (MQ=255)
                                                                                  caacCAGATGGCTAATCTGCCGCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACt                                                                                                                                                                    <  2:107717/251‑1 (MQ=255)
                                                                                                    gCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTc                                                                                                                                                     <  1:68862/247‑1 (MQ=255)
                                                                                                       tCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAACGGCAAGCACGCCCCCTATAAACGGAATCCAATTGGCTAAT‑GTTTTATTACTGCGAAAAAAA‑TTCCATCACGCCAAACAGGCTTcctgggtcgggcccaccggg                                                                                                                                               >  2:379930/1‑233 (MQ=255)
                                                                                                          tAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAACGGCAAGCACGTCCTCTTTAAGCGGAATCAATTTTGCTAAT‑GGTTTACTACTGCGAAAAATA‑TCCCCGCACGCCAAACGGACTTCCCCGTCC‑GGCTCACCAAGGtt                                                                                                                                            >  1:556078/1‑250 (MQ=255)
                                                                                                           aaGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGAATCTTCAGGAGGGGCCGGCGGTTG‑TATTTCTCCAGCGGGAAGCCCCTCCTCCATAAACGGAATCAATTTGGCTAAT‑GTTTTACTACTTCGCAAAATA‑TCCCTGCCCGCCCAACGGACTTCCCCGTCC‑CGCTCACCCGGGGTac                                                                                                                                          >  1:61355/1‑251 (MQ=255)
                                                                                                                  gAGGTACATTTTCAGTGTCCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGAATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCCTCTTCAGGAGGGGCCGGCGGTTG‑TGTTTCTCCCGCGGCAACGAAGCCCTCTCTCAGCGGGATCTATTTTGCTAAT‑GTTTTACTCCTGCGCAAAATA‑TACCTGCAACCCAAACCGACTTTCCCGTCC‑GGcacccccggggtacaccgcgc                                                                                                                                   >  2:297974/1‑231 (MQ=255)
                                                                                                                     gTACATGTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTAt                                                                                                                                <  1:613997/251‑1 (MQ=255)
                                                                                                                             tCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTA‑TATTTTTTCCGCGGCAAGCACGGACTCTTTTAGGGGAAACAATTTTGGTAAA‑ATTTTACTACTGCGAAAAATA‑TTCCTGCACGACAAACGGAATTTCCGGTCC‑GGCCCAACAAGGGTACGAAGCGCtcactccaaa                                                                                                                         >  1:682225/1‑241 (MQ=255)
                                                                                                                                    ccACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCgg                                                                                                                <  1:49858/251‑1 (MQ=255)
                                                                                                                                        gACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCAGCGGCAAGCAAGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTTCTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGGTAAGCAGCGCTATCTCCAAA‑GCTTTCGGTac                                                                                                             >  1:454229/1‑250 (MQ=255)
                                                                                                                                         aCCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TGTTTTTCCAGCGGGAAGCACGGCCTCTATAAGCGGAATCAATTTTGCTAAT‑GTTTTAATACTGGGGAAAATA‑TCCCTTCACCCCAAAAGGACTTCCCCGGCC‑CGCTCACCAGGGGTACGACGCGCCATCTCAAAA‑GCTTTCGGTaca                                                                                                            >  2:613958/1‑250 (MQ=255)
                                                                                                                                         aCCAACATACTCATTTCAGCTACTCCGATGGCCCGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑ATTTTTTTCAGAGGCAAGCACGTCCTCTATAAAAGGAATCAATTTAGCTAAT‑GTTTTACTACTTCGAAAAATA‑TCCCTGCCAGCCCAACTGACTTTCCAGTCCTGGC‑CACCACGGTTACAGCCCCCTATCTCAAAA‑GCTTTCGGTaacc                                                                                                           >  1:584978/1‑248 (MQ=255)
                                                                                                                                               aTACTCATTTCAGCTACGTCGATGGCCTGCATCTCCCAGGGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGTGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCGAAACGGACTTCCCGGTCC‑GGCGCACCAGGGTTACGACGCGCTATCTCAAAA‑GCGTTCGGTAAAACAgtgt                                                                                                     <  2:79273/251‑1 (MQ=255)
                                                                                                                                                       ttCAGCGACTCCGCTGCCCAGAATCTCCCAGTGAAGAAGGAGTCGGAATGGTACCAGGCTGTCCATCTTCAGGTGGTCTCGGCATTTG‑TATTTCTCCAGCTTCAAGCTCGTCTTCTATAAGCGGAACCAATTTGGCTAAT‑GTTTTAATATTGCGAAAAATA‑TCACTGCTCGCCAAACGTTCTTCCCGGCCC‑GGCTAACCAGGCTTACGACGCGCTATCTCGAAA‑GCTTTCGGTACAACAGTGTCAAAttt                                                                                              <  1:484881/250‑1 (MQ=255)
                                                                                                                                                         cAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TGTTTTTCCAGCGGCAAGCACGTCCTCTATAAACGGGATCAATTTTGCTAAT‑GTTTTTATACTGCGAAAAATA‑TTCCTGCCCCCCAAACGGACTTTCCGGGCC‑GGCCCACCCGGGTTACGACCCGCTATTTCCAAA‑CCTTTCGGTACAACAGAGGCAAATTTAat                                                                                           >  1:64108/1‑251 (MQ=255)
                                                                                                                                                                        ccTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATc                                                                            <  1:143998/251‑1 (MQ=255)
                                                                                                                                                                        ccTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGGGGGGCCGGCGGTTT‑GATTTTTCCAGCGGCCAGCCCGGCCTCTTTAAACGGGATCAATTTTGGTAATTGTTTT‑CTACTGCGGAAAAAA‑ACCCTGCACGCCCAAAGGACTTCCCGGGCC‑GGCCCCCCAGGGGTAACACGCGCCATCTCCAAA‑ACTTTTGGGTCAACAGGGTCCAATTTAAAGATGTCTGCAATATc                                                                            >  2:209372/1‑251 (MQ=255)
                                                                                                                                                                             aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAATTTGGGTAAT‑GTTTTTCTCCTGCGCAAAATA‑TCCCTGCACGCCCCACGGACTTCCCGGTCC‑GGCTCACCAGGGTTAACACCCGCTATCTCAAAA‑GCTTTTGGTGCAACAGTGTCAAATTTAATTATGTCTGCCATATCGTaa                                                                        >  1:350482/1‑250 (MQ=255)
                                                                                                                                                                             aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGGGGGGCCGGCGGTTG‑TATTTCTCCAAAGGCAAACCCGTGCTCTTTAAACGGAATCAATTTGGCTAAT‑GTTTTACTACTGTGGAAAATA‑TCCCTGCACGCCAAACGGGATTCCCGGTCC‑GGCTAACCAGGGGTAACACGCGGCATCTCAAAAAGCTT‑CTGGACAACCGTGGCAAATTTTATAATTTCTGCGATATCGTaa                                                                        >  1:128574/1‑250 (MQ=255)
                                                                                                                                                                                 cccAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCGTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑GCCCTGCACGCCAAACTGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACa                                                                     <  2:98706/249‑1 (MQ=255)
                                                                                                                                                                                 cccAGTGAAACAGGAAGCGGAATGACAACAGGCTGTGCATCTGCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTTCTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GTCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTCTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACaaa                                                                   <  1:377166/251‑1 (MQ=255)
                                                                                                                                                                                   cAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATGTGGCTAAT‑GTTGTACTACTGCGAAAAATA‑TCCCTGCACGGCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCaaaa                                                                                                                         <  1:683502/196‑1 (MQ=255)
                                                                                                                                                                                    aGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCCACGGCAAGCACGTCCTCTATAAACGGGATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGGCC‑GGCTCACCAGGGTTACGACGCGCTATCTCCAAA‑ACTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACaaa                                                                   >  1:129081/1‑248 (MQ=255)
                                                                                                                                                                                       gAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCAGAGGGAAGCCCCTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTTCTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGGCC‑GGCTCACCAGGGGTACGGCGCGCCATCTCAAAA‑GCTTTCGGTCCAACCGTGTCAAATTTAATGATGTCTGCCATATCGTAAACAAAGGaaa                                                              >  1:465823/1‑250 (MQ=255)
                                                                                                                                                                                        aaaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCCGCGGGAAGCACCGCCTCTATAAGCGGAATCAATTTTGCTAAT‑GTTTTACTCCTGCGAAAATAA‑TCCCTGCACGCCCAACGGACTTTCCCGTCC‑CGCTCACCCGGGTTACGACGCGCTATCTCAAAA‑GGTTTCGGTACAACAGTGTCCAATTTAATGATGTCTGCCatagacaaaacaaaagaa                                                               >  2:185305/1‑233 (MQ=255)
                                                                                                                                                                                            aGGAAGCGGAATGGCAACAGGCTGTTCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATAGCGTAAACAAAGGAAAGAGGCt                                                        <  1:599002/251‑1 (MQ=255)
                                                                                                                                                                                            aGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTTTCCAGCGGCCAGCACGTCCTCTATAAACGGGATCAATTTGGCTAAT‑GTTTTTCTTCTGCGAAAAAAA‑ACCCTGGACGCCAAACGGGCTTCCCGGGCC‑GGCCCCCCAGGGTTTCGACGCGCTATCTCCCAA‑GCTTTCGGTACCACAGGGTCCAATTTAATTGTGTCTTCAATATTCTAAACAAAgggaaaagggt                                                        >  2:310855/1‑242 (MQ=255)
                                                                                                                                                                                            aGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAACACGTCCTCTCTAAGCGGGATCAATTTGGCTAAT‑GTTTTTCTACTGCGAAAAATA‑TCCCTGCACGCCCAACGGACTTCCCGGTCC‑GGCTCACCCGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAAAGATGTCTGCCATATCCTAAAAAAAGGAAAGAgggt                                                        >  1:614076/1‑249 (MQ=255)
                                                                                                                                                                                                aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTcc                                                    <  1:325829/251‑1 (MQ=255)
                                                                                                                                                                                                 gCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TCTTTTTCCAGCAGGAAGCCCGGTCTCTATAAGGGGGAACAAATTGGCTAAT‑GTTTTACTACTGCGGAAAATAATCCCT‑CCCCCCAAAACGACTTCCCCGTCC‑GGGTCCCCCGGGGTACGCAGCCCCTTCTCCACA‑ACTTTTGGTACACAAAAGTCAAATTTAATGATGTATGCACTAATGTAAACAAAGGaaaaggggttttct                                                   >  2:636526/1‑240 (MQ=255)
                                                                                                                                                                                                    gAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCaa                                                                                                                           <  1:491025/177‑1 (MQ=255)
                                                                                                                                                                                                        ggCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAACGGCAAGCACGTCCCCTATAAGCGGGATCAATTTGGCTAAA‑GTTTTTCTACTTCGAAAAATA‑TCCCTGCCCGCCCAACGGACTTCCCGGTCC‑CGCTCACCAGGGGTACGACGCCCCATCTCCAAA‑GCTTTCGGTACAACAATGGCAAATTTAATGATGTCTGCCATATCGTAAACAAAGGAAAAAAGCTTTCCttttttcc                                            >  1:273424/1‑244 (MQ=255)
                                                                                                                                                                                                         gCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑AATTTTTCCACCGACAAACACATCCCCTTTAAGCGGGATCAATTTTGCCAAT‑GTTTTAATACTTCGAAAAAAA‑TCCCTTCACCCCAAAACGACTTCCCCGTCC‑GGCTCACCCAGGTTACGACCCCCTTTCTCAAAA‑ACTTTCGGTTCAACAGTGTCAAATTTTATGATTTCTTCAATATCTTTAACAAAGGAAAAATGCTTTCCTTTATGccc                                           >  1:67742/1‑249 (MQ=255)
                                                                                                                                                                                                                   gtgCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCaa                                   <  1:636526/249‑1 (MQ=255)
                                                                                                                                                                                                                   gtgCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAAt                                  <  1:379930/250‑1 (MQ=255)
                                                                                                                                                                                                                           tCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTg                             <  1:310855/247‑1 (MQ=255)
                                                                                                                                                                                                                                tGGGGCCGGCGGTTG‑TTTTTTTCCCGCGGCAAGAACATCCTCTTTAAACGGAATCAATTTTGTTAAT‑GTTTTTCTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGGCTTCCCGGTCC‑GGCGCCCCAGGGGTACGACGCGCTATCTCACAA‑GATTTTGGGACAAAAATGGCAAAATTAATGATTTCTGCAATATCGTAAAAAAAGGAAAGAGGCATTCCTGTATGCACAAACCCAAAAGCGGGTGCAAAAc                    >  2:285346/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccGGCGGTTG‑TATTGCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCTCGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACgcgc                                                                                                                                   <  1:337645/136‑1 (MQ=255)
                                                                                                                                                                                                                                        gaggTTG‑TATTTATCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GATTTACTACTGCGAAACATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCA‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGACAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCGTTCCTATATGCACAAACCAAATAGCTGGAGCATAACCAGCTGc             >  2:578263/3‑250 (MQ=255)
                                                                                                                                                                                                                                           gTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGt         <  2:350482/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    tctcCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgccg  <  1:41706/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
CCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTG‑TATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAT‑GTTTTACTACTGCGAAAAATA‑TCCCTGCACGCCAAACGGACTTCCCGGTCC‑GGCTCACCAGGGTTACGACGCGCTATCTCAAAA‑GCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  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: ‑