Predicted mutation
evidence seq id position mutation annotation gene description
RA 173 30 2 bp→GC intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*173300AG86.2% 376.6 / 42.2 145intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base A (11/9);  new base G (67/58);  total (78/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01
*173310AC86.2% 431.5 / 45.5 145intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base A (11/9);  new base C (67/58);  total (78/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 151 total aligned reads displayed.

‑‑‑‑‑‑‑‑‑ACATCTAACGCCCGAGCCGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGG  >  173/1‑131
                                      ||                                                                                                    
ccgcatccgacatctaacgCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggat                                         <  1:597331‑M1/92‑41 (MQ=255)
ccgcatccgacatccaacgTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggat                                         <  2:16588‑M1/92‑41 (MQ=255)
ccgcatccgacatccaacgTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggat                                         <  1:1307640‑M1/92‑41 (MQ=255)
 cgcatccgacatctaacGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggata                                        <  1:374840‑M1/93‑42 (MQ=255)
 cgcatccgacatctaacGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggata                                        >  1:1366277‑M1/9‑60 (MQ=255)
  gcatccgacatctaaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataa                                       >  1:429336‑M1/8‑59 (MQ=255)
  gcatccgacatctaaCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataa                                       >  1:59531‑M1/8‑59 (MQ=255)
  gcatccgacatccaaCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataa                                       <  1:50944‑M1/94‑43 (MQ=255)
  gcatccgaCATCTAACGCCCGAGCCGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCAcat                                       >  2:243210/8‑101 (MQ=255)
  gcatccgaCATCTAACGCCCGAGCCGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCAcat                                       <  1:1306973/94‑1 (MQ=255)
   catccgacatctaACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataag                                      <  1:706363‑M1/95‑44 (MQ=255)
   catccgacatctaACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataag                                      >  2:1701354‑M1/7‑58 (MQ=255)
   catccgacatctaACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataag                                      <  2:1098797‑M1/95‑44 (MQ=255)
   catccgacatctaACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataag                                      <  1:152270‑M1/95‑44 (MQ=255)
   catccgacatccaACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataag                                      >  1:1385401‑M1/7‑58 (MQ=255)
    atccgacatctAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataagg                                     >  1:501776‑M1/6‑57 (MQ=255)
    atccgacatctAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataagg                                     >  2:788388‑M1/6‑57 (MQ=255)
    atccgacatctAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataagg                                     >  2:886949‑M1/6‑57 (MQ=255)
    atccgacatctAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataagg                                     >  1:1562459‑M1/6‑57 (MQ=255)
    atccgacatctAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataagg                                     >  2:707788‑M1/6‑57 (MQ=255)
    atccgacatctAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataagg                                     >  1:631831‑M1/6‑57 (MQ=255)
    atccgacatccAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataagg                                     >  1:1152555‑M1/6‑57 (MQ=255)
    atccgacatccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataagg                                     >  2:1407147‑M1/6‑57 (MQ=255)
      ccgacatCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcg                                   >  1:1078607‑M1/4‑55 (MQ=255)
       cgacaTCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcgg                                  >  2:1540420‑M1/3‑54 (MQ=255)
       cgacaTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcgg                                  >  2:1435195‑M1/3‑54 (MQ=255)
        gacATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggt                                 >  1:264893‑M1/2‑53 (MQ=255)
        gacATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggt                                 <  1:1172263‑M1/100‑49 (MQ=255)
        gacATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggt                                 <  2:957283‑M1/100‑49 (MQ=255)
         aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggtt                                >  2:1446590‑M1/1‑52 (MQ=255)
         aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggtt                                <  2:7773‑M1/101‑50 (MQ=255)
         aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggtt                                >  1:1659466‑M1/1‑52 (MQ=255)
         aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggtt                                <  2:1464262‑M1/101‑50 (MQ=255)
         aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggtt                                >  2:1302019‑M1/1‑52 (MQ=255)
         aCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggtt                                <  2:997255‑M1/101‑50 (MQ=255)
          cATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctcctccccgtaggttggataaggcggtta                               <  1:1416048‑M1/101‑51 (MQ=255)
           aTCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttac                              <  2:1232721‑M1/101‑52 (MQ=255)
           aTCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttac                              <  2:915805‑M1/101‑52 (MQ=255)
           aTCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttac                              <  1:881807‑M1/101‑52 (MQ=255)
           aTCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttac                              <  2:1548528‑M1/101‑52 (MQ=255)
           aTCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttac                              >  2:1265991‑M1/1‑50 (MQ=255)
           aTCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttac                              >  1:1232726‑M1/1‑50 (MQ=255)
           aTCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttac                              <  2:1338398‑M1/101‑52 (MQ=255)
            tccaacgtCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacg                             >  2:1416042‑M1/4‑49 (MQ=255)
            tccaacgcCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacg                             <  2:875184‑M1/98‑53 (MQ=255)
            tccaacgcCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacg                             <  1:761505‑M1/98‑53 (MQ=255)
            tCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacg                             >  1:741102‑M1/1‑49 (MQ=255)
            tCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacg                             >  1:17588‑M1/1‑49 (MQ=255)
             ccaacgTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgc                            >  1:1015533‑M1/3‑48 (MQ=255)
             cTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgc                            >  2:209807‑M1/1‑48 (MQ=255)
              tAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttacgcc                           <  2:631826‑M1/101‑55 (MQ=255)
              caacGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgcc                           >  1:16588‑M1/2‑47 (MQ=255)
              caacGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgcc                           <  1:181715‑M1/100‑55 (MQ=255)
                aCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttacgccgc                         >  2:1072190‑M1/1‑45 (MQ=255)
                aCGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgc                         >  1:206288‑M1/1‑45 (MQ=255)
                 cgtccgagCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggccgcgtaaggcggttacgccgca                        >  2:1141430‑M1/4‑44 (MQ=255)
                 cGCCCGAGCCGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAata                        <  1:8453/101‑1 (MQ=255)
                 cGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttacgccgca                        >  2:8453‑M1/1‑44 (MQ=255)
                  gtccgaGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgccgcat                       >  1:790425‑M1/3‑43 (MQ=255)
                  gCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcat                       <  2:839630‑M1/101‑59 (MQ=255)
                   cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggatacgccgcatc                      >  1:784511‑M1/1‑42 (MQ=255)
                   cccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagagctccaccccgtaggtcggataaggcggttacgccgcatc                      >  1:1643375‑M1/1‑42 (MQ=255)
                   cccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatc                      <  1:387235‑M1/101‑60 (MQ=255)
                   cccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatc                      <  2:1468091‑M1/101‑60 (MQ=255)
                   cccAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatc                      >  1:338545‑M1/1‑42 (MQ=255)
                    ccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgccgcatcc                     <  1:737581‑M1/101‑61 (MQ=255)
                    ccGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgccgcatcc                     <  1:85755‑M1/101‑61 (MQ=255)
                    ccGAGCCGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATcc                     >  1:620272/1‑101 (MQ=255)
                     caagcCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttacgccgcatccg                    <  1:943997‑M1/99‑62 (MQ=255)
                      gAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgccgcatccga                   >  1:759234‑M1/1‑39 (MQ=255)
                       aGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgac                  >  1:1300825‑M1/1‑38 (MQ=255)
                        gCCGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCa                 >  2:217581/1‑101 (MQ=255)
                         ccGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttacgccgcatccgacat                <  1:209807‑M1/101‑66 (MQ=255)
                         ccGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttacgccgcatccgacat                <  2:1166489‑M1/101‑66 (MQ=255)
                         ccGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacat                >  1:839635‑M1/1‑36 (MQ=255)
                         ccGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgccgcatccgacat                <  2:1650175‑M1/101‑66 (MQ=255)
                          cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatc               >  2:1590649‑M1/1‑35 (MQ=255)
                          cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatc               >  2:1343812‑M1/1‑35 (MQ=255)
                          cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggttaaggcggttacgccgcatccgacatc               >  1:1440888‑M1/1‑35 (MQ=255)
                          cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgccgcatccgacatc               >  2:1583706‑M1/1‑35 (MQ=255)
                          cGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCAct               >  2:419751/1‑101 (MQ=255)
                          cGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCAct               >  1:1328399/1‑101 (MQ=255)
                           ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatct              <  2:949637‑M1/101‑68 (MQ=255)
                           ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatcc              <  2:759229‑M1/101‑68 (MQ=255)
                           ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatcc              <  2:624361‑M1/101‑68 (MQ=255)
                           ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatcc              >  2:737576‑M1/1‑34 (MQ=255)
                           ggTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCActt              >  1:1616747/1‑101 (MQ=255)
                            gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctgcaccccgtaggttggataaggcggttacgccgcatccgacatcta             <  2:264893‑M1/101‑69 (MQ=255)
                            gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttatgccgcatccgacatcta             >  1:41989‑M1/1‑33 (MQ=255)
                            gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatcta             <  2:564597‑M1/101‑69 (MQ=255)
                            gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatcta             <  2:361177‑M1/101‑69 (MQ=255)
                              tGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatctaac           >  2:1021887‑M1/1‑31 (MQ=255)
                              tGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatctaac           <  1:608991‑M1/101‑71 (MQ=255)
                              tGCCTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggtcggataaggcggttacgccgcatccgacatctaac           <  2:1241047‑M1/101‑71 (MQ=255)
                              tGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCActtctt           <  1:1446596/101‑1 (MQ=255)
                                ccTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatctaacgc         >  2:375608‑M1/1‑29 (MQ=255)
                                 cTGATGCGACGCTGGCGCGTCTTATCAGgcctacaccgctgtgaagtgctccaccccgtaggttggataaggcggttacgccgcatccgacatctaacgcc        <  2:17588‑M1/101‑74 (MQ=255)
                                 cTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCACTTCttttt        <  2:906719/101‑1 (MQ=255)
                                   gATAAGACGCTGGCGCGTCTTATCAGCATACACCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCACTTCTTTTTGa      >  2:760961/1‑101 (MQ=255)
                                      aaGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCACTTCTTTTTGAACg   >  1:929103/1‑101 (MQ=255)
                                      ||                                                                                                    
‑‑‑‑‑‑‑‑‑ACATCTAACGCCCGAGCCGGTTGCCTGATAAGACGCTGGCGCGTCTTATCAGCATACGCCACATCCCCGGCATACCATGCCGGATGTGGCACATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGG  >  173/1‑131

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