Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*8719,4780CT73.4% 258.3 / 109.3 158intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (19/23);  new base T (65/51);  total (84/74)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 161 total aligned reads displayed.

GTAGCAGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGC‑‑‑‑‑‑‑‑‑‑  >  87/19378‑19524
                                                                                                    |                                                        
gtAGCAGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGc                                                          <  1:1322989/101‑1 (MQ=255)
 tAGCAGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCt                                                         >  2:480557/1‑101 (MQ=255)
  aGCAGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCtc                                                        <  2:208046/101‑1 (MQ=255)
  aGCAGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCtc                                                        >  2:458843/1‑101 (MQ=255)
   gCAGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCtca                                                       <  2:1639545/101‑1 (MQ=255)
     aGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCtcatc                                                     >  1:1678240/1‑101 (MQ=255)
        aCAAATTTACGAAAGCCTGCTCCCTGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAgg                                                  <  1:480319/101‑1 (MQ=255)
              ttACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATc                                            >  2:1093489/1‑101 (MQ=255)
               tACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCg                                           <  1:1160855/101‑1 (MQ=255)
                     gttctgctattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACa                                     >  1:1007001‑M2/66‑101 (MQ=255)
                     gttctgctattatggaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACa                                     >  2:417029‑M2/68‑101 (MQ=255)
                      ttctgctattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACaa                                    >  1:967183‑M2/65‑101 (MQ=255)
                      ttctgctattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACaa                                    >  2:1326004‑M2/65‑101 (MQ=255)
                      ttctgctattatggaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACaa                                    >  1:757095‑M2/67‑101 (MQ=255)
                        ctgctattattgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAata                                  >  2:239098‑M2/65‑101 (MQ=255)
                        ctgctattattgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAata                                  >  1:780940‑M2/65‑101 (MQ=255)
                        ctgctattatggaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAata                                  <  2:631414‑M2/37‑1 (MQ=255)
                        cTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAata                                  <  1:834419/101‑1 (MQ=255)
                        cTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAata                                  >  1:231767/1‑101 (MQ=255)
                         tgctattattgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAatat                                 >  1:1423832‑M2/64‑101 (MQ=255)
                         tgctattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAAcaatgt                                 <  1:6125‑M2/40‑3 (MQ=255)
                         tgctattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAAcaatgt                                 <  2:117965‑M2/40‑3 (MQ=255)
                         tgctattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAAcaatgt                                 <  1:1576187‑M2/40‑3 (MQ=255)
                         tgctattatggaccgaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAatat                                 <  2:665675‑M2/38‑1 (MQ=255)
                          gctattattgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATg                                >  1:290046‑M2/63‑101 (MQ=255)
                          gctattattgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATg                                <  1:1436831‑M2/39‑1 (MQ=255)
                           ctattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGa                               >  2:66369‑M2/60‑101 (MQ=255)
                           cTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGa                               >  1:415049/1‑101 (MQ=255)
                           cTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGa                               >  2:1368493/1‑101 (MQ=255)
                           cTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGa                               <  2:644885/101‑1 (MQ=255)
                            tattattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGaa                              >  2:1645022‑M2/59‑101 (MQ=255)
                             ccccGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAAGCAATAAGCTCATCAGGGACATCGCAAACAATATGAAg                             <  1:1432197/101‑1 (MQ=255)
                             attatggaccgaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAg                             >  1:436454‑M2/60‑101 (MQ=255)
                             attatggaccgaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAg                             <  1:779115‑M2/42‑1 (MQ=255)
                             attatggaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAg                             <  1:564971‑M2/42‑1 (MQ=255)
                              ttattgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATAAATAAGTTCATCAGGGACATCGCAAACAATATGAAGa                            >  1:200196‑M2/59‑101 (MQ=255)
                              ttattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGa                            >  1:965026‑M2/57‑101 (MQ=255)
                              ttatggaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGa                            >  2:1098102‑M2/59‑101 (MQ=255)
                              cccGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGa                            >  2:1279438/1‑101 (MQ=255)
                              atattgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGa                            >  2:1595826‑M2/59‑101 (MQ=255)
                               tattgaccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGAc                           <  1:5948‑M2/46‑1 (MQ=255)
                               tatggaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGAc                           <  1:391079‑M2/44‑1 (MQ=255)
                                 tggaccgaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACgg                         <  1:262603‑M2/46‑1 (MQ=255)
                                  tgaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGAcgga                        >  2:1493826‑M2/55‑100 (MQ=255)
                                  cAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACggg                        <  2:1326787/101‑1 (MQ=255)
                                   gaccgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGacggat                       <  2:1657447‑M2/48‑3 (MQ=255)
                                   gaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGacggat                       <  1:156720‑M2/48‑3 (MQ=255)
                                   gaccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGacggat                       <  2:265585‑M2/48‑3 (MQ=255)
                                    accgaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGAtt                      <  1:483648‑M2/49‑1 (MQ=255)
                                    accgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGAtt                      >  1:1128387‑M2/53‑101 (MQ=255)
                                     ggACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGttt                     <  2:922821/101‑1 (MQ=255)
                                     ccgaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGAttt                     <  2:39195‑M2/52‑1 (MQ=255)
                                     ccgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGAttt                     >  2:871277‑M2/52‑101 (MQ=255)
                                      cgaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttc                    >  1:697420‑M2/51‑101 (MQ=255)
                                      cgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttc                    <  1:870625‑M2/51‑1 (MQ=255)
                                      cgaaaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttc                    >  2:792193‑M2/51‑101 (MQ=255)
                                      cgaaaggagactgccacgcataatttctccagccctttccattcaaacgtGAAATCAATAAGTACATCAGGGACATCGCAAACAATATGAAGACGGATttc                    <  2:533406‑M2/51‑1 (MQ=255)
                                       gaaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttct                   <  2:1473867‑M2/52‑1 (MQ=255)
                                       gaaaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttct                   >  1:866784‑M2/50‑101 (MQ=255)
                                       gaaaggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATttct                   >  2:629765‑M2/48‑101 (MQ=255)
                                        cTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTttctt                  >  1:752943/1‑101 (MQ=255)
                                        aaaggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttctt                  >  1:1077777‑M2/49‑101 (MQ=255)
                                         tGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTttcttc                 <  2:1512896/101‑1 (MQ=255)
                                         aaggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 >  1:1374260‑M2/48‑101 (MQ=255)
                                         aaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 <  1:1035360‑M2/54‑1 (MQ=255)
                                         aaggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttc                 >  2:607614‑M2/48‑101 (MQ=255)
                                          aggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:206541‑M2/55‑1 (MQ=255)
                                          aggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                >  2:1304349‑M2/47‑101 (MQ=255)
                                          aggagactgccacgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:1156362‑M2/55‑1 (MQ=255)
                                          aggagactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:1665166‑M2/55‑1 (MQ=255)
                                          aggagactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTtcttct                >  2:1150829‑M2/45‑101 (MQ=255)
                                          aggagactgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttct                <  1:1536446‑M2/55‑1 (MQ=255)
                                            cGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTctct              <  1:322428/101‑1 (MQ=255)
                                            cGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTctct              <  2:1612132/101‑1 (MQ=255)
                                             gTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTg             <  1:1003762/101‑1 (MQ=255)
                                             agactgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTg             >  1:1501981‑M2/44‑101 (MQ=255)
                                             agactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTg             >  1:1275765‑M2/42‑101 (MQ=255)
                                              gactgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTGc            >  2:1650876‑M2/41‑101 (MQ=255)
                                               actgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTCTgcc           <  2:79660‑M2/62‑2 (MQ=255)
                                               actgccacgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTgcc           >  2:97019‑M2/42‑100 (MQ=255)
                                                atgccacgcataatttcttcagtcctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCtgccg          <  1:523366‑M2/61‑3 (MQ=255)
                                                aGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccg          >  1:239988/1‑99 (MQ=255)
                                                aGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccg          >  2:1436001/1‑99 (MQ=255)
                                                 tgccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCTctgccgt         >  1:1032233‑M2/38‑98 (MQ=255)
                                                  gccacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTTCtctgccgtg        <  2:305946‑M2/65‑5 (MQ=255)
                                                    cacgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCTtctctgccgtgac      >  2:1238466‑M2/35‑95 (MQ=255)
                                                     gTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGccgtgact     <  1:184086/101‑8 (MQ=255)
                                                     acgcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCttctctgccgtgact     <  2:1035343‑M2/66‑8 (MQ=255)
                                                     acgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCttctctgccgtgact     >  1:1435330‑M2/34‑94 (MQ=255)
                                                     acgcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTTCttctctgccgtgact     >  1:1201447‑M2/34‑94 (MQ=255)
                                                      cgcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTcttctctgccgtgactc    <  1:1468830‑M2/67‑9 (MQ=255)
                                                       gcataatttctttagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttctctgccgtgactct   >  2:565134‑M2/34‑92 (MQ=255)
                                                       gcataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATTtcttctctgccgtgactct   >  1:1241548‑M2/32‑92 (MQ=255)
                                                       gcataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTtcttctctgccgtgactct   >  1:1518419‑M2/34‑92 (MQ=255)
                                                        cataatttctttagccctttcctttcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATatcttctctgccgtgactctt  >  2:1653265‑M2/33‑91 (MQ=255)
                                                        cataatttcttcagtcctttccattcaaacgtAAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  >  1:1314948‑M2/33‑91 (MQ=255)
                                                        cataatttcttcagccctttccattcaaatGTGAAATCAATAAGTTCATCAGGGACATCGCAAACAATGTGAAGACGGATttcttctctgccgtgactctt  >  1:518901‑M2/31‑91 (MQ=255)
                                                        cataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  >  2:1156989‑M2/33‑91 (MQ=255)
                                                        cataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  >  1:1117617‑M2/33‑91 (MQ=255)
                                                        cataatttcttcagccctttccattcaaacgtGAAATCAATAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATttcttctctgccgtgactctt  <  2:911084‑M2/69‑11 (MQ=255)
                                                                                                    |                                                        
GTAGCAGGACAAATTTACGAAAGCCTGCTCCCCGCAAGGACTGACGTCAGGTAGTTTCTAACCATAGCTGCCCTTTCTATTCAAACGTGAAATCAATAAGCTCATCAGGGACATCGCAAACAATATGAAGACGGGTTTCTTCTCTGC‑‑‑‑‑‑‑‑‑‑  >  87/19378‑19524

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