Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*13312,6980CG63.5% 60.1 / 190.9 204intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base G (44/57);  total (90/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 208 total aligned reads displayed.

GTCATCATGTAGTGGAAGTGCCTGGAGGAGGTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATT  >  133/12604‑12798
                                                                                              |                                                                                                    
        gTAGTGGAAGTGCCTGGAGGAGGTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCttt                                                                                        >  2:202925/1‑101 (MQ=255)
        gTAGTGGAAGTGCCTGGAGGAGGTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCttt                                                                                        <  2:784266/101‑1 (MQ=255)
                     cTGGAGGAGGTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGt                                                                           >  1:140575/1‑101 (MQ=255)
                         aggaggTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTa                                                                       <  2:1155505/101‑1 (MQ=255)
                            aggTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAAcc                                                                    >  2:296910/1‑101 (MQ=255)
                                   aaagcaggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCtt                                                             >  1:227341/4‑101 (MQ=11)
                                      gcaggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGaa                                                          <  2:766571/101‑1 (MQ=11)
                                       caggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAAt                                                         <  2:885817/101‑1 (MQ=11)
                                       caggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAAt                                                         <  2:1183399/101‑1 (MQ=11)
                                        aggGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATa                                                        >  1:1033775/1‑101 (MQ=255)
                                        aggGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATa                                                        >  2:1078537/1‑101 (MQ=255)
                                        aggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATa                                                        <  2:878644/101‑1 (MQ=11)
                                         gggtattactACGAAAGCCCGCTCCCCGCGAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATAc                                                       <  1:583040/92‑1 (MQ=11)
                                         gggCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATAc                                                       <  1:325971/101‑1 (MQ=255)
                                          ggCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATAcc                                                      >  2:908685/1‑101 (MQ=255)
                                          ggCAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATAcc                                                      >  1:288691/1‑101 (MQ=255)
                                           gCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCt                                                     >  2:776717/1‑101 (MQ=255)
                                           gCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCt                                                     >  2:278387/1‑101 (MQ=255)
                                             aaGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTa                                                   >  2:1321788/1‑101 (MQ=14)
                                              aGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTAt                                                  <  1:947310/101‑1 (MQ=255)
                                              aGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTAt                                                  <  2:1706495/101‑1 (MQ=255)
                                              aGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTAt                                                  <  2:1348768/101‑1 (MQ=14)
                                               gATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATc                                                 <  1:1236896/101‑1 (MQ=255)
                                               gATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATc                                                 <  1:1703886/101‑1 (MQ=255)
                                                actACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCa                                                >  1:760160/3‑101 (MQ=14)
                                                aTTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCa                                                >  2:658400/1‑101 (MQ=14)
                                                aTTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCa                                                >  2:436573/1‑101 (MQ=14)
                                                aTTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCa                                                <  2:112080/101‑1 (MQ=255)
                                                 ttACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCAt                                               <  1:1665516/101‑1 (MQ=255)
                                                 ttACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCAt                                               <  2:902329/101‑1 (MQ=17)
                                                 ttACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCAt                                               <  1:1161970/101‑1 (MQ=18)
                                                  tACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCAtt                                              >  2:515816/1‑101 (MQ=255)
                                                  tACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCAtt                                              <  1:990288/101‑1 (MQ=14)
                                                   aCGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTc                                             <  1:1258838/101‑1 (MQ=255)
                                                   aCGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTc                                             <  2:582519/101‑1 (MQ=255)
                                                    cGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCg                                            >  1:1558206/1‑101 (MQ=17)
                                                     gAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGt                                           >  1:892097/1‑101 (MQ=255)
                                                     gAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATGCCTTATCATTCGt                                           <  1:1021957/101‑1 (MQ=18)
                                                     gAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGt                                           <  2:1340638/101‑1 (MQ=14)
                                                      aaaGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGtc                                          >  1:1425019/1‑100 (MQ=255)
                                                       aaGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGtca                                         >  1:1001687/1‑99 (MQ=255)
                                                       aaGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGtca                                         >  1:45408/1‑99 (MQ=255)
                                                       aaGCCCGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGtca                                         <  2:1636565/101‑3 (MQ=14)
                                                          cccGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTaaaa                                      >  1:634115/1‑101 (MQ=18)
                                                          cccGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCaaaa                                      <  1:219242/101‑1 (MQ=14)
                                                            cgctCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAata                                    <  1:387409/98‑1 (MQ=18)
                                                            cGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAata                                    >  2:165943/1‑101 (MQ=255)
                                                             gctCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAatat                                   >  1:1249687/3‑101 (MQ=17)
                                                              ctCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATAtt                                  >  2:617379/2‑101 (MQ=255)
                                                               tCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTa                                 <  1:1647818/101‑1 (MQ=255)
                                                               tCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTa                                 <  2:43086/101‑1 (MQ=255)
                                                               tCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTa                                 <  2:880844/101‑1 (MQ=17)
                                                                ccccGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTaa                                >  2:1231926/1‑101 (MQ=255)
                                                                 cccGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAAt                               <  1:250270/101‑1 (MQ=18)
                                                                 cccGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAAt                               <  1:505390/101‑1 (MQ=255)
                                                                 cccGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAAt                               <  2:145381/101‑1 (MQ=255)
                                                                  ccGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATa                              <  2:1620640/101‑1 (MQ=255)
                                                                  ccGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATa                              <  2:322346/101‑1 (MQ=18)
                                                                  ccGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATa                              <  2:81138/101‑1 (MQ=255)
                                                                  ccGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATa                              <  1:912815/101‑1 (MQ=255)
                                                                   cGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAg                             <  1:1462481/101‑1 (MQ=255)
                                                                   cGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAg                             >  1:409585/1‑101 (MQ=255)
                                                                   cGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAg                             >  1:1581542/1‑101 (MQ=255)
                                                                    gCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGc                            >  2:569250/1‑101 (MQ=255)
                                                                    gCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGc                            <  1:1343493/101‑1 (MQ=255)
                                                                      aaGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGa                          >  1:423810/1‑101 (MQ=255)
                                                                       aGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGAt                         <  1:780320/101‑1 (MQ=255)
                                                                        ggACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGAtg                        <  1:1490713/101‑1 (MQ=255)
                                                                        ggACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAAtg                        <  1:1282837/101‑1 (MQ=18)
                                                                         gACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAAtgt                       <  1:1063781/101‑1 (MQ=37)
                                                                         gACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAAtgt                       <  2:1045660/101‑1 (MQ=37)
                                                                          aCTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAAtgtg                      >  2:467826/1‑101 (MQ=21)
                                                                           cTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGc                     >  1:819258/1‑101 (MQ=21)
                                                                            tGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGcc                    >  1:205913/1‑101 (MQ=255)
                                                                            tGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGcc                    <  1:257859/101‑1 (MQ=21)
                                                                            tGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGcc                    >  1:1121841/1‑101 (MQ=21)
                                                                             gACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCg                   >  2:1577873/1‑101 (MQ=255)
                                                                             gACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCg                   <  1:644959/101‑1 (MQ=255)
                                                                              aCGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGt                  >  1:349413/1‑101 (MQ=21)
                                                                               cGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTa                 >  1:1486435/1‑101 (MQ=21)
                                                                                gCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTAt                <  2:352188/101‑1 (MQ=38)
                                                                                gCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTAt                >  2:697765/1‑101 (MQ=21)
                                                                                 ccAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATc               <  2:288691/101‑1 (MQ=25)
                                                                                  gaGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATcc              >  1:8016/2‑101 (MQ=255)
                                                                                   agttAGGTTCTGTCAATGGCTGCTTTTCGCATGTTACGTCTTCACCCTGCCTTGAATACTTGCTCTTTCGTCAAAATATTAGTAGCAAGGTGCCGTATccc             <  2:860938/98‑1 (MQ=16)
                                                                                    ggtAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCt            <  2:1004923/99‑1 (MQ=25)
                                                                                    ggtAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCt            >  1:168364/3‑101 (MQ=25)
                                                                                     gtAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTg           <  1:569445/100‑1 (MQ=38)
                                                                                     aTAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTg           <  1:1388264/101‑1 (MQ=255)
                                                                                      tAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGa          >  2:536223/1‑101 (MQ=37)
                                                                                      tAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGa          >  2:1516773/1‑101 (MQ=255)
                                                                                       aGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCGTATCCCTGaa         <  2:230740/101‑1 (MQ=255)
                                                                                       aGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCGTATCCCTGaa         <  1:1322122/101‑1 (MQ=255)
                                                                                        gTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGaaa        <  1:1098020/101‑1 (MQ=255)
                                                                                        gTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGaaa        >  1:530526/1‑101 (MQ=25)
                                                                                        gTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGaaa        >  2:537492/1‑101 (MQ=25)
                                                                                         tttCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGaaaa       <  1:703232/101‑1 (MQ=255)
                                                                                          ttCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAaat      <  2:1467805/101‑1 (MQ=37)
                                                                                          ttCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAaat      >  2:870620/1‑101 (MQ=37)
                                                                                          ttCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAaat      >  2:1282710/1‑101 (MQ=255)
                                                                                              |                                                                                                    
GTCATCATGTAGTGGAAGTGCCTGGAGGAGGTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATT  >  133/12604‑12798

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