Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*21180TC48.1% 67.1 / 19.4 54intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base T (1/27);  new base C (26/0);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 2.88e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

‑‑‑‑‑‑‑‑‑‑‑‑GATGCAAGGCATTAATAGATAAAGCCAGATCAACTGCCGACGAAATCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  2/1‑196
                                                                                                                                   |                                                                                                                                                                                                                                                      
aaaaaagaaaaggatgcaaggaattAATAGAAAAAGCCAGATCAACTGCCGAAGAAATCATAAGAGACCAAAAAAATTAAAAA‑GATCCTATTTTGTCCGTTGATACA‑GAAAAAGCTTTTCTAATAATTTCCCCACAGTGGAAACCGCAGCTACTGGAAAGTTTGCCGCCCCGAAAAGGGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaacc                                                                                                                               >  1:770469‑M2/13‑208 (MQ=255)
aaaaaagaaaaggatgaaaggaattAATAGATAAAGCCAGATCAACTGCCGACGAAATCATAAGAGACCAAAAAAATTAAAAA‑GATACTATTTGGCCCGTTGATACA‑GACAAAGCTTTTTTAATAATTTCCCCACAGTGGGAACCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGAACTCCCTGCTACCCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaacc                                                                                                                               >  1:928384‑M2/13‑208 (MQ=255)
            gaggaaaggaatTAATAAAAAAACCCAGATCATCTGCCGACGAAATCATTAGAGACATAAAAAATTTAAAA‑AATACTATTTTGTCCGTTGATACA‑GACAAAGCTTTTCTAATAATTTCCCCACAGTGGAAACCGCAGCTACTGAAAAGTTTGCCCCCCGGTAAAGCGGGACTACCTCCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagcc                                                                                                                    >  1:389194‑M2/6‑196 (MQ=255)
            gATGCAAGGCATTAATAGATAAAGCCAGATCAACTGCCGACGAAATCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagac                                                                                                                                                                <  1:919177‑M2/206‑11 (MQ=255)
                   ggCATTAATAGATAAAGCCAGATCAACTGCCGACGAAATCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagc                                                                                                                     <  1:3065‑M2/242‑54 (MQ=255)
                     aataaataaaaaaagccagatcaacTGCCGCCGAAATCATTAGAGCCATAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GAAAAAGCGTTTCTAATAATTTCCCCACAGTGGGAACCGCAGCTACTGAAAAGTTTGCCGCCCCGAAAAGCGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtg                                                                                                           >  1:388126‑M2/12‑187 (MQ=255)
                      aTTAATAAAAAAAGCCAGATCACCTCCCGACAAAATAATTAGAGACAAAAAAAATTAAAAA‑GATATTATTTTGTCCGTTGATAAA‑GAAAAAGCGTTTCTAATAATTTCCCCACAGGGGGAACCGCAGCTACTGGAAAGTTTGCCCCCCCGAAAAGCGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgag                                                                                                         >  1:622878‑M2/1‑186 (MQ=255)
                      aTTAATAAAAAAAGCCAGATCAACTGCCGACAAAATCATTAGAGCCCAAAAAAATTAAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCTTTTCTAATAATTTCCCCACAGTGGGAATCGCAGCTACTGAAAATTTTGCCGCCCCGTAAAGCGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgag                                                                                                         >  1:545698‑M2/1‑186 (MQ=255)
                       taaataGATAAAGCCAGATCAACTGCCGACAAAATCATTAGAGACCAAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTCCCCACAGGGGAAACCGCAGCTACTGGAAAGTTTGCCCCCCCGTAAAGCGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagg                                                                                                        >  1:494321‑M2/3‑185 (MQ=255)
                                 aaGCCAGATCAACTGCCGACAAAATCATTAGAGACATAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCTTTTCTAATATTTTCCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGAAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttat                                                                                              >  1:299609‑M2/1‑175 (MQ=255)
                                 aaGCCAAATCAACTGCCGACAAAATCTTTAGAGACCTAAAAAATTAAAAA‑GATTCTATTTTGTCCGTTGATACA‑GACAAGGCGTTTCTAATAATTTCCCCACAGGGGGAACCGCAGCTACTGGAAAGTTGGCCCCCCGGTAAGGCGGACCTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttat                                                                                              >  1:239617‑M2/1‑175 (MQ=255)
                                  aGCCAGATCACCTGCCGACGAAATCATAAGAGACCAAAAAAATTAAAAA‑GATCCTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTCCCCACAGGGGGAATCGCAGCTACGGAAAAGTTTGCCGCCCCGTAAAGGGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaacca                                                                                                                              >  1:1039353‑M2/1‑174 (MQ=255)
                                   cccaGATCACCTGCCGACAAAATCATAAGAGACA‑AAAAAATTTAAAAAGATACTATTTTGTCGGTTGATACA‑GACAAAGCTTTTCAAATATTTTCCCCACAGTGGAAACCGCAGCTACTGAAAAGTTTGCCCCCCCGTAAAGCGGGACTACCTGCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttat                                                                                              >  1:520745‑M2/2‑173 (MQ=255)
                                                       aaTCATAAGAGACCAAAAAAATTTAAAA‑GATCCTATTTTGTCCGTTGATACA‑GACAAAGCTTTTCTAATATTTTCCCCACAGTGGGAACCGCAGCTACTGGAAAGTTTGCCCCCCCGAAAAGCGGGACTCCCTGCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttat                                                                                              >  1:502875‑M2/1‑153 (MQ=255)
                                                        ataattaGAGACCAAAAAAATTAAAAA‑GATCCTATTTTGTCCGTTGATACA‑GAAAAGGCGTTTCTAATATTTTCCCCACAGGGGAAACCGCAGCTACGGAAAAGTTTGCCCCCCCGTAAAGCGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagc                                                                                                                     >  1:420817‑M2/4‑152 (MQ=255)
                                                        aTCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCAAATAATTTCCCCACAGTGGGAACCGCAGCTACTGGAAATTTTGCCCCCCCGAAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccc                                                                                                                          >  1:417251‑M2/1‑152 (MQ=255)
                                                         tCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatg                                                                      <  1:643880‑M2/251‑101 (MQ=255)
                                                         tCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaa                                                                                    <  1:170076‑M2/237‑87 (MQ=255)
                                                         tCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttat                                                                                              <  1:790988‑M2/227‑77 (MQ=255)
                                                          cATTAGAGACCTAAAAATTTTAAAA‑GATTTTATTTTGTCCGTTGATACA‑GACAAAGTGTTTCTAATAATTTCCCCACAGTGGGAATCGCAGCTTCTGGAAATTTTGCCGCCCCGTAAAGGGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgt                                                                     >  1:423133‑M2/1‑150 (MQ=255)
                                                          cATTAGAGACCAAAAAAATTAAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATATTTTCCCCACAGGGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGGGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgt                                                                     >  1:79453‑M2/1‑150 (MQ=255)
                                                          aataagagaccaAAAAAATTAAAAA‑GATCTTATTTTGCCCGTTAA‑ACAAGACAAAGCGTTTCTAATAATTTCCCCACAGGGGGAACCGCAGCTACTGAAAAGTTTGCCGCCCCGTAAAGGGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgt                                                                     >  1:303126‑M2/5‑150 (MQ=255)
                                                          aatTAGAGACCAAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑AAAAAAGCTTTCCTAATAATTTCCCCACAGGGGGAACCGCAGCTACTGAAAAGTTTGCCCCCCCGAAAAGCGGGACTCCCTGCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgt                                                                     >  1:257110‑M2/2‑150 (MQ=255)
                                                               gaagatgtgtaaaagaaaca‑gatcctattttgtccgttaATACA‑GAAAAAGCGTTTCTAATAATTTCCCCACAGTGGGAATCGCAGCTACTGAAAAGTTTCCCCCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaa                                                                 >  1:897697‑M2/20‑145 (MQ=255)
                                                                 aacCTAAAAAATTAAAAA‑GATCCTATTTTGTCCTTTGATACA‑GAAAAAGCTTTTCTAATAATTTCCCCCCAGGGGGAACCGCAGCTCCTGAAAAGTTTGCCGCCCCGTAAAGCGGAACTCCCTGCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatga                                                              >  1:969394‑M2/2‑143 (MQ=255)
                                                                        aaaaTTAAAAA‑AATACTATTTTGTCCGTTGATACA‑GACAAAGCTTTTCAAATAATTTCCCCACAGGGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGGGGGACTCCCTGCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaa                                                       >  1:547878‑M2/1‑136 (MQ=255)
                                                                           aTTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaat                                                                       <  1:451243‑M2/232‑100 (MQ=255)
                                                                                     atatTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATATTTTCCCCACAGTGGGATCCGCAGCAACTGAAAAGTTTCCCGCCCCGAAAAGCGGGACTCCCTCCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatg                                                                      >  1:903742‑M2/1‑124 (MQ=255)
                                                                                     aTACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttat                                                                                              <  1:85818‑M2/200‑77 (MQ=255)
                                                                                     aTACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcc                                                                                                     <  1:546260‑M2/193‑70 (MQ=255)
                                                                                     aTACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtg                                                                                                           <  1:370111‑M2/187‑64 (MQ=255)
                                                                                     aTACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtg                                                                                                           <  1:279680‑M2/187‑64 (MQ=255)
                                                                                        cTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgttttta                                        <  1:412338‑M2/251‑131 (MQ=255)
                                                                                          tttttgTCCGTTGATACA‑GAAAAAGCGTTTCTAATATTTTCCCCACAGGGGGAATCGCAGCTACTGAAAAGTTTCCCGCCCCGTAAAGCGGGACTCCCTGCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaacc                                                                                                                                     >  1:902128‑M2/2‑119 (MQ=255)
                                                                                          aTTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaa                                      <  1:753411‑M2/251‑133 (MQ=255)
                                                                                          aTTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtg                                                                                                           <  1:502030‑M2/182‑64 (MQ=255)
                                                                                          aTTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtg                                                                                                           <  1:793834‑M2/182‑64 (MQ=255)
                                                                                          aTTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccc                                                                                                                                    <  1:810617‑M2/157‑39 (MQ=255)
                                                                                          aTTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtat                                                                                                                                                                      <  1:613462‑M2/123‑5 (MQ=255)
                                                                                                 ccGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaat                               <  1:705001‑M2/251‑140 (MQ=255)
                                                                                                 ccGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaat                               <  1:196681‑M2/251‑140 (MQ=255)
                                                                                                  cGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaatt                              <  1:1064724‑M2/251‑141 (MQ=255)
                                                                                                  cGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaatt                              <  1:511000‑M2/251‑141 (MQ=255)
                                                                                                  cGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacg                                                                                                             <  1:479763‑M2/172‑62 (MQ=255)
                                                                                                      gATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacgg                          <  1:618366‑M2/251‑145 (MQ=255)
                                                                                                      gATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacgg                          <  1:847246‑M2/251‑145 (MQ=255)
                                                                                                      gATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatct                                                                                                                                                                    <  1:713584‑M2/113‑7 (MQ=255)
                                                                                                             gACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacggtatga                     <  1:477176‑M2/250‑150 (MQ=255)
                                                                                                             aaaaaagctTTTCTAATATTTTCCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAGGCGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaagg                                                                                 >  1:706067‑M2/4‑101 (MQ=255)
                                                                                                                 aaGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattac                            <  1:440602‑M2/239‑143 (MQ=255)
                                                                                                                    cGTTTCAAATAATTTCCCCACAGGGGGAATCCCAGCTACTGGAAAGTTTGCCCCCCCGAAAAGCGGGACTCCCTGCTACCCCCGCCGATGTATCgtatctagacaaatacaaggcaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacggtatgacgcaatcg             >  1:369850‑M2/1‑94 (MQ=255)
                                                                                                                        tttaataattTCCCCCCAGGGGGAACCCCACCTACTGAAAATTTTGCCCCCCCGAAAAGCGGACCTACCTCCTCCCCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaacccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacggtatgacgcaatcggggga         >  1:598660‑M2/3‑90 (MQ=255)
                                                                                                                         cTAATAATTTTCCCACAGGGGGAACCGCAGCTACGGGAAAGTTGCCCGCCCCGAAAAGGGGGACTCCCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaat                                                                       >  1:209824‑M2/1‑89 (MQ=255)
                                                                                                                               ttttccCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCCCCCCGTAAAGGGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacggtatgacgcaatcgggggataagtg  >  1:1015369‑M2/2‑83 (MQ=255)
                                                                                                                               ttttccCCACAGTGGGAACCGCAGCTACTGGAAAGTTTGCCCCCCCGTAAAGCGGGACTACCTGCTCCTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacggtatgacgcaatcgggggataagtg  >  1:604657‑M2/2‑83 (MQ=255)
                                                                                                                               atTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagacaaatacaagccaaagaaataaacaaccctaaccaaccccaagccctccacgtgagggcctttttattactcaaaaaaggacaataaaatgttgaatgaacaagaaaaacacgcgtttttaaatggcaattacggtatgacgcaatcgggggataagtg  <  1:474476‑M2/251‑169 (MQ=255)
                                                                                                                                  ttCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATCgtatctagac                                                                                                                                                                <  1:329015‑M2/90‑11 (MQ=255)
                                                                                                                                   |                                                                                                                                                                                                                                                      
‑‑‑‑‑‑‑‑‑‑‑‑GATGCAAGGCATTAATAGATAAAGCCAGATCAACTGCCGACGAAATCATTAGAGACCTAAAAAATTTAAAA‑GATACTATTTTGTCCGTTGATACA‑GACAAAGCGTTTCTAATAATTTTCCCACAGTGGGAATCGCAGCTACTGGAAAGTTTGCCGCCCCGTAAAGCGGGACTACCTGCTACTCCCGCCGATGTATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  2/1‑196

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