Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*5914,3840TC67.0% 111.3 / 72.5 91intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base T (18/12);  new base C (40/21);  total (58/33)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGG  >  59/14285‑14481
                                                                                                   |                                                                                                 
cgttccctgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTc                                                                                                  >  1:1102105‑M2/10‑101 (MQ=255)
   tccctgAAGGCCGGGGCAGCCCACATGGATGTGGTCTGGGGGCGCGTTTTAGAGGGATTTTACCTCGCGCCGGCCCCGGTAGCCGTAAGGGATAAGCCGAg                                                                                               <  1:328786‑M2/95‑1 (MQ=255)
    ccctgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGGGGGCGCGTTTTACAGGGATGTTCCCTCGCGCCCGCCCGGGTAGCCGTAAGGGATAAGTCGAgg                                                                                              <  1:658790‑M2/96‑1 (MQ=255)
    ccctgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAgg                                                                                              >  1:1028428‑M2/6‑101 (MQ=255)
      ctgAAGGCGGGGGCAGCCCACATGGATGGGGGCTGAGGGCGCGTTTTACAGGGATGTTCCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGc                                                                                            <  1:1002279‑M2/98‑1 (MQ=255)
      ctgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGc                                                                                            >  2:5995‑M2/4‑101 (MQ=255)
      ctgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGc                                                                                            >  2:942093‑M2/4‑101 (MQ=255)
      ctgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGc                                                                                            >  1:801400‑M2/4‑101 (MQ=255)
       tgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCa                                                                                           >  1:709893‑M2/3‑101 (MQ=255)
       tgAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCa                                                                                           >  2:611857‑M2/3‑101 (MQ=255)
         aaGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCAcc                                                                                         >  2:772369/1‑101 (MQ=255)
           ggCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACcgc                                                                                       >  2:39763/1‑101 (MQ=255)
            gcgggggCAGCCCACAGGGTTGTGGGCTGAGGGCGCTTTTTACAGGGATGTTCCCCCGCCCGCGACCCGCTAGCCGTAAGGGATAAGTCGAGGGCACcgcg                                                                                      <  2:1323402/98‑1 (MQ=255)
              cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCa                                                                                    >  1:1599204/1‑101 (MQ=255)
               ggggCAGCCCACATGTGTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTCCCTCGCCCCCGCCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAg                                                                                   <  1:1145246/101‑1 (MQ=255)
               ggggCACCCCACATGGATGTGGGCTGAGGGCGCGTTTCACGGGGATGTCACCTCGCCCCCGACCCGTTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAg                                                                                   <  1:733337/101‑1 (MQ=255)
                 ggCACCCCACATGGATGTGGGTTGAGGGCGCTTTTTAGGGGGATGTTCCCTCGCGCCGGCCCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAgcg                                                                                 <  2:202476/101‑1 (MQ=255)
                     gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGa                                                                             >  1:1669814/1‑101 (MQ=255)
                     gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGa                                                                             >  1:650422/1‑101 (MQ=255)
                     gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGa                                                                             >  1:966733/1‑101 (MQ=255)
                      cccctgtGGTTGGGGGATGAGGGCGCGTTTAACGGGGTTGTTCCCCGGCGCCCGACCGGTTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGAt                                                                            <  2:326238/95‑1 (MQ=255)
                        cacaGGGATGTGGGCTGGGGGCGCGTTTTACAGGGATGTTCCCTCGCGCCCGCCCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGAttt                                                                          <  2:1530274/101‑1 (MQ=255)
                         gcaTGGATGGGCGGCGAGGGCGCATTTTACAGGGCCGTTCCCTTGCCCCCGGCCGGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGAtttt                                                                         <  2:24711/100‑1 (MQ=255)
                         agagggttgtgGGCGGGGGGCGCTTTTTACGGGGATGTTCCCCCGCGCCCGCCCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGAtttt                                                                         <  2:801395/94‑1 (MQ=255)
                         acagggtggggggCTGGGGGCGTTTTTTACAGGGTTGTTCCCCCGCGCCCGCCCGGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGAtttt                                                                         <  2:1491742/91‑1 (MQ=255)
                         acaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGAtttt                                                                         >  1:1071898/1‑101 (MQ=255)
                           aTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGt                                                                       >  2:629832/1‑101 (MQ=255)
                            tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTCCCTCCCGCCCGCCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGtt                                                                      <  1:749615/101‑1 (MQ=255)
                             ggtggggggCGGAGGGCGCGTTTTACGGGGATGTTCCCTCGCCCCCGCCCCGTTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTc                                                                     <  2:1562694/95‑1 (MQ=255)
                              gatgcgcggcgAGGGCGCATCTTACAGGGACGTTCCTTTGCGCCCGGCCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCg                                                                    <  2:1629697/91‑1 (MQ=12)
                              gATGTGGGCTGGGGGCGCGTTTTACAGGGATGTTCCCTCGCCCCCGCCCGGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCg                                                                    <  2:1270496/101‑1 (MQ=255)
                              gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCg                                                                    >  1:12451/1‑101 (MQ=255)
                               aTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGc                                                                   >  2:478807/1‑101 (MQ=255)
                                tgtgGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGcc                                                                  >  1:1648765/1‑101 (MQ=255)
                                  ggggCGGAGGGCGCGTTTTAGAGGGATGTTCCCTCGCCCCCGCCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCag                                                                <  1:810963/100‑1 (MQ=255)
                                  cggggggAGGGCCCATCTTACAGGGACGTTCCCTGGCCCCCGACCGGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCag                                                                <  2:1433442/95‑1 (MQ=255)
                                     gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGc                                                             >  1:1086617/1‑101 (MQ=255)
                                     gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGc                                                             >  1:738273/1‑101 (MQ=255)
                                       tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGccc                                                           >  2:1577895/1‑101 (MQ=255)
                                          gggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCggg                                                        >  1:934905/1‑101 (MQ=255)
                                           ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCgggg                                                       >  1:1382860/1‑101 (MQ=255)
                                            gcgttttttACGGGGATTTTCCCCCGCGCCCGCCCGGTTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCggggg                                                      <  2:736152/96‑1 (MQ=255)
                                            gcgttttttACAGGGATGTTCCCCCGCGCCCGCCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCggggg                                                      <  2:970586/96‑1 (MQ=255)
                                            gcgcGTTTTACAGGGATGTTCCCTCGCGCCCGACCCGTTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCggggg                                                      <  1:420565/101‑1 (MQ=255)
                                            gcgcGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCggggg                                                      >  1:137079/1‑101 (MQ=255)
                                             cgcgTTTTACAGGGATGTTCCCTCGCGCCCGCCCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGt                                                     <  2:405038/101‑1 (MQ=255)
                                              gcttcttACAGGGACGTTCCCTTGCGCCCGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTg                                                    <  1:1476339/96‑1 (MQ=255)
                                              gcatcttACAGGGACGTTCCCTTGCGCCCGCCCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTg                                                    <  2:452122/96‑1 (MQ=25)
                                                atcttACAGGGACGTCCCCTTGCGCCCGCCCGGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCa                                                  <  2:232986/98‑1 (MQ=255)
                                                   ttagggggATGTTCCCCCGCCCCCGCCCGGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggg                                               <  2:133748/96‑1 (MQ=255)
                                                   ttACAGGGTTGTTCCCCCCCCCCCGCCCCGGTAGCCGTAAGGGATAAGACGCGGGGACCGCCCAGCGGCGATTTTGTTCGCTAGACGCCGGGGGTGCAggg                                               <  2:515581/101‑1 (MQ=255)
                                                    tACAGGGTTGTCCCCCCCCCCCCGCCCTGTTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAgggg                                              <  2:1205119/101‑1 (MQ=255)
                                                    tACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAgggg                                              >  2:706908/1‑101 (MQ=255)
                                                    tACAGGGACGTTCCCTTGCGCCCGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAgggg                                              >  2:150983/1‑101 (MQ=255)
                                                     aCAGGGATGTTCCCGCGCGCCCGCCCCGGTAGCTGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                                             <  2:141768/101‑1 (MQ=255)
                                                     aCAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                                             >  1:304573/1‑101 (MQ=255)
                                                     aCAGGGACGTTCCCTTGCGCCCGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCACAGCCTGGGGGTGCAggggg                                             >  1:473100/1‑101 (MQ=255)
                                                          gacgTTCCCTTGCCCCCGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGggcggcg                                        <  2:214565/98‑1 (MQ=255)
                                                            tGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcggc                                      >  2:167413/1‑101 (MQ=255)
                                                             tttcccccgcgCCCGACCCGTTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGgcggcggcg                                     <  2:1423976/94‑1 (MQ=255)
                                                             gTTCCCTTGCGCCCGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGgcggcggcg                                     >  1:1295600/1‑101 (MQ=255)
                                                              ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGa                                    >  2:1399518/1‑101 (MQ=255)
                                                               tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt                                   >  2:657148/1‑101 (MQ=255)
                                                               tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt                                   >  2:338899/1‑101 (MQ=255)
                                                                cccTTGCGCCCGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGAtt                                  >  2:122855/2‑101 (MQ=255)
                                                                 ccccgcgCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTg                                 <  2:1701162/98‑1 (MQ=255)
                                                                 cattgcgcCGGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTg                                 <  2:1181320/97‑1 (MQ=255)
                                                                  cTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                >  2:294958/1‑101 (MQ=255)
                                                                    tgcgcCCGCCCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGgcc                              <  2:879598/100‑1 (MQ=255)
                                                                     gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccg                             >  2:765813/1‑101 (MQ=255)
                                                                      cgcCCGACCCGGTAGCCGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccgc                            >  2:546290/1‑101 (MQ=255)
                                                                       gcCCGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGgccgcc                           >  2:1620909/1‑101 (MQ=38)
                                                                        cccGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                          >  2:789891/1‑101 (MQ=255)
                                                                        cccGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGGCCGccc                          >  1:1147036/1‑101 (MQ=38)
                                                                         ccGACCCGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGGCCGcccc                         >  2:1712877/1‑101 (MQ=38)
                                                                              ccGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                    >  1:286997/1‑101 (MQ=37)
                                                                              ccGGTAGCCGCAAGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcg                    >  2:47935/1‑101 (MQ=37)
                                                                                 gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgcgc                 >  1:367712/1‑101 (MQ=255)
                                                                                 gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTgcgcgc                 >  2:134319/1‑101 (MQ=255)
                                                                                     ccGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCt             >  2:728538/1‑101 (MQ=255)
                                                                                      cgcaaGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGTTTGGCCGCCCCCTGCGCGCTccct            >  1:1560693/4‑99 (MQ=21)
                                                                                      cGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCtt            >  2:1035790/1‑101 (MQ=255)
                                                                                      cGTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCtt            >  2:1664667/1‑101 (MQ=255)
                                                                                       gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg           >  1:1363511/1‑101 (MQ=255)
                                                                                       gTAAGGGATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTg           >  2:567487/1‑101 (MQ=255)
                                                                                         aaGGGATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCCTgcg         >  2:333439/1‑101 (MQ=25)
                                                                                            ggATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCa      >  1:1485215/1‑101 (MQ=255)
                                                                                            ggATAAGCTGAGGGAATCGCGCAGCGGCGATTTTGTTCGCCAGAGCCTGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCCTGCGCCa      >  2:1212643/1‑101 (MQ=21)
                                                                                             gATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg     >  2:1343699/1‑101 (MQ=255)
                                                                                             gATAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAg     >  2:1093422/1‑101 (MQ=255)
                                                                                               tAAGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTg   >  2:1308234/1‑101 (MQ=255)
                                                                                                aaGCCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTga  >  1:927150/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGG  >  59/14285‑14481

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