Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,9470CT26.1% 88.0 / 13.9 47E276K (GAA→AAA) cas1multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (5/29);  new base T (12/0);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 1.59e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  NC_000913/2878699‑2879195
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
accGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCTGATACAGCTCCTACATATACCCCTGCACGTACCTCCAACAGCCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATGGCAGCTTCTCCGATGGCCGGCATCTCCCAGTGAAACAGGAAGAGGAATGGCAACGGGCGGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTtct                                                                                                                                                                                                                                                            <  2:107415/249‑1 (MQ=255)
       gTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGACGCGGAATGGCAACAGGTTGTGCATCTTCAGGTGGGGCCGGCGGCTGTTTTtctaccacggc                                                                                                                                                                                                                                                    >  2:248316/1‑243 (MQ=255)
            gCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTGTTTCTcccccggccagcc                                                                                                                                                                                                                                               >  2:251579/1‑240 (MQ=255)
             cTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCCACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGCACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTTCCAGCGGCCAGCAc                                                                                                                                                                                                                                             >  1:334344/1‑251 (MQ=255)
                    ttCCCACAGCATTTCACGAAGTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAACCCTCCTCTTAAGCGCGGAGGTGCATTTTCATTGACCACGACCAACACACTCATCTCAGCTACTCCGATGGCCTGCATCGCCCAGTGAAACATGAAGCTGAATGGCAGCAGGCTTTGCACCTTCAGTTTGGGCCGGCGGTTGTATGTCGCCTGCGTCAAGCATGTCctc                                                                                                                                                                                                                                       <  2:526330/250‑1 (MQ=255)
                      cccAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTAGTTCTCCAGCGGCAAGCACGTCCTCtat                                                                                                                                                                                                                                    <  1:422848/251‑1 (MQ=255)
                                 tCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATc                                                                                                                                                                                                                          <  2:131375/250‑1 (MQ=255)
                                    ccaaTTTTTGCGGATTCATCCCCTACCAATACCCCCGCACCTTCCTCCAACAACCAAAATGCTTATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAAtt                                                                                                                                                                                                                     <  1:432225/249‑1 (MQ=255)
                                         ttttGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCTGAATCAATTTGGc                                                                                                                                                                                                                 <  1:457693/250‑1 (MQ=255)
                                         ttttGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCAGCGGGAAGCACGCCCTCTATAAGCGGAAACAAATTGGCt                                                                                                                                                                                                                >  1:51051/1‑251 (MQ=255)
                                          tttGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCACCGGGAAAAACGTCCACCATTAGCGGGATCAAGTTTggt                                                                                                                                                                                                                 >  2:149135/1‑249 (MQ=255)
                                               ggATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGtt                                                                                                                                                                                                          <  1:261393/251‑1 (MQ=255)
                                                   aCATCACCTACATACACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTGCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGGGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGGCCTCTATAAGCGGAATCAATTTGGCTAATGTTTta                                                                                                                                                                                                       <  2:96235/250‑1 (MQ=255)
                                                          cTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCCGCGGCAAGCACGTCCTCTCTAAACGGAATCAATTTGGGTAATGTTTTACTACTGCg                                                                                                                                                                                               >  1:149211/1‑251 (MQ=255)
                                                             cATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTGTCCAGCGGCAAGCACGGTCGCTATAAACGGGATCAATTTGGGTAATGtttttctcgtgctca                                                                                                                                                                                             >  1:315983/1‑239 (MQ=255)
                                                                  aCCCCTGCAAGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGAGGGGCCGGCGGTTGGATTTTACCCGCGGGCAACACACCCTCCATAAGCGGAAACAATTTGGCTAATGGTTTACAACTGCGacaaacat                                                                                                                                                                                      >  2:418780/1‑248 (MQ=255)
                                                                      cTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCCGCGGCAAGCACGTCCCCTATAAGCGGAATCa                                                                                                                                                                                                                         >  1:208103/1‑214 (MQ=255)
                                                                      cTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCa                                                                                                                                                                                                                         <  2:208103/214‑1 (MQ=255)
                                                                              aCCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGcc                                                                                                                                                                           <  1:299282/250‑1 (MQ=255)
                                                                                    aacaacCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGGACATTGTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGa                                                                                                                                                                    <  2:130398/251‑1 (MQ=255)
                                                                                        acCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTc                                                                                                                                                                <  1:426832/251‑1 (MQ=255)
                                                                                                        gCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCAc                                                                                                                                                  <  1:312594/249‑1 (MQ=255)
                                                                                                             gTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAATTTTGCTAATGTTTTTCTACTGCGAAAAATATACCTGCACGCCAAACGGACTTCCCCGTCCGGcccccccgggtt                                                                                                                                           >  1:452713/1‑240 (MQ=255)
                                                                                                               aaGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTTCACCGGCCAGCACGTACTCTTTAAACGGGAGCAATTTGGGTAATGTTTTTTTACTGTGCAAAAATTCCCTGTACGCCCAACGGGTTTACCCGTCCGGCGCAccccgggtac                                                                                                                                         >  2:307856/1‑243 (MQ=255)
                                                                                                               aaGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTa                                                                                                                                          <  2:149211/250‑1 (MQ=255)
                                                                                                                aGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGGGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGCGTTacg                                                                                                                                        <  2:524046/251‑1 (MQ=255)
                                                                                                                                   aGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCACCGGCGCGCACGTCCCCTATAAGCGGAATTAAATTTGCTTATGTTTTACTACTGCGAAAAATATCCCTGCACGCCCAACAGACTTCTCGGTCCGGCTCACCCGGGGTACGACGCGCCATCTCcaaaact                                                                                                                     >  1:327681/1‑248 (MQ=255)
                                                                                                                                   aGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGGCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCt                                                                                                                     <  2:364673/251‑1 (MQ=255)
                                                                                                                                      gACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTc                                                                                                                  <  2:95542/251‑1 (MQ=255)
                                                                                                                                       aCCACGACCAACATACTCATTTCAGCTACTCCGATGGCATGCCTCTCCCAGAGAAACAGGAAGCGGAATGGCAACAGGATGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGGAAGCCCGTCCTCTATAAGAGGAATCAATTTGGCTACTGTTTTTCTACTGCGCAAAATATACCTGCACGCCAAAAGGACTTCCCCGTCCGGGTCCCCCCGGTTACGACGCGCTATCTCTAAAgacttcc                                                                                                                 >  1:63070/1‑245 (MQ=255)
                                                                                                                                            gACCAACATACTCATTTCAGCTACTCCTATGGACTGCATCACCCAGGGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCGCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAGATCCCTGCACGCCAAACGGCCTTCCCGGTCCGGCGCACCAGGGTGACGACGCGCTATCACAAAAGCTTTCGGtgca                                                                                                            <  2:45549/251‑4 (MQ=255)
                                                                                                                                                    aaCTCATGTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTGTTACTACTGCGAAAAATATCCCTGCGCGCCAATCGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCGCTCGGTACAACAGTGTc                                                                                                    <  2:334344/250‑1 (MQ=255)
                                                                                                                                                             cAGCAACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTGCCCTGTCCGGCTCACCAGGGTTACGACGCGCTATCGCAAAAGCTTTCGGTACAACAGTTTCAAATTTAat                                                                                           <  2:219308/251‑1 (MQ=255)
                                                                                                                                                                  acgccGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAatgat                                                                                        <  2:124053/246‑1 (MQ=255)
                                                                                                                                                                     ccGATGGCCTGCATCTCCCAGTGAAACATGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAGATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTGTCGGTACAACAGTGTCAAATTTAATGATGTCTg                                                                                   <  2:535071/251‑1 (MQ=255)
                                                                                                                                                                                        aGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAatgatg                                                                                       <  2:111410/228‑1 (MQ=255)
                                                                                                                                                                                        aGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATATTCAGGTGGGGCCGGCGGTTGTATTTTTCCACGGGCCAACACGGCCTCTATAAACGGAATCAATTTGGGTAATGGTTTTCTACTGCGAAAAATATCACTGCACGCCCAAAGGGCTTCCCGGTCCGGCTCACCAGGGGTACGACGCGCTATCTCAAAAACTTTCGGTACAACAGTGTCAAATTTTatgatgagatcgggagagcacacgccta                                                                 >  1:111410/1‑228 (MQ=255)
                                                                                                                                                                                            aaaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAACACGTCCTCTTTAAGCGGAATCAATTTTGCTAATGTTTTTCTACTGCGAAAAATATCCCCTCACGCCCAACGGAATTCCCCGTCCGGGTCACCCAGGTTTCGACGCGCTATCTCAAAAGCTTTCGGTACAAAAGGGTCAAATTTAATGATGTCTGCAATATCGTAAAAAAAAGaaaa                                                             >  1:38795/1‑249 (MQ=255)
                                                                                                                                                                                                     gCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCt                                                   <  2:403205/251‑1 (MQ=255)
                                                                                                                                                                                                       ggAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTATTTTTTCCGCGGCAAGCACGTCCTCTATAAGCGGAAACAATTTTGATAATGTTTTACTACTGCGAAAAAAATCCCTGCACGCCAAACGGACTTCCCGGTCCGGGCCACCCGGGGTACGACGCGCTATCTCAAAAACTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTgg                                                 >  1:238987/1‑250 (MQ=255)
                                                                                                                                                                                                          aTGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATg                                              <  1:8429/251‑1 (MQ=255)
                                                                                                                                                                                                               aaCAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCACCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATTCCTTCACGCCAAACGGACTTCCCGGTCCGGCTCACCCGGGTTTCGACGCGCTATCTCAAAAAATTTCGGTACCACAGGGTCAAAATTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCCTTCCTGTATTcaca                                          >  1:171350/1‑250 (MQ=255)
                                                                                                                                                                                                                  aGGCTGTGCATCTTCAGGTGGTGCCGGCTGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAGATTCCTGCACGCCAAACGGTCGTCTCGTTCCTGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTGCGGTACAACAGTGTCAAATTTAATTATGTCGGCGATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAAcc                                      <  1:100639/251‑1 (MQ=255)
                                                                                                                                                                                                                                cAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGc                         <  1:149135/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTg              <  1:191540/251‑1 (MQ=255)
                                                                                                                                                                                                                                           ggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGc             >  1:112791/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ggcggTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGATAAGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAACTGGTGCATAACCAGCTGc             >  1:399859/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGGAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGt         <  2:171350/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTAGCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCGGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATc      <  1:248316/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       tctcCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTGCCTTTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgccg  <  1:541598/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
TCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCG  >  NC_000913/2878699‑2879195

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