Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*918,4130TG77.0% 503.0 / 152.5 252intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base G (110/84);  total (140/112)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 268 total aligned reads displayed.

AAACGCCAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTACAGAACAATCGTTAGC  >  91/8315‑8511
                                                                                                   |                                                                                                  
aaaCGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGcc                                                                                                  <  1:210864/101‑1 (MQ=9)
 accgcaaaaatgcctgatgcgctacgCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCg                                                                                                 <  2:568506‑M2/75‑1 (MQ=255)
 aaagcaaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGTATCCTGCAATATATTGAATTTGCAGGTTTTTGTAGGCCGGATAAGGCGTTCACGCCg                                                                                                 >  2:1639094/7‑101 (MQ=0)
  aagcaaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGTATCCTGCAATATATTGAATTTGCAGGTTTTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                >  2:1327361/6‑101 (MQ=0)
    gCCAAACTGCCTGAAGCGCTACGATTATCAGGACTACGCCGCTCCTGCAATATATTGAATTTTCATGCTTTTGTAGGCCGGATACGGCGATCACGCCgaaa                                                                                              >  2:583753/1‑98 (MQ=9)
     caaaattgcctgaTGCGCTACGCTTATCAGGCCTACATAATCTCTGCAATCTATTGAATTTGCATGTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                             <  2:607837‑M2/88‑1 (MQ=255)
     caaaaTTGCCTGATGCGCTACGCTTATCATGCCTACCCGGATCCTGCAATATATTTAATTTGCAGGTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                             >  2:820253/3‑101 (MQ=2)
      aaatcTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATca                                                                                            <  2:474854/97‑2 (MQ=2)
      aaaattgcctgaTGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATcc                                                                                            <  2:368073‑M2/89‑1 (MQ=255)
      aaaattgcctgaTGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATcc                                                                                            <  2:787625‑M2/89‑1 (MQ=255)
      aaaattgcctgaTGCGCTACGCTTATCAGGCCTACACGTATCCTGCAATATATTGAATTTGCAGGTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                            <  2:1526067‑M2/89‑1 (MQ=255)
      aaaaatgcctgatgcgctacgCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATcc                                                                                            <  2:327409‑M2/80‑1 (MQ=255)
      aaaaatgcctgatgcgctacgCTTATCAGGCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATcc                                                                                            <  2:1355008‑M2/80‑1 (MQ=255)
            gCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCatt                                                                                      <  2:1184852/101‑3 (MQ=14)
             ccTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCAGGCatga                                                                                     >  1:1502388/1‑98 (MQ=2)
             ccTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCAGGCatga                                                                                     <  2:137365/101‑4 (MQ=2)
               tGATGCGCTACGCTTATCAGGCCTACACGTATCCTGCAATATATTGAATTTGCAGGTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCatgaca                                                                                   <  1:393762/101‑6 (MQ=11)
                gATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCAGGCagatcgg                                                                                  >  2:696807/1‑95 (MQ=2)
                    cgcTACGCTTATCAGGCCTACACGTATCCTGCAATATATTGAATTTGCAGGTTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCatgacaaagcg                                                                              >  2:1676388/1‑91 (MQ=0)
                        aCGCTTATCAGGCCTACATGATCTCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAagaaacagcaaaca                                                                          <  1:1316905‑M2/101‑15 (MQ=255)
                           cTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                                                                       <  1:563016/101‑1 (MQ=255)
                           cTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                                                                       <  1:370749/101‑1 (MQ=255)
                            ttATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTg                                                                      >  2:250067/1‑101 (MQ=255)
                                  ggCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTc                                                                >  2:1432516/1‑101 (MQ=255)
                                   gCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTcc                                                               >  2:63395/1‑101 (MQ=255)
                                     cTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCa                                                             >  1:420186/1‑101 (MQ=255)
                                        acatgatctctGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcg                                                           <  2:974362‑M2/96‑30 (MQ=255)
                                         gCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACgtt                                                         <  2:155587/101‑1 (MQ=255)
                                          gaggatggtgcaatataTTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgag                                                         <  1:283861‑M2/93‑38 (MQ=255)
                                          atgatctCTGCAATATATTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgtt                                                         <  1:835933‑M2/98‑32 (MQ=255)
                                          aaatctattgcaacatgTTGAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgctt                                                         >  1:53818‑M2/9‑70 (MQ=255)
                                          aaatctattgcaacatgTTGAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgctt                                                         >  1:600201‑M2/9‑70 (MQ=255)
                                           tgatcTCTGCAATATATTGAATTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaacc                                                        >  1:1515774‑M2/3‑69 (MQ=255)
                                            gatCTCTGCAATATATTGAATTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccc                                                       >  2:1661215‑M2/2‑68 (MQ=255)
                                             gatggtgcaatATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtag                                                      >  1:954302‑M2/6‑61 (MQ=255)
                                               tggtgcaATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGGTTCCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagct                                                    >  2:1093380‑M2/4‑59 (MQ=255)
                                                attgcAACATGTTGAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcagg                                                   >  1:151006‑M2/3‑64 (MQ=255)
                                                 gtgCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctca                                                  >  1:1181892‑M2/2‑57 (MQ=255)
                                                  tGCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaa                                                 <  1:1432110‑M2/101‑46 (MQ=255)
                                                  tGCAACATGTTGAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcc                                                 <  1:1156842‑M2/101‑40 (MQ=255)
                                                   gCAATATATTGAATTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttctt                                                >  1:1546572‑M2/1‑61 (MQ=255)
                                                   gCAATATATTGAATTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttctt                                                <  1:580972‑M2/101‑41 (MQ=255)
                                                   gCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcagcaaactgagtagctcaag                                                >  1:1449413‑M2/1‑55 (MQ=255)
                                                   gCAATATATTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaag                                                >  2:1047781‑M2/1‑55 (MQ=255)
                                                    cAATATAGTAAATTTACGTACTTTTGTAGGCCGGATAAGTCGTTCACGCCGCATCCGGCACGCTGTGCCAGTAATCTGAAACTTCCACCTTGTTCgggtgt                                               <  2:1174356/101‑4 (MQ=21)
                                                     aaTATATTGAGTTTGCGTTTTTTTGTAGGCCGTATAAGGGTTTCCGGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgttt                                              >  2:905751‑M2/1‑59 (MQ=255)
                                                      atatatTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgtttt                                             <  1:1364807‑M2/101‑44 (MQ=255)
                                                      atatatTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgtttt                                             >  2:835928‑M2/1‑58 (MQ=255)
                                                      atatatTGAGTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgtttt                                             >  1:1148350‑M2/1‑58 (MQ=255)
                                                      atatatTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaa                                             <  2:1535637‑M2/101‑50 (MQ=255)
                                                      acatgttgAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctaca                                             >  1:1112565‑M2/3‑58 (MQ=255)
                                                       catgttGAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaa                                            >  1:932528‑M2/2‑57 (MQ=255)
                                                        atatTGAATTTAAGCGATTTTGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaat                                           >  1:1056476‑M2/1‑50 (MQ=255)
                                                         tatTGAATTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttcttttaaag                                          >  2:397713‑M2/1‑55 (MQ=255)
                                                          gttgAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaat                                         >  1:1046170‑M2/2‑54 (MQ=255)
                                                          gttgAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaat                                         >  1:356571‑M2/2‑54 (MQ=255)
                                                          gttgAATCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaat                                         >  1:682498‑M2/2‑54 (MQ=255)
                                                           ttGAATTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttcttttaaagag                                        >  1:1206624‑M2/1‑53 (MQ=255)
                                                            tGAATTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttcttttaaagagg                                       <  2:1717404‑M2/101‑50 (MQ=255)
                                                              aaTCTTCGTGCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaatctat                                     >  2:776018‑M2/1‑50 (MQ=255)
                                                               aTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGc                                    <  2:211743/101‑1 (MQ=255)
                                                               aTTTGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttcttttaaagagggta                                    >  1:1015122‑M2/1‑49 (MQ=255)
                                                                tttGCGTGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttcttttaaagagggtat                                   >  1:516850‑M2/1‑48 (MQ=255)
                                                                tcttcgtgcGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaatctattg                                   >  1:817447‑M2/5‑48 (MQ=255)
                                                                 ttaagcgattttgtaggccgGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaa                                  >  2:1570431‑M2/9‑41 (MQ=255)
                                                                 cttcgtgCGTTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaatctattgc                                  >  1:325732‑M2/4‑47 (MQ=255)
                                                                 aatattgctttTGTAGGCCGGATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcat                                  >  2:459781‑M2/6‑39 (MQ=255)
                                                                   tattgctTTTGTAGGCCGGATAAGCCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattg                                >  1:706681‑M2/4‑37 (MQ=255)
                                                                   aagcgattttgtaggcCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaaga                                >  1:1405161‑M2/7‑39 (MQ=255)
                                                                   aagcgattttgtaggcCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaaga                                <  1:78389‑M2/95‑63 (MQ=255)
                                                                     ttgCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcg                              >  1:17222‑M2/2‑35 (MQ=255)
                                                                      tGATTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtgaacaaagcgcattttgtcagcaatctaaccctcttcttttaaagagggtattgatca                             >  2:281386‑M2/1‑42 (MQ=255)
                                                                        cTTTTGTAGGCCGGATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggca                           <  2:1074338‑M2/101‑70 (MQ=255)
                                                                         ttttGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaagaccaaca                          >  1:772363‑M2/1‑33 (MQ=255)
                                                                         ttttGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaagaccaaca                          <  1:461715‑M2/101‑69 (MQ=255)
                                                                          tttGTAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaagaccaacag                         <  2:650367‑M2/101‑70 (MQ=255)
                                                                          tttGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaatctattgcaacatgttg                         >  2:1601882‑M2/1‑38 (MQ=255)
                                                                           ttGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACAttt                        >  1:376397/1‑101 (MQ=255)
                                                                             gTAGGCCGGATAAGGTGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgttttttctgcttttcttcgcgaattaa                      >  1:1269384‑M2/1‑35 (MQ=255)
                                                                             gTAGGCCGGATAAGGCGTTCACGACGCATCCGGCAttaggtgctcaatgcctgatgcgctacgcttatcaggcctacaaaatctattgcaacatgttgaat                      >  2:1122715‑M2/1‑35 (MQ=255)
                                                                              tAGGCCGGATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagc                     <  2:353591‑M2/101‑76 (MQ=255)
                                                                              tAGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaagaccaacaggagt                     >  1:574524‑M2/1‑28 (MQ=255)
                                                                               aGGCCGGATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagca                    <  1:1662556‑M2/101‑77 (MQ=255)
                                                                               aGGCCGGATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagca                    <  1:1578486‑M2/101‑77 (MQ=255)
                                                                               aGGCCGGATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagca                    <  2:863684‑M2/101‑77 (MQ=255)
                                                                               aGGCCGGATAAGGCGTTCACGCCGCATccggcaaaaacaacgaacactttgtcaacaaactgagtagctcaaggaaatcccaatgaagaccaacaggagtc                    <  1:1563139‑M2/101‑75 (MQ=255)
                                                                                ggCCGGATAAGGCGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgttttttctgcttttcttcgcgaattaattc                   >  2:1316899‑M2/1‑32 (MQ=255)
                                                                                ggCCGAATAAGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagcac                   >  2:550093‑M2/1‑24 (MQ=255)
                                                                                 gCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTa                  <  2:1158730/101‑1 (MQ=255)
                                                                                 gCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTa                  <  2:737865/101‑1 (MQ=255)
                                                                                   cGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTACa                <  2:6206/101‑1 (MQ=255)
                                                                                   cGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTACa                <  2:1281025/101‑1 (MQ=255)
                                                                                    ggATAAGGCGTGCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagcacggca               >  1:1096103‑M2/1‑20 (MQ=255)
                                                                                       aaagGCGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgttttttctgcttttcttcgcgaattaattccgcttcg            >  1:75620‑M2/2‑25 (MQ=255)
                                                                                        aaGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagcacggcagcaa           <  2:52285‑M2/101‑86 (MQ=255)
                                                                                         aGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTACAGAACaa          <  1:1646252/101‑1 (MQ=255)
                                                                                         aGGCGTTCACGCCGCatccggcatgaacaaagagcacttcatcagcaatctgaaactcaactaaacagcccaatcattgcggcacacagcacggcagcaac          >  2:1396586‑M2/1‑15 (MQ=255)
                                                                                          ggCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTACAGAACAAt         >  1:99211/1‑101 (MQ=255)
                                                                                            cGTTCACGCCGCATCCGGCAagaaacagcaaacaatccaaaacgccgcgttcagcggcgttttttctgcttttcttcgcgaattaattccgcttcgcaatt       >  2:1001390‑M2/1‑20 (MQ=255)
                                                                                             gTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTACAGAACAATCGt      >  1:243458/1‑101 (MQ=255)
                                                                                                   |                                                                                                  
AAACGCCAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCTGGACATTTTCTCTACAGAACAATCGTTAGC  >  91/8315‑8511

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