Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*6221,2320CA59.7% 58.9 / 170.1 124intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base A (38/36);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 131 total aligned reads displayed.

ATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAG  >  62/21136‑21331
                                                                                                |                                                                                                   
atatTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGaagca                                                                                                 >  2:1656308/1‑101 (MQ=255)
 tatTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGaagcag                                                                                                <  1:680087/101‑1 (MQ=255)
  atTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGa                                                                                               <  1:1344842/101‑1 (MQ=255)
  atTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGa                                                                                               >  2:997655/1‑101 (MQ=255)
   tttATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGAc                                                                                              <  1:853199/101‑1 (MQ=255)
    ttATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACa                                                                                             >  1:181964/1‑101 (MQ=255)
       tgaaatatatatAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAAt                                                                                          >  1:361004/9‑101 (MQ=17)
        gaaatatatatAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATa                                                                                         >  1:474282/8‑101 (MQ=18)
          tGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAAc                                                                                       <  1:126693/101‑1 (MQ=255)
          tGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAAc                                                                                       >  1:807548/1‑101 (MQ=255)
          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                                                                                  >  1:1083292/1‑101 (MQ=255)
               atatAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGAt                                                                                  >  1:1200613/1‑101 (MQ=38)
               atatAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGAt                                                                                  >  2:1185458/1‑101 (MQ=255)
               atatAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGAt                                                                                  <  2:270367/101‑1 (MQ=255)
                tataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATa                                                                                 >  1:829674/1‑101 (MQ=255)
                 ataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaa                                                                                >  1:200495/1‑101 (MQ=38)
                 ataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaa                                                                                >  1:1180817/1‑101 (MQ=38)
                 ataAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaa                                                                                <  1:1502439/101‑1 (MQ=38)
                  taAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATaaa                                                                               <  2:1581668/101‑1 (MQ=255)
                     atatTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACa                                                                            >  1:389643/1‑101 (MQ=37)
                       atTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAg                                                                          >  1:881829/1‑101 (MQ=25)
                         tttCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGa                                                                        >  2:994050/1‑101 (MQ=255)
                         tttCCTAATTGTTTCTATCTAATGGATATCTCACTTAAGGCTTTCTTATGAATCTGTAGGGGTTTGCCTGCAAGCGGACAAATAACCCGATAAAACAAGGa                                                                        >  2:554045/1‑101 (MQ=255)
                         tttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGa                                                                        >  2:640716/1‑101 (MQ=37)
                          ttCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGAt                                                                       <  2:498738/101‑1 (MQ=255)
                          ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGAt                                                                       >  2:851083/1‑101 (MQ=37)
                            ccTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGa                                                                     >  1:525054/1‑101 (MQ=38)
                            ccTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGa                                                                     <  1:320679/101‑1 (MQ=38)
                             cTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAg                                                                    >  1:212909/1‑101 (MQ=255)
                               aaTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCa                                                                  <  1:1114914/101‑1 (MQ=38)
                                aTTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                                 >  2:1311635/1‑101 (MQ=255)
                                aTTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                                 >  1:1345569/1‑101 (MQ=255)
                                aTTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                                 <  2:804413/101‑1 (MQ=255)
                                aTTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                                 <  2:963729/101‑1 (MQ=255)
                                aTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAg                                                                 >  2:497866/1‑101 (MQ=38)
                                 ttGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGa                                                                >  2:1120685/1‑101 (MQ=255)
                                 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:560927/1‑101 (MQ=255)
                                       cTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCg                                                          >  2:321660/1‑101 (MQ=255)
                                        tATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCgg                                                         >  1:1566591/1‑101 (MQ=255)
                                        tATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCgg                                                         <  1:1614779/101‑1 (MQ=255)
                                         atctgtcagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGa                                                        >  1:1671467/9‑101 (MQ=255)
                                         aTCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGa                                                        >  2:550133/1‑101 (MQ=255)
                                           ctgacagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaa                                                      >  1:411702/7‑101 (MQ=255)
                                           cTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaa                                                      >  1:1273717/1‑101 (MQ=255)
                                           cTGATGGATATCTCACTTAAGGCTTTCTTATAAACCTGTAGGGTTTTGCCTGGAAACAGACAAATCACCCGATAAAACAACGCTCGGCAGATGCGCGGaaa                                                      >  2:517444/1‑101 (MQ=255)
                                            tgacagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaaa                                                     >  2:380882/6‑101 (MQ=255)
                                            tGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGaaaa                                                     <  2:439779/101‑1 (MQ=255)
                                             gATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGGAAATGGGCGGAAAAc                                                    >  2:877384/1‑101 (MQ=255)
                                             gATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAAc                                                    <  2:263568/101‑1 (MQ=255)
                                              aTGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAAcc                                                   <  2:574041/101‑1 (MQ=255)
                                               cagATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATACCCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCg                                                  <  2:285957/99‑1 (MQ=255)
                                                 gATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGc                                                >  1:757216/1‑101 (MQ=255)
                                                 gATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGGTGAGCGGAAAACCGGc                                                <  2:1655210/101‑1 (MQ=255)
                                                 gATATCTCACTTAAGGCTTTCTTATAAATCTGCAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGc                                                >  2:1101191/1‑101 (MQ=255)
                                                  atatCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGcg                                               >  2:1053507/1‑100 (MQ=255)
                                                  atatCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGcg                                               <  1:1693276/101‑2 (MQ=255)
                                                   tatCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAg                                              <  2:1191257/101‑1 (MQ=255)
                                                         aCTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTc                                        >  1:444797/1‑101 (MQ=255)
                                                          cTTAAGGCTTTCTTATAAATCTGTAGTGTTTCGCCTGTCACCAGACAAATAACCCGATAAAACAAGGATGAGCAGATAAGCGAACAACCGGcgggcggtca                                       >  2:1115793/1‑92 (MQ=255)
                                                            tAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAgg                                     <  1:497871/101‑1 (MQ=255)
                                                              aGGCTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGc                                   <  2:200495/101‑1 (MQ=255)
                                                               ggCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCg                                  >  1:402055/1‑101 (MQ=255)
                                                                 cTTTCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGAc                                >  2:1087516/1‑101 (MQ=38)
                                                                  tttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACa                               >  1:1604558/1‑101 (MQ=255)
                                                                   ttCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACAt                              <  2:1581280/101‑1 (MQ=255)
                                                                   ttCTTATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACAt                              >  1:774596/1‑101 (MQ=255)
                                                                   ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACAt                              <  2:1671461/101‑1 (MQ=255)
                                                                      ttATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGAc                           <  2:952957/101‑1 (MQ=255)
                                                                      ttATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGAc                           >  1:1308271/1‑101 (MQ=255)
                                                                      ttATAAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGAc                           <  1:656826/101‑1 (MQ=255)
                                                                      ttATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGAc                           <  2:1357470/101‑1 (MQ=255)
                                                                      ttATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGAc                           <  1:1236819/101‑1 (MQ=255)
                                                                        ataAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGc                         <  1:628608/101‑1 (MQ=255)
                                                                        ataAATCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGc                         >  2:319251/1‑101 (MQ=38)
                                                                        ataAATCTGTAGGATTTTGCCTGGAAGAAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGc                         >  2:782159/1‑101 (MQ=38)
                                                                         taAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCa                        <  2:240478/101‑1 (MQ=255)
                                                                           aaTCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGt                      <  2:1247775/101‑1 (MQ=255)
                                                                           aaTCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGt                      >  1:341716/1‑101 (MQ=255)
                                                                             tCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt                    >  2:705279/1‑101 (MQ=38)
                                                                             tCTGTAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt                    <  1:9390/101‑1 (MQ=255)
                                                                               tGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATgg                  <  2:881824/101‑1 (MQ=255)
                                                                                 tAGGGTTTTGCCTGGAAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCg                <  2:1566585/101‑1 (MQ=255)
                                                                                 tAGGGATTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCg                <  2:416431/101‑1 (MQ=255)
                                                                                  aGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCgg               <  1:1653287/101‑1 (MQ=255)
                                                                                    ggTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTa             <  2:1363425/101‑1 (MQ=255)
                                                                                            cTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGtc     >  1:309626/1‑100 (MQ=255)
                                                                                            cTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGtc     <  2:173814/101‑2 (MQ=255)
                                                                                               tCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAg  <  2:444792/101‑1 (MQ=255)
                                                                                                |                                                                                                   
ATATTTATCATGGTGATATAAATATTTTCCTAATTGTTTCTATCTGATGGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAG  >  62/21136‑21331

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: ‑