Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,9470CT25.6% 72.8 / 10.1 44E276K (GAA→AAA) cas1multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (7/25);  new base T (11/0);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 5.53e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATA  >  NC_000913/2878704‑2879175
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
cAGTCCAGCTATTTGGTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAGCAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTGCTCCAGCg                                                                                                                                                                                                                                   <  2:47645/251‑1 (MQ=255)
            ttGTTCCCAGATCATTTCACGAATTGTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTGCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCAGCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc                                                                                                                                                                                                                        <  2:106737/250‑1 (MQ=255)
                tCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCctct                                                                                                                                                                                                                    <  1:241964/250‑1 (MQ=255)
                           ttCACGAATTTTCAAGGACACACCACCTACATATACCACTGCACGCACCTCCACCAACCAGATGGCTACTCTGCCTCGTAAGAGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGCCCGGCGGTTTTTTTTTTCCAGCGGCCAGCACgctcatactaaacgcact                                                                                                                                                                                                         >  2:626799/1‑233 (MQ=255)
                              aCGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAAt                                                                                                                                                                                                   <  1:332794/251‑1 (MQ=255)
                                    ttttGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACACCCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGAACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGCAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCCGGTGGGGCCGGCGGTTGGGTTTTTCCCGAGGCCAGAACGCCCCCTATTAGCGGGATCAAttttggc                                                                                                                                                                                             >  1:355937/1‑247 (MQ=255)
                                          ggATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGCTGTTTTTTTCCTGCCGCCAGCACCTCCTCtatacggcgagtctatatggttactgtt                                                                                                                                                                                         >  2:334258/1‑227 (MQ=255)
                                           gATACATCACCTACATCTACCCCTGCACGTACCTCCAACAACCAGATGGCGAATCTGCCTCGTAGGCGCGGAGGTACATTTTCAGTGACCCCGACCAACATACTCATTTCAGCTACTCCGATGGCCGGCAGCTCCCAGTGAAACAGGAAGCGGAAGGGCAACAGTCTGTGCCTCCGCAGGTGGGTCCGGCGGTGGTATTTCTCCAGCGTCGAGCTCTTCCTCGATAAGCGGAATCAATTGGGCTAATGttt                                                                                                                                                                                      <  2:415398/251‑1 (MQ=255)
                                            aTACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGtttt                                                                                                                                                                                     <  1:513945/251‑1 (MQ=255)
                                              aCATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTtac                                                                                                                                                                                   <  2:135980/251‑1 (MQ=255)
                                               cATCACCTACATATCCCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTACGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCACCGGCAAGCACCTCCTCTCTAAGACGAAACAATATGTCTAATGTTTtact                                                                                                                                                                                  >  2:307656/1‑251 (MQ=255)
                                                   aCCTACAGATTCCCCTGCACGTACCCCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACAGACTCATGTCAGCTCCTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGTCCGGCTGTTGTATTTCTCCAGCGGCAATCCTGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTATACTACTg                                                                                                                                                                              <  1:425240/251‑1 (MQ=255)
                                                       aCATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGTTTGTATTTCTCCAGCGGc                                                                                                                                                                                                                                 <  2:429029/198‑1 (MQ=255)
                                                           ataCCCCTGCGCGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCGCTATAAGCGGAATCAATTTGGCTAATGTTTCACTACTGCTAAAAat                                                                                                                                                                      <  2:105693/251‑1 (MQ=255)
                                                                ccTGCACGTCCCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATccc                                                                                                                                                                 <  1:3598/251‑1 (MQ=255)
                                                                            tCCAACAACCAGCTGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTTCTACTTCGAAAAATATACCTGCACGCCCAAc                                                                                                                                                    >  1:515404/1‑251 (MQ=255)
                                                                            cccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACGCAGTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCTGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAAc                                                                                                                                                    <  2:299453/250‑1 (MQ=255)
                                                                               aacaacCCGATGGCTAATCTGCCTCGTAAGCTCGGAGGTACATTTTCAGTGACCACGACCAACATACTCACTGCAGCGACGCCTATTGCCTGCATCTCCCAGTTAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGTTGCCAGCGGTTGTATATCTCCAGCCTCGAGCACGTCCGCTAGAAACGGAATCAATTTGGCTAATGTTGTACCACTGTGAATAAGACCCCTGCGCTCCAAACGGa                                                                                                                                                <  2:428822/251‑1 (MQ=255)
                                                                               aacaacCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCCCGTCCTCTATTAACGGAAACAATTTGGCTAATGTTTTTCTTCTTCGAAAAAAATCCCTTCACCCCCAACggc                                                                                                                                                 >  1:290384/1‑250 (MQ=255)
                                                                                  aacCAGCTGGATAATCTGCCTCGTAAGCGCGGCGGTACATTTTCAGTGACCACGACCAACCTACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTTCCAACGGGAAGCACGTCCTCTATTAGAGGTGTTCAATTGGCTATTGTTTTTCTATAGCGAAAAataatcccgcgcgccacacagacat                                                                                                                                               >  2:295243/1‑229 (MQ=255)
                                                                                                          aaGCACGGAGGTACATTTTCAGTGACCACCACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTCTGCATCTTCAGGTGGGGCCGGAGGTTGTATTTCTCTACCGGCAAGCAAGTCCCCTATTAAAGGCATCAAATTTTCTAATTTTTTTTTACTTCGACAAATATTCATCCACGCCCAACGGGCCTTTCCGTCCGGCTCACCCGGGttct                                                                                                                     >  1:293345/1‑249 (MQ=255)
                                                                                                           aGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTACCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTTCCAGCCGCAAGCACCTCCTCTCTAAACGGAAACAATTTGGCTAAAGTTTTACTACTGCGAAAAAAATCCCTGCACGCCCAACGGACTTCCCCGTCCGGCTCACCCGGGTTacg                                                                                                                    >  1:352466/1‑251 (MQ=255)
                                                                                                             cgcgGAGGTACGTTTTCAGTGACCACGACCAACAGACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTGCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCGCTTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAGATCCCTGCACGCCAAACGGACTTCCGGGTCCGGCTCACCAGGGTTacgac                                                                                                                  <  2:541423/251‑1 (MQ=255)
                                                                                                                     tACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTTTTTTTCCCGCGGCAAGCACGTCCTCCATAAGCGGAATCAATTTTGCTAATGTTTTTCTACTGCGAAAAATATCCCTGCACGCCAAACGGACGTCCCCGTCCCGCTCACCAGGGGTACGAcgccccatc                                                                                                          >  1:561606/1‑245 (MQ=255)
                                                                                                                       cATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTTCAGCGGCCAGCAAGTCCTCTATAAGCGGGATCAATTTTGCTTATGGTTTTCTACTTCGAGAAATATCCCTGCAACCCCAACGGACTTCCCGGTTCGGCTCCCCAGGGGTacgaaccgccatccc                                                                                                        >  2:467195/1‑240 (MQ=255)
                                                                                                                               gTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCACAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCCCGTTCTCTATAAACGGGATCAAATTTGGTAATGTTTTACTACTGCGAAAAATATCCCTTCACGCCAAACCGACTTCCCGGGCCCGCTCCCCACGGGTACCGCGAGCTATCTCAAAAGCtt                                                                                                >  1:88416/1‑251 (MQ=255)
                                                                                                                                           aaCATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCCTCTTCCGGTGGGGCCGGCGGTTGTTATTTTCCCGCGGCAAGCACGGCCTCTATCAGCGGGAACAATTTGGGTCATGTTTTTCTACTTCGCAAAATATCCCCGCACGCCAAACGGACTTCCCGGCCCCGCTCACCCCGGGTACGACCCGCCATTTCCAAAGCtgtctgggcaac                                                                                      >  2:465408/1‑238 (MQ=255)
                                                                                                                                           aaCATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAg                                                                                    <  1:275928/251‑1 (MQ=255)
                                                                                                                                                            tACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCCCCAACGGCAAGCACGTCCTCTATAAGGGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAAAAACCCTGCACCCCAAACGGACTTCCCGGGCCGGCTCACCAGGGGTACCACGCGCTATCTCAAAAACTTTTGGTACAAACGGGTCAAAttttaagatt                                                                   >  1:553212/1‑244 (MQ=255)
                                                                                                                                                                ccGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTg                                                               <  1:123980/251‑1 (MQ=255)
                                                                                                                                                                 cGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGc                                                              <  1:468844/251‑1 (MQ=255)
                                                                                                                                                                       ccagCATCTACCATTGAAACAGGAAGCGGAATTGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAATCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTTCTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCTGTCCGGCTCACCAGGGTTACGACGCTCTATCTCACAATCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTTCAATATc                                                        <  2:70974/248‑1 (MQ=255)
                                                                                                                                                                       ccTGCATCTCCCAGGGAAACAGGAATCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGTGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTGTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATc                                                        <  2:202790/251‑1 (MQ=255)
                                                                                                                                                                            aTCTCCCAGTGAAACAGGCAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGCTTTATTTTTCCCACGGCAAGCAAGTCCTCTATAAACCGAATCAATTTTGCTAATGTTTTTCTCCTGCGACAAAAATCCCTGCACGCCAAACAGGCTTCCCTGTCCGGATCACCAGGGTTACGACGCGCTCTTTCCAAAGCTTTTTGTACAACAATGTCAAAATTTatgatttttgccataaccata                                                    >  2:578713/1‑234 (MQ=255)
                                                                                                                                                                                cccAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCCACGGCCAACACCTCCTCTATAAACGGAATCAATTTTGGTTATGGTTTTCTACTGCCACAAATATCCCTTCACGCCAAACGGACTTTCCCGTCCGGCCCACCCCGGGTACGAAGCGCTATCTCAAAAACTTTCGGTGCCACAGTGTCAAATTTTATTATTTCTGCAataatgtaaaaaa                                                >  2:334912/1‑240 (MQ=255)
                                                                                                                                                                                    gTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAgt                                                                                   <  1:128163/211‑1 (MQ=255)
                                                                                                                                                                                        aaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGGGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATAGCGTAAACAAAGGAAagag                                       <  1:172685/251‑1 (MQ=255)
                                                                                                                                                                                                   gAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGGCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAATCTTTCGGTACAACAGTGTCAAAt                                                                            <  2:428383/203‑1 (MQ=255)
                                                                                                                                                                                                   gAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTTTTTCTCCAGCGGGAAGCACGTCCCCTCTAAGAGGAATATATTTTGGTTATATTTTACTCCTGCGAAAAATATCCATGACAGCCAAACAGACTTCCCGGTCCCGCTCACCAGGGTTACGACGCGCTATCTCCAAAAATTTTGGTACAACAATGTCaaaa                                                                            >  1:428383/1‑202 (MQ=255)
                                                                                                                                                                                                                cTGTGCATCTTCAGGTGGGGCCGGCGGTTGGTTTTCCTCCACGGGCAGCCACTCCTCCATAAGCGGGAACAATTTTGATAATGGTTTACTACTGCGAAAAAAAACCCTGCACGCCAAACGGACTTCCCGGGCCCGGTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTTGGTACCACAGTGTCAAATTTAATGATGTCTGCACAATCGCAAACAAAGGGAAAAGGCTTTTCTGTATGCACAAACCCa                >  2:233910/1‑250 (MQ=255)
                                                                                                                                                                                                                 tgtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTCCAGCGGCAAGCACGTCCTCTATAAGCGGGATCAAATTGGCTCATGTTTTTCTACTCCGGAAAAAATCCCCGCACGCCAAACGGACTTTCCGGGCCCGCCCACCACGGGTACGACGCGCTATCTCAAAAACTTTCGGTACAACAGAGGCAAATTTAAAGATGTTTGCAATATCGTAAAAAAAAGAAAAAGGGTTTCCTGTATGCACAAAACCAAt              >  1:405182/1‑251 (MQ=255)
                                                                                                                                                                                                                          tCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTg      <  2:515404/250‑1 (MQ=255)
                                                                                                                                                                                                                          tCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACGACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTg      <  2:290384/250‑1 (MQ=255)
                                                                                                                                                                                                                             ggTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATa  <  1:490837/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
CAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATA  >  NC_000913/2878704‑2879175

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: