Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*6221,2250CT64.8% 132.4 / 123.3 128intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (44/39);  total (67/61)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 133 total aligned reads displayed.

TAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAT  >  62/21128‑21324
                                                                                                 |                                                                                                   
taaatCACATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGcc                                                                                                  >  2:656821/1‑101 (MQ=255)
 aaatCACATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCt                                                                                                 >  2:106502/1‑101 (MQ=255)
   atCACATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTgg                                                                                               <  1:539083/101‑2 (MQ=255)
    tCACATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTgga                                                                                              >  1:291241/1‑99 (MQ=255)
    tCACATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTgga                                                                                              <  1:539549/101‑3 (MQ=255)
       caTATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTggaagc                                                                                           >  1:77269/1‑96 (MQ=255)
        atatTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGaagca                                                                                          >  2:1656308/1‑101 (MQ=255)
         tatTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGaagcag                                                                                         <  1:680087/101‑1 (MQ=255)
         tatTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGaagcag                                                                                         <  1:683695/101‑1 (MQ=255)
          atTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGa                                                                                        >  2:997655/1‑101 (MQ=255)
          atTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGa                                                                                        <  1:1344842/101‑1 (MQ=255)
            ttATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACa                                                                                      >  1:181964/1‑101 (MQ=255)
               tgaaatatatatAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAAt                                                                                   >  1:361004/9‑101 (MQ=17)
                  tGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAAc                                                                                >  1:807548/1‑101 (MQ=255)
                  tGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAAc                                                                                <  1:126693/101‑1 (MQ=255)
                  aatatatatAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAAc                                                                                <  2:527989/96‑1 (MQ=21)
                  aatatatatAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAAc                                                                                <  1:482056/96‑1 (MQ=18)
                    gTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAAccc                                                                              <  1:82890/101‑1 (MQ=38)
                    gTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAAccc                                                                              <  1:1057908/101‑1 (MQ=38)
                     tGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCg                                                                             <  1:83784/101‑1 (MQ=38)
                       atatAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGAt                                                                           <  2:270367/101‑1 (MQ=255)
                       atatAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGAt                                                                           >  1:1083292/1‑101 (MQ=255)
                       atatAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGAt                                                                           >  2:1185458/1‑101 (MQ=255)
                       atatAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGAt                                                                           >  1:1200613/1‑101 (MQ=38)
                        tataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATa                                                                          >  1:829674/1‑101 (MQ=255)
                         ataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaa                                                                         <  1:1502439/101‑1 (MQ=38)
                         ataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaa                                                                         >  1:1180817/1‑101 (MQ=38)
                         ataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaa                                                                         >  1:200495/1‑101 (MQ=38)
                          taAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaaa                                                                        <  2:1581668/101‑1 (MQ=255)
                             atatTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACa                                                                     >  1:389643/1‑101 (MQ=37)
                               atTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAACACAAg                                                                   >  2:1252474/1‑101 (MQ=39)
                               atTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAg                                                                   >  1:480282/1‑101 (MQ=37)
                               atTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAg                                                                   >  1:881829/1‑101 (MQ=25)
                                 tttCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGa                                                                 >  2:994050/1‑101 (MQ=255)
                                  ttCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGAt                                                                <  2:1061470/101‑1 (MQ=255)
                                  ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGAt                                                                >  2:851083/1‑101 (MQ=37)
                                    ccTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGa                                                              <  1:320679/101‑1 (MQ=38)
                                     cTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAg                                                             >  1:212909/1‑101 (MQ=255)
                                       aaTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCa                                                           <  1:1114914/101‑1 (MQ=38)
                                        aTTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                          <  2:963729/101‑1 (MQ=255)
                                        aTTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                          <  2:804413/101‑1 (MQ=255)
                                        aTTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                          >  1:1345569/1‑101 (MQ=255)
                                        aTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCTTGTCAGCAGACAAATGACCCGATATAACAAGGATGAGCAg                                                          >  1:855110/1‑101 (MQ=37)
                                        aTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                          >  2:497866/1‑101 (MQ=38)
                                         ttGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGa                                                         >  2:1164779/1‑101 (MQ=38)
                                          tGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGAt                                                        >  2:1200273/1‑101 (MQ=255)
                                           gTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATg                                                       >  1:1045900/1‑101 (MQ=255)
                                            tttCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGa                                                      <  2:1083287/101‑1 (MQ=255)
                                            ttcttATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGa                                                      <  1:1262258/97‑1 (MQ=255)
                                            ttcttATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGa                                                      >  1:416436/5‑101 (MQ=255)
                                             ttCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAg                                                     <  2:654984/101‑1 (MQ=255)
                                               ttATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCg                                                   <  2:205521/100‑1 (MQ=255)
                                               cTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCg                                                   >  2:321660/1‑101 (MQ=255)
                                               cTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCg                                                   >  2:560927/1‑101 (MQ=255)
                                                tATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCgg                                                  >  1:1566591/1‑101 (MQ=255)
                                                tATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAACGATGAGCAGATCAGCgg                                                  >  2:619620/1‑101 (MQ=255)
                                                 atctgtcagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGa                                                 >  1:1671467/9‑101 (MQ=255)
                                                 aTCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGa                                                 >  2:550133/1‑101 (MQ=255)
                                                   cTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaa                                               >  1:1273717/1‑101 (MQ=255)
                                                   cTGATGGATATCTCACTTAAGGCTTTCTTATAAACCTGTAGGGTTTTGCCTGGAAACAGACAAATCACCCGATAAAACAACGCTCGGCAGATGCGCGGaaa                                               >  2:517444/1‑101 (MQ=255)
                                                    tgacagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaaa                                              <  2:1552962/96‑1 (MQ=255)
                                                    tgacagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaaa                                              >  2:380882/6‑101 (MQ=255)
                                                    tGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaaa                                              <  1:414521/101‑1 (MQ=255)
                                                    tGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaaa                                              <  2:439779/101‑1 (MQ=255)
                                                     gacagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAAc                                             <  2:316655/97‑1 (MQ=255)
                                                     gATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGGAAATGGGCGGAAAAc                                             >  2:877384/1‑101 (MQ=255)
                                                      aTGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAAcc                                            <  2:574041/101‑1 (MQ=255)
                                                       cagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATACCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCg                                           <  2:285957/99‑1 (MQ=255)
                                                        ggATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCgg                                          <  2:1110001/101‑1 (MQ=255)
                                                         gATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGGTGAGCGGAAAACCGGc                                         <  2:1655210/101‑1 (MQ=255)
                                                         gATATCTCACTTAAGGCTTTCTTATAAATCTGCAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGc                                         >  2:1101191/1‑101 (MQ=255)
                                                          atatCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGcg                                        <  1:1693276/101‑2 (MQ=255)
                                                                 aCTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTc                                 >  1:444797/1‑101 (MQ=255)
                                                                  cTTAAGGCTTTCTTATAAATCTGTAGTGTTTCGCCTGTCACCAGACAAATAACCCGATAAAACAAGGATGAGCAGATAAGCGAACAACCGGcgggcggtca                                >  2:1115793/1‑92 (MQ=255)
                                                                   ttAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAg                               <  1:625136/101‑1 (MQ=255)
                                                                    tAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAgg                              <  1:497871/101‑1 (MQ=255)
                                                                      aGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGc                            <  2:200495/101‑1 (MQ=255)
                                                                       ggCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCg                           >  1:402055/1‑101 (MQ=255)
                                                                          tttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACa                        >  1:1604558/1‑101 (MQ=255)
                                                                           ttCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACAt                       >  1:774596/1‑101 (MQ=255)
                                                                           ttCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACAt                       <  2:1581280/101‑1 (MQ=255)
                                                                           ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACAt                       <  2:1671461/101‑1 (MQ=255)
                                                                              ttATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGAc                    <  1:656826/101‑1 (MQ=255)
                                                                              ttATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGAc                    >  1:1308271/1‑101 (MQ=255)
                                                                              ttATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGAc                    <  2:952957/101‑1 (MQ=255)
                                                                              ttATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGAc                    <  1:1236819/101‑1 (MQ=255)
                                                                              ttATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGAc                    <  2:1357470/101‑1 (MQ=255)
                                                                                ataAATCTGTAGGATTTTGCCTGGAAGAAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGc                  >  2:782159/1‑101 (MQ=38)
                                                                                 taAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCa                 <  2:240478/101‑1 (MQ=255)
                                                                                  aaaTCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAg                >  2:458115/1‑101 (MQ=255)
                                                                                   aaTCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGt               >  1:341716/1‑101 (MQ=255)
                                                                                   aaTCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGt               <  2:1247775/101‑1 (MQ=255)
                                                                                     tCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt             <  1:9390/101‑1 (MQ=255)
                                                                                     tCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt             >  2:705279/1‑101 (MQ=38)
                                                                                       tGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATgg           <  2:881824/101‑1 (MQ=255)
                                                                                        gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGc          >  2:1166578/1‑101 (MQ=255)
                                                                                         tAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCg         <  2:1566585/101‑1 (MQ=255)
                                                                                          aGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCgg        <  1:1653287/101‑1 (MQ=255)
                                                                                            ggTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTa      <  2:1363425/101‑1 (MQ=255)
                                                                                                tCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAt  <  2:389638/101‑1 (MQ=255)
                                                                                                 |                                                                                                   
TAAATCACATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAT  >  62/21128‑21324

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: