Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,9470CT28.3% 79.8 / 12.9 48E276K (GAA→AAA) cas1multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (4/29);  new base T (13/0);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 2.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCG  >  NC_000913/2878700‑2879192
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
ccGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTcc                                                                                                                                                                                                                                                           >  1:83449/1‑251 (MQ=255)
               tttGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGAGGGGCCGGCGGTTGTATTTTTTCATCGGCAAGAACGTc                                                                                                                                                                                                                                            >  1:227875/1‑251 (MQ=255)
                 tGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTTTTTTTCCCGCGGCaaaccccagcc                                                                                                                                                                                                                                          >  2:463144/1‑242 (MQ=255)
                  gTTCCCAGATCATTTCACGAATTTTTTCGGATACATCACCTACATACACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCct                                                                                                                                                                                                                                          <  2:289058/250‑1 (MQ=255)
                    tCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGAGGCGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTAAGGTGGGGCCGGCGGTTGGATTTATcccccgcgcagcacgctaccc                                                                                                                                                                                                                                        >  2:541885/1‑231 (MQ=255)
                               ttCACGAATTTTTGCGGATACAGCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCAGCCGCGTAAGCGCGGAGTTCCATTTTCAGTGACCACGACCAACATACTCATGTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTGCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAAt                                                                                                                                                                                                                             <  2:481575/250‑1 (MQ=255)
                                tCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCAGGCGGTTGGATTTTTCCAACGGCacacaccacctctcatcgccgaatc                                                                                                                                                                                                                            >  2:86964/1‑226 (MQ=255)
                                tCACGAATTTTTGCGGATACATCACCTACATAAACCCCTGCACGTACCGCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATc                                                                                                                                                                                                                            <  2:164483/250‑1 (MQ=255)
                                                  aCATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGGCCTCTATAAGCGGGATAACTTTGGCTAATGtttttt                                                                                                                                                                                                         >  1:63914/1‑249 (MQ=255)
                                                   cATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtact                                                                                                                                                                                                        <  1:72769/251‑1 (MQ=255)
                                                   cATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtact                                                                                                                                                                                                        <  1:463144/251‑1 (MQ=255)
                                                    aTCACCTACATATACCCCTGCCCGTACCTCCAACACCCAGATGGCTCATCTGCCTCGTAAGCGCGGAGGTACATTTTCATTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCAGCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtacta                                                                                                                                                                                                       <  2:381685/251‑1 (MQ=255)
                                                       cccTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCATCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTg                                                                                                                                                                                                    <  1:258477/250‑1 (MQ=255)
                                                           aCATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACGCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGaa                                                                                                                                                                                                <  2:124573/251‑1 (MQ=255)
                                                                   cccTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATCTcc                                                                                                                                                                                        <  2:83449/251‑1 (MQ=255)
                                                                    ccTGCACGTACCTACAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACGCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGTGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATggtt                                                                                                                                                                                                            <  2:157172/230‑4 (MQ=255)
                                                                       gCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTTCATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAATAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTg                                                                                                                                                                                     <  2:296402/250‑1 (MQ=255)
                                                                          cGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACg                                                                                                                                                                                 <  2:312454/251‑1 (MQ=255)
                                                                             ccccccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTGCAGTGACCACGACCAACATACTCATGTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAGGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑a                                                                                                                                                                             <  2:124232/247‑1 (MQ=255)
                                                                                   aacaacCCGATGGCTAATCTGCCTCGGAAGCGCGGAGGTACATTTTCAGTTACCACGTCCAACATACTCAGTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCGCCAGCGGCAAGCACGTCCTCTATAATCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACgta                                                                                                                                                                       <  2:20409/251‑3 (MQ=255)
                                                                                        ccAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACCGCCTCTATAAAAGGAATCAATTTTGCTTATGTTTTACTACTTCGAAAAAAATTCCCGCCCCCC‑AAACGGA‑C‑tttc                                                                                                                                                                >  1:370185/1‑249 (MQ=255)
                                                                                            aTGGCTAATTTGCCGCGTAAGCGCGGAGGTACATTTTCAGTTACCTCGACCAACATACTCATTGCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGTTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACTGA‑C‑TTCCCGGt                                                                                                                                                            <  1:383390/251‑1 (MQ=255)
                                                                                             gggCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCATGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTc                                                                                                                                                           <  1:103156/250‑1 (MQ=255)
                                                                                               gCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCg                                                                                                                                                         <  1:363098/251‑1 (MQ=255)
                                                                                                        ccTCGTAAGCGCGGGAGTACATTTTCATTGACCTAGACCAGCATGCTCATTTCAGCTACTCCGATGGCCGGCCTCTCCCAGTGAAACAGGTAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCAGGTGATTGTATTTCTCCAGCGGCAAGCACGTCCCCTATAAGCGGAATCAATTTGGCTAATGTTTTACAACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGACCGGCTCACCA‑g                                                                                                                                               <  2:26272/251‑1 (MQ=255)
                                                                                                         cTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTTTTTTTCCAGCGGGAAACACGTTCTCTATAAAAGGAGTCAATTTGTGTCATGTTTTTCTACTGCGAAAAAAATACCCGCACGCC‑CAACAGA‑C‑CTCCCgggcccggtcacccgg                                                                                                                                               >  1:235591/1‑237 (MQ=255)
                                                                                                               aGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTCTCCAGCGGCGAACACCTCCCCTTTAAGCGGGAACAATTTGGGTTAAGTTTTTATACTTCGAAAAATATCCCTGCCCGCC‑AAACGGAAC‑TTCCCGGTC‑GGCTCACCC‑GGG‑GTacg                                                                                                                                       >  1:324307/1‑251 (MQ=255)
                                                                                                                   cgGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCACCGGCAAACACGTCCTCCATAAGCGGAATCAATTTTGTTCATTTTTTTTTAATACGGAAAATATCCCTGCCCCCC‑AAAAGGA‑CATTCCCgttcgcctcaaccagggtacgaacc                                                                                                                                      >  1:127401/1‑227 (MQ=255)
                                                                                                                                       ccACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCgg                                                                                                              <  1:85757/251‑1 (MQ=255)
                                                                                                                                                   tACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCCAACACGTCCTCTATAAGCGGAATCAAATTTGCTAATGTTTTACTTCTGCGAAAAAAATTCCTGCACGCC‑AAAAGGA‑C‑TTCCCGGGCCCGACCACCC‑GGG‑GTACGACGC‑GCTCTCTCAAAAGCTTTCTGTacaaacgagtc                                                                                                  >  1:39399/1‑244 (MQ=255)
                                                                                                                                                                tACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTTTTTTTCCAGCGGCCAGGACGTCCTCTTTAAGCGGGAACAATTTGGGTAATGTTTTACTTCTGGGCAAAATATCCCTGCACGCC‑AAACGGG‑C‑TTCCCCGTCCGGCCCACCC‑GGG‑GTAACACGC‑GCTATCTCCAAAACTTTCGGTACAAAAATGG‑CCAATTTAatggtt                                                                                    >  1:473952/1‑248 (MQ=255)
                                                                                                                                                                                gtctcttccgatctcAGGAAGCGGAATGGCAACAGGCTGTGCATCTTGAGGTGGGGCCGGCGGTTTTATTTCTCCAGCGGCAAGCACGTCCTCGATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAAATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCCCCA‑GGG‑GTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTaaa                                                                    <  1:374582/237‑1 (MQ=255)
                                                                                                                                                                                aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTTTTTGTCCCACGGCAAGCAACTCCTCTCTAAGCCGAATCCATTTTGCTAATGGTTTTCTTATACTCAAAAAATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCGCACAA‑GGG‑ATAAGCAGC‑CCCATCCCATAAGATTTCGGTACAACAGTGT‑CAAATATTATGTTGTCTGCAATAACGTaa                                                                     >  2:176190/1‑250 (MQ=255)
                                                                                                                                                                                  ctcCCAGTTAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCCGCCGGAAGCACGTCCTCTATTAGCGGAATCAATTTTGCTCATGTTTTTATTCTGCGGAAAATATCCCCGCACGCCCAAACG‑A‑A‑TTTCCGGGCCGGCTCACCA‑GGGGTTA‑GACGCCGC‑ATCTCAAAAGCTTTCGGTACAACAGTGT‑CCAATTTAATGGTATCTTCCATATCGTaaaa                                                                   >  1:258787/1‑249 (MQ=255)
                                                                                                                                                                                      cAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCAAGCACGTCCTCTATAAGCGGGAACCATTTTGCTAATTTTTTACTACTGCGGAAAATATTCCTGCACCCC‑AAACCGA‑C‑TTCCCGGGCCGGCTCCCCACGGG‑TT‑CGACGC‑GCTATCTCAAAAACTTTCCGTACAACAGAGT‑CCAATTTTATTATGTCTGccaaattgtaaaacaa                                                                >  1:201940/1‑234 (MQ=255)
                                                                                                                                                                                            aaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAagag                                                        <  2:201940/251‑1 (MQ=255)
                                                                                                                                                                                                ggAAGCGGAATGGCAACCGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCATCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TGCCCGGTCCGGCTCACCA‑GTG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTTT‑CAAATTTAATGATGTCTTCAATATCGTAAACAAAGGATAGAGGCtt                                                    <  2:470516/251‑1 (MQ=255)
                                                                                                                                                                                                   aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTTGGGCCGGCGGTTGTATTTCGCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCTGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTcc                                                 <  2:389925/251‑1 (MQ=255)
                                                                                                                                                                                                         aTGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTATATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTAACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAAc                                                                   <  1:496587/227‑1 (MQ=255)
                                                                                                                                                                                                                cAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCCGCAAGCAAGTCCTCTATAAACGGGATCAATTTTGCTAAAGTTTTACTACTGAGGCAAATATCCCTCCCCGCC‑AAACGGA‑A‑TTTCCGGTCCCGC‑CCCCACGGG‑TTACGACGC‑GCTAACTCCAAAGCCTTTGGTACAAAAAGGTGCAAATTTA‑TGATGTTTGCAATATCGTAAAAAAAAGGAAGAAGCGTTCCTGTATgaaaaaa                                     >  1:261688/1‑244 (MQ=255)
                                                                                                                                                                                                                      gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGATGCTTTCCTGTATGCACAAACCCaa                                <  1:433477/249‑1 (MQ=255)
                                                                                                                                                                                                                      gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATa                              <  2:114129/251‑1 (MQ=255)
                                                                                                                                                                                                                        gtagcTTCAGGGGGGGCCGGCGGGTGTATTTCTCCAGCGGCAAGCACGTCCTCGATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAGATCCCTGCACGCC‑AAACGGA‑C‑TGCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTGTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAg                             <  2:520408/246‑1 (MQ=255)
                                                                                                                                                                                                                          aTCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCAAACACGTCCTCTATAAGAGGGATCAATTTGGGTAAAGTTTTTCTAATTCGAAAAATATCCCCTCACGCC‑AAACGGG‑C‑TTCCCCGTCCGGCTCACCC‑GGG‑GTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATATAATAATGTCTGCACTATCGTAAACAAAAGAAAGCGGCTTTCCTGTATGCACCAAACCAATAGCTg                          >  1:34504/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ggTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATa                   <  1:105600/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTg           <  2:63914/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAATTTGGCTCATGTTTTTCTACTGCGAAAAAAATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGGCCGGCTCACCC‑GGG‑TTACGACGC‑GCTATCTCCAAAGCTTTCGGGCCAACAGGGT‑CAAATTTAATGATTTCTGCACTATCCTAAAACAAGGAAAGAGGCATTCCTGTTACCCCAAACCCAATACCGGGTGCATAACCAGCTGc          >  2:247943/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcc  >  2:272084/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
CCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCC‑AAACGGA‑C‑TTCCCGGTCCGGCTCACCA‑GGG‑TTACGACGC‑GCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCG  >  NC_000913/2878700‑2879192

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