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.

CCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCG  >  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)
                                                                             ccccccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTGCAGTGACCACGACCAACATACTCATGTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAGGTTTTACTACTGCGAAAAATATCCCTGCACGCCa                                                                                                                                                                             <  2:124232/247‑1 (MQ=255)
                                                                                   aacaacCCGATGGCTAATCTGCCTCGGAAGCGCGGAGGTACATTTTCAGTTACCACGTCCAACATACTCAGTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCGCCAGCGGCAAGCACGTCCTCTATAATCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACgta                                                                                                                                                                       <  2:20409/251‑3 (MQ=255)
                                                                                        ccAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACCGCCTCTATAAAAGGAATCAATTTTGCTTATGTTTTACTACTTCGAAAAAAATTCCCGCCCCCCAAACGGACtttc                                                                                                                                                                >  1:370185/1‑249 (MQ=255)
                                                                                            aTGGCTAATTTGCCGCGTAAGCGCGGAGGTACATTTTCAGTTACCTCGACCAACATACTCATTGCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGTTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACTGACTTCCCGGt                                                                                                                                                            <  1:383390/251‑1 (MQ=255)
                                                                                             gggCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCATGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTc                                                                                                                                                           <  1:103156/250‑1 (MQ=255)
                                                                                               gCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCg                                                                                                                                                         <  1:363098/251‑1 (MQ=255)
                                                                                                        ccTCGTAAGCGCGGGAGTACATTTTCATTGACCTAGACCAGCATGCTCATTTCAGCTACTCCGATGGCCGGCCTCTCCCAGTGAAACAGGTAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCAGGTGATTGTATTTCTCCAGCGGCAAGCACGTCCCCTATAAGCGGAATCAATTTGGCTAATGTTTTACAACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGACCGGCTCACCAg                                                                                                                                               <  2:26272/251‑1 (MQ=255)
                                                                                                         cTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTTTTTTTCCAGCGGGAAACACGTTCTCTATAAAAGGAGTCAATTTGTGTCATGTTTTTCTACTGCGAAAAAAATACCCGCACGCCCAACAGACCTCCCgggcccggtcacccgg                                                                                                                                               >  1:235591/1‑237 (MQ=255)
                                                                                                               aGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTCTCCAGCGGCGAACACCTCCCCTTTAAGCGGGAACAATTTGGGTTAAGTTTTTATACTTCGAAAAATATCCCTGCCCGCCAAACGGAACTTCCCGGTCGGCTCACCCGGGGTacg                                                                                                                                       >  1:324307/1‑251 (MQ=255)
                                                                                                                   cgGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCACCGGCAAACACGTCCTCCATAAGCGGAATCAATTTTGTTCATTTTTTTTTAATACGGAAAATATCCCTGCCCCCCAAAAGGACATTCCCgttcgcctcaaccagggtacgaacc                                                                                                                                      >  1:127401/1‑227 (MQ=255)
                                                                                                                                       ccACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCgg                                                                                                              <  1:85757/251‑1 (MQ=255)
                                                                                                                                                   tACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCCAACACGTCCTCTATAAGCGGAATCAAATTTGCTAATGTTTTACTTCTGCGAAAAAAATTCCTGCACGCCAAAAGGACTTCCCGGGCCCGACCACCCGGGGTACGACGCGCTCTCTCAAAAGCTTTCTGTacaaacgagtc                                                                                                  >  1:39399/1‑244 (MQ=255)
                                                                                                                                                                tACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTTTTTTTCCAGCGGCCAGGACGTCCTCTTTAAGCGGGAACAATTTGGGTAATGTTTTACTTCTGGGCAAAATATCCCTGCACGCCAAACGGGCTTCCCCGTCCGGCCCACCCGGGGTAACACGCGCTATCTCCAAAACTTTCGGTACAAAAATGGCCAATTTAatggtt                                                                                    >  1:473952/1‑248 (MQ=255)
                                                                                                                                                                                gtctcttccgatctcAGGAAGCGGAATGGCAACAGGCTGTGCATCTTGAGGTGGGGCCGGCGGTTTTATTTCTCCAGCGGCAAGCACGTCCTCGATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAAATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCCCCAGGGGTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTaaa                                                                    <  1:374582/237‑1 (MQ=255)
                                                                                                                                                                                aTCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTTTTTGTCCCACGGCAAGCAACTCCTCTCTAAGCCGAATCCATTTTGCTAATGGTTTTCTTATACTCAAAAAATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCGCACAAGGGATAAGCAGCCCCATCCCATAAGATTTCGGTACAACAGTGTCAAATATTATGTTGTCTGCAATAACGTaa                                                                     >  2:176190/1‑250 (MQ=255)
                                                                                                                                                                                  ctcCCAGTTAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCCGCCGGAAGCACGTCCTCTATTAGCGGAATCAATTTTGCTCATGTTTTTATTCTGCGGAAAATATCCCCGCACGCCCAAACGAATTTCCGGGCCGGCTCACCAGGGGTTAGACGCCGCATCTCAAAAGCTTTCGGTACAACAGTGTCCAATTTAATGGTATCTTCCATATCGTaaaa                                                                   >  1:258787/1‑249 (MQ=255)
                                                                                                                                                                                      cAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCAAGCACGTCCTCTATAAGCGGGAACCATTTTGCTAATTTTTTACTACTGCGGAAAATATTCCTGCACCCCAAACCGACTTCCCGGGCCGGCTCCCCACGGGTTCGACGCGCTATCTCAAAAACTTTCCGTACAACAGAGTCCAATTTTATTATGTCTGccaaattgtaaaacaa                                                                >  1:201940/1‑234 (MQ=255)
                                                                                                                                                                                            aaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAagag                                                        <  2:201940/251‑1 (MQ=255)
                                                                                                                                                                                                ggAAGCGGAATGGCAACCGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCATCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTGCCCGGTCCGGCTCACCAGTGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTTTCAAATTTAATGATGTCTTCAATATCGTAAACAAAGGATAGAGGCtt                                                    <  2:470516/251‑1 (MQ=255)
                                                                                                                                                                                                   aGCGGAATGGCAACAGGCTGTGCATCTTCAGGTTGGGCCGGCGGTTGTATTTCGCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCTGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTcc                                                 <  2:389925/251‑1 (MQ=255)
                                                                                                                                                                                                         aTGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTATATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTAACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAAc                                                                   <  1:496587/227‑1 (MQ=255)
                                                                                                                                                                                                                cAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCCGCAAGCAAGTCCTCTATAAACGGGATCAATTTTGCTAAAGTTTTACTACTGAGGCAAATATCCCTCCCCGCCAAACGGAATTTCCGGTCCCGCCCCCACGGGTTACGACGCGCTAACTCCAAAGCCTTTGGTACAAAAAGGTGCAAATTTATGATGTTTGCAATATCGTAAAAAAAAGGAAGAAGCGTTCCTGTATgaaaaaa                                     >  1:261688/1‑244 (MQ=255)
                                                                                                                                                                                                                      gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGATGCTTTCCTGTATGCACAAACCCaa                                <  1:433477/249‑1 (MQ=255)
                                                                                                                                                                                                                      gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATa                              <  2:114129/251‑1 (MQ=255)
                                                                                                                                                                                                                        gtagcTTCAGGGGGGGCCGGCGGGTGTATTTCTCCAGCGGCAAGCACGTCCTCGATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAGATCCCTGCACGCCAAACGGACTGCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTGTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAg                             <  2:520408/246‑1 (MQ=255)
                                                                                                                                                                                                                          aTCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCCGCGGCAAACACGTCCTCTATAAGAGGGATCAATTTGGGTAAAGTTTTTCTAATTCGAAAAATATCCCCTCACGCCAAACGGGCTTCCCCGTCCGGCTCACCCGGGGTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATATAATAATGTCTGCACTATCGTAAACAAAAGAAAGCGGCTTTCCTGTATGCACCAAACCAATAGCTg                          >  1:34504/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ggTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATa                   <  1:105600/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTg           <  2:63914/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAATTTGGCTCATGTTTTTCTACTGCGAAAAAAATCCCTGCACGCCAAACGGACTTCCCGGGCCGGCTCACCCGGGTTACGACGCGCTATCTCCAAAGCTTTCGGGCCAACAGGGTCAAATTTAATGATTTCTGCACTATCCTAAAACAAGGAAAGAGGCATTCCTGTTACCCCAAACCCAATACCGGGTGCATAACCAGCTGc          >  2:247943/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcc  >  2:272084/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
CCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCG  >  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: