Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*13312,7140TC74.9% 350.4 / 165.8 259intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base T (26/39);  new base C (87/107);  total (113/146)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 272 total aligned reads displayed.

GTGGAAGTGCCTGGAGGAGGTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCTGCTGCATTTCCTCT  >  133/12615‑12814
                                                                                                   |                                                                                                    
              aggaggTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTa                                                                                       <  2:1155505/101‑1 (MQ=255)
               ggaggTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTaa                                                                                      <  2:1314692/101‑1 (MQ=255)
               ggaggTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTaa                                                                                      <  1:457791/101‑1 (MQ=255)
                 aggTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAAcc                                                                                    >  2:296910/1‑101 (MQ=255)
                        caggcaggGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCtt                                                                             >  2:1322116/1‑101 (MQ=255)
                        aaagcaggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCtt                                                                             >  1:227341/4‑101 (MQ=11)
                            caggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAAt                                                                         <  2:1183399/101‑1 (MQ=11)
                             aggGCAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATa                                                                        <  1:133411/101‑1 (MQ=11)
                             aggGAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATa                                                                        <  2:878644/101‑1 (MQ=11)
                              gggtattactACGAAAGCCCGCTCCCCGCGAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATAc                                                                       <  1:583040/92‑1 (MQ=11)
                              gggtattactACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATAc                                                                       <  1:333930/92‑1 (MQ=11)
                              gggAAAGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATAc                                                                       >  1:1422608/1‑101 (MQ=11)
                               ggCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATAcc                                                                      >  2:908685/1‑101 (MQ=255)
                                gCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCt                                                                     >  2:776717/1‑101 (MQ=255)
                                  attactACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCttc                                                                   >  1:515364/6‑100 (MQ=14)
                                  aaGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTa                                                                   <  1:1572301/101‑1 (MQ=11)
                                   aGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTAt                                                                  <  1:397532/101‑1 (MQ=255)
                                   aGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTAt                                                                  <  2:1348768/101‑1 (MQ=14)
                                   aGATTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTAt                                                                  >  1:683584/1‑101 (MQ=255)
                                    tactACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTtctc                                                                 <  2:828567/98‑1 (MQ=11)
                                     actACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCa                                                                >  1:760160/3‑101 (MQ=14)
                                     aTTACGAAAGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCa                                                                >  2:658400/1‑101 (MQ=14)
                                     aTTACGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCa                                                                <  2:112080/101‑1 (MQ=255)
                                       tACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCAtt                                                              >  2:515816/1‑101 (MQ=255)
                                        aCGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCAtta                                                             <  2:585743/101‑2 (MQ=255)
                                        aCGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTc                                                             <  1:1258838/101‑1 (MQ=255)
                                         cGAAAGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCg                                                            >  2:272823/1‑101 (MQ=255)
                                           aaaGCCCGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGtt                                                          >  2:1719482/1‑101 (MQ=17)
                                            aaGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGtca                                                         >  1:1001687/1‑99 (MQ=255)
                                            aaGCCCGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGtca                                                         <  1:1460238/101‑3 (MQ=14)
                                            aaGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGtca                                                         <  1:1395930/101‑3 (MQ=21)
                                             aGCCCGCTCCCCGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCaa                                                        >  2:868523/1‑101 (MQ=17)
                                               cccGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTaaaa                                                      >  1:634115/1‑101 (MQ=18)
                                               cccGCTCCCCGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTaaaa                                                      >  2:1388258/1‑101 (MQ=18)
                                               cccGCTCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCaaaa                                                      <  2:1072498/101‑1 (MQ=14)
                                                 cGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAata                                                    >  2:165943/1‑101 (MQ=255)
                                                  gTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAatat                                                   >  1:209574/1‑101 (MQ=255)
                                                    tCCCCGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTa                                                 <  2:880844/101‑1 (MQ=17)
                                                      cccGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAAt                                               <  2:448018/101‑1 (MQ=255)
                                                      cccGCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAAt                                               <  1:868528/101‑1 (MQ=255)
                                                       ccGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATa                                              <  2:1620640/101‑1 (MQ=255)
                                                       ccGCAAGGACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATa                                              <  2:322346/101‑1 (MQ=18)
                                                       ccGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTATCCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATa                                              <  1:1661513/101‑1 (MQ=255)
                                                       ccGCAAGGACTGACCCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATa                                              <  1:1158279/101‑1 (MQ=18)
                                                        cGCAAGGACTGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAg                                             <  1:1683723/101‑1 (MQ=255)
                                                         gCAAGGACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGc                                            >  2:569250/1‑101 (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)
                                                             ggACTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAAtg                                        <  1:801726/101‑1 (MQ=18)
                                                             ggACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAAtg                                        <  2:1266727/101‑1 (MQ=21)
                                                              gACTGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAAtgt                                       <  1:1063781/101‑1 (MQ=37)
                                                               aCTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAAtgtg                                      <  1:367475/101‑1 (MQ=37)
                                                                cTGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGc                                     >  1:1032862/1‑101 (MQ=21)
                                                                 tGACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGcc                                    >  1:1362450/1‑101 (MQ=255)
                                                                 tGACGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGcc                                    <  1:257859/101‑1 (MQ=21)
                                                                 tGACGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGcc                                    >  1:1121841/1‑101 (MQ=21)
                                                                  gACGCGAGATAGTTTCTATCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCg                                   >  2:432070/1‑101 (MQ=37)
                                                                  gACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCg                                   >  2:1577873/1‑101 (MQ=255)
                                                                    cGCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTa                                 >  1:1486435/1‑101 (MQ=21)
                                                                    cGCCAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTa                                 <  2:899036/101‑1 (MQ=25)
                                                                     gCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTAt                                <  2:352188/101‑1 (MQ=38)
                                                                     gCCAGGTAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTAt                                >  2:697765/1‑101 (MQ=21)
                                                                      ccAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATc                               <  2:288691/101‑1 (MQ=25)
                                                                      ccAGATAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATc                               <  2:1465222/101‑1 (MQ=25)
                                                                       cAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCGTATcc                              >  2:798125/1‑101 (MQ=255)
                                                                         ggtAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCt                            <  2:1004923/99‑1 (MQ=25)
                                                                          gtAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTg                           <  1:569445/100‑1 (MQ=38)
                                                                           tAGTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGa                          >  2:536223/1‑101 (MQ=37)
                                                                             gTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGaaa                        >  1:530526/1‑101 (MQ=25)
                                                                             gTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGaaa                        >  1:195636/1‑101 (MQ=25)
                                                                             gTTTCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGaaa                        >  2:537492/1‑101 (MQ=25)
                                                                              tttCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGaaaa                       >  1:1010777/1‑101 (MQ=37)
                                                                               ttCTGTCCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAaat                      >  1:1245672/1‑101 (MQ=37)
                                                                                   gtcCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAAtt                  <  1:1082869/100‑1 (MQ=255)
                                                                                   ctcCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAAtt                  >  1:917711/1‑101 (MQ=255)
                                                                                   atcCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAAtt                  >  1:208784/2‑101 (MQ=255)
                                                                                    tcCATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTc                 <  2:728347/101‑1 (MQ=255)
                                                                                     ccATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTct                <  1:1586442/101‑1 (MQ=255)
                                                                                     ccATGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTct                <  1:895295/101‑1 (MQ=255)
                                                                                       aTGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTctgc              >  2:755944/1‑101 (MQ=255)
                                                                                       aTGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTctgc              >  1:1172934/1‑101 (MQ=255)
                                                                                        tGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTctgct             >  2:15074/1‑101 (MQ=38)
                                                                                        tGGCTGCTTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTctgct             >  1:910938/1‑101 (MQ=255)
                                                                                          gctgctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCtgctgc           <  2:260169/101‑1 (MQ=39)
                                                                                           ctgctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCtgctgcg          >  2:1525436/1‑100 (MQ=38)
                                                                                           ctgctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCTGCTGCa          >  1:104474/1‑101 (MQ=255)
                                                                                            tgctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCtgctgcgt         <  1:167792/101‑3 (MQ=255)
                                                                                            tgctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCtgcggcgt         >  1:286291/1‑99 (MQ=255)
                                                                                            tgctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCtgcggcgt         >  1:958892/1‑99 (MQ=255)
                                                                                             gctTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCTGCTGCAtt        <  2:455512/101‑1 (MQ=255)
                                                                                             gctTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTAGTCAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCTGCTGCAtt        <  2:140575/101‑1 (MQ=255)
                                                                                             gctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCTGCGGCGtt        <  1:755949/101‑1 (MQ=255)
                                                                                              ctTTTCGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCTGCTGCAttt       >  2:626419/1‑101 (MQ=255)
                                                                                               ttttcGCATCTTACGTCTTAACCCTGCCTTGAATACCTTCTCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCTGCTGCGTTTc      >  2:582295/1‑101 (MQ=255)
                                                                                               ttttcGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCTGCTGCATTTc      <  1:294903/101‑1 (MQ=255)
                                                                                                tttcGCATCTTACGTTTTAACCCTGCATTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCTGCTGCATTTcc     <  2:1637282/101‑1 (MQ=255)
                                                                                                tttcGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCTGCGGCGTTTcc     >  2:225780/1‑101 (MQ=255)
                                                                                                tttcGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTCAAAATATTAATAGCAATGTGCCGTATCCCTGAAAATAATTCTGCTGCATTTcc     >  2:397864/1‑101 (MQ=255)
                                                                                                   |                                                                                                    
GTGGAAGTGCCTGGAGGAGGTGAGCAGGCAGGGCAAGATTACGAAAGCCCGTTCCCCGCAAGGACTGACGCCAGATAGTTTCTCTCCATGGCTGCTTTTTGCATCTTACGTCTTAACCCTGCCTTGAATACCTTATCATTCGTTAAAATATTAATAGCGATGTGCCGTATCCCTGAAAATAATTCTGCTGCATTTCCTCT  >  133/12615‑12814

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