Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AC_000091	526,302	(G)6→7	coding (3818/4281 nt)	rhsD →	rhsD element protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AC_000091	526,296	1	.	G	98.7%	296.7 / ‑2.2	75	coding (3812/4281 nt)	rhsD	rhsD element protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (48/26);  total (49/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

TATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑GGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAA  >  AC_000091/526209‑526380                                                                                          |                                                                                     tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gg                                                                                    >  2:1443833/1‑90 (MQ=255) tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gg                                                                                    >  1:1566385/1‑90 (MQ=255) tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gg                                                                                    >  1:1700025/1‑90 (MQ=255)  atatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑ggg                                                                                   <  1:1063354/90‑1 (MQ=255)  atatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑ggg                                                                                   >  2:2243730/1‑90 (MQ=255)   tatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggg                                                                                  <  2:1227816/90‑1 (MQ=255)   tatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggg                                                                                  <  1:2187260/90‑1 (MQ=255)    atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑ggggg                                                                                 >  2:349193/1‑90 (MQ=255)     tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggggg                                                                                <  1:111582/90‑1 (MQ=255)     tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggggg                                                                                >  1:1202211/1‑90 (MQ=255)     tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggggg                                                                                >  2:2755251/1‑90 (MQ=255)     tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggggg                                                                                >  1:1993321/1‑90 (MQ=255)      cAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑ggggggg                                                                               <  1:563648/90‑2 (MQ=255)        gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggggggag                                                                             <  2:954361/90‑4 (MQ=255)        gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggggggag                                                                             >  1:3033211/1‑87 (MQ=255)        gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCAAGATCTGGAGCATGTA‑C‑gggggggag                                                                             >  1:1913536/1‑87 (MQ=255)         tATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑gggggggagt                                                                            >  2:611893/1‑86 (MQ=255)            ccTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTtg                                                                          >  1:1852747/1‑90 (MQ=255)             cTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTtgt                                                                         <  2:899347/90‑1 (MQ=255)             cTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTtgt                                                                         >  1:224022/1‑90 (MQ=255)              tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTtgtg                                                                        >  1:1333107/1‑90 (MQ=255)                tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTgtggt                                                                      >  1:478552/1‑90 (MQ=255)                 aaaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTgtggtg                                                                     >  2:884352/1‑90 (MQ=255)                  aaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGgtgt                                                                    >  1:2077956/1‑90 (MQ=255)                  aaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGgggt                                                                    >  2:2854270/1‑90 (MQ=255)                   aTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGtt                                                                   >  2:3449295/1‑90 (MQ=255)                   aTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGtt                                                                   >  1:2977215/1‑90 (MQ=255)                    tCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTc                                                                  >  1:3083648/1‑90 (MQ=255)                     ccACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCt                                                                 >  1:1045366/1‑90 (MQ=255)                     ccACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCt                                                                 <  2:317674/90‑1 (MQ=255)                      cACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCtt                                                                <  2:1333658/90‑1 (MQ=255)                       aCTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGGGTTTGTGGTGTTCttt                                                               >  2:2459347/1‑90 (MQ=255)                        cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTc                                                              >  2:1090617/1‑90 (MQ=255)                          acaacaAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCAc                                                            >  2:2118117/1‑90 (MQ=255)                          acaacaAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCAc                                                            >  1:1767416/1‑90 (MQ=255)                           caacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACg                                                           >  2:22423/1‑90 (MQ=255)                           caacaaATTGACCCTATGCGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACg                                                           <  1:1284024/90‑1 (MQ=255)                             acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACAta                                                         <  2:2211158/90‑1 (MQ=255)                               aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGtata                                                       <  1:2010428/90‑1 (MQ=255)                                aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTataa                                                      >  2:1576439/1‑90 (MQ=255)                                aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTataa                                                      >  2:557393/1‑90 (MQ=255)                                 aTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTataat                                                     >  1:896145/1‑90 (MQ=255)                                 aTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTataat                                                     >  2:1192559/1‑90 (MQ=255)                                  ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTataata                                                    >  2:2335723/1‑90 (MQ=255)                                   tGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAg                                                   <  2:3028352/90‑1 (MQ=255)                                    gACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAgg                                                  >  2:431821/1‑90 (MQ=255)                                     aCCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGa                                                 <  1:245597/90‑1 (MQ=255)                                     aCCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGa                                                 >  2:344066/1‑90 (MQ=255)                                       ccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGAcc                                               <  1:1412538/90‑1 (MQ=255)                                       ccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGAcc                                               >  2:3455078/1‑90 (MQ=255)                                       ccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGAcc                                               <  2:344346/90‑1 (MQ=255)                                        cTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCa                                              >  1:2656418/1‑90 (MQ=255)                                        atATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCa                                              <  1:2375348/89‑1 (MQ=255)                                         tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTATG‑GGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCaa                                             >  1:2511156/1‑90 (MQ=255)                                          aTGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACTAAg                                            >  2:1865227/1‑90 (MQ=255)                                          aTGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAg                                            >  1:3297067/1‑90 (MQ=255)                                          aTGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAg                                            <  1:1623592/90‑1 (MQ=255)                                           tGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGt                                           <  1:1982557/90‑1 (MQ=255)                                            gggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                          <  2:3455121/90‑1 (MQ=255)                                              gATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                        >  2:203863/1‑90 (MQ=255)                                              gATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                        >  1:1908618/1‑90 (MQ=255)                                                ttattGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                      <  2:3167145/90‑1 (MQ=255)                                                  attGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTg                                    >  1:2492761/1‑90 (MQ=255)                                                    tGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGAt                                  >  1:1816668/1‑90 (MQ=255)                                                     gCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                 <  1:1267221/90‑1 (MQ=255)                                                      cAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                >  2:2762760/1‑90 (MQ=255)                                                       aGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGt                               <  1:3010676/90‑1 (MQ=255)                                                       aGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGt                               >  2:2764024/1‑90 (MQ=255)                                                       aGACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGt                               >  1:3254906/1‑90 (MQ=255)                                                        gACTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                              >  2:3034662/1‑90 (MQ=255)                                                         aCTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                             <  2:747814/90‑1 (MQ=255)                                                          cTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                            >  2:1742844/1‑90 (MQ=255)                                                          cTTGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                            >  1:3385444/1‑90 (MQ=255)                                                           ttGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                           <  1:486585/90‑1 (MQ=255)                                                           ttGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                           >  2:2225437/1‑90 (MQ=255)                                                            tGGGATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGc                          >  2:2805791/1‑90 (MQ=255)                                                              ggATGATGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                        >  1:2039430/1‑90 (MQ=255)                                                               gatgatGCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGa                       >  1:2390159/1‑90 (MQ=255)                                                                   atgCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                   <  2:3259139/90‑1 (MQ=255)                                                                   atgCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                   >  2:219851/1‑90 (MQ=255)                                                                    tgCCAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctg                  <  1:1213/90‑1 (MQ=255)                                                                       cAGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGt               <  2:577803/90‑1 (MQ=255)                                                                        aGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGtt              <  2:1126957/90‑1 (MQ=255)                                                                        aGATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGtt              <  1:186949/90‑1 (MQ=255)                                                                         gATCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa             >  1:863456/1‑90 (MQ=255)                                                                          aTCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg            <  2:1626499/90‑1 (MQ=255)                                                                           tCTGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa           >  2:1536678/1‑90 (MQ=255)                                                                             tGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg         >  2:2984905/1‑90 (MQ=255)                                                                             tGGAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATg         >  1:2695622/1‑90 (MQ=255)                                                                               gAGCATGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCt       >  1:925558/1‑90 (MQ=255)                                                                                   aTGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTga   >  1:1192913/1‑90 (MQ=255)                                                                                    tGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaa  <  1:1762069/90‑1 (MQ=255)                                                                                    tGTA‑CGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaa  >  2:2135499/1‑90 (MQ=255)                                                                                          |                                                                                     TATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTA‑C‑GGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAA  >  AC_000091/526209‑526380

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