Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,871,089	G→A	6.2%	A235A (GCG→GCA)	yeaI →	putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,871,089	0	G	A	6.2%	198.4 / 2.5	65	A235A (GCG→GCA)	yeaI	putative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base A (2/2);  total (32/33)
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.35e-01

CTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTC  >  NC_000913/1870943‑1871234                                                                                                                                                   |                                                                                                                                                  cTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGa                                                                                                                                                  <  7:168401/148‑1 (MQ=255)     acaacaATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTgg                                                                                                                                             >  8:312879/1‑149 (MQ=255)      caacaaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTggg                                                                                                                                            <  1:430239/149‑1 (MQ=255)         caaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCTCACAATCTGAGCAGTTATAAAGCTGACTATTCTTTGTATGTGGTCGATTACTGGCCGGACAACCATACTGCGACCTGGGTaa                                                                                                                                         >  5:127070/1‑149 (MQ=255)          aaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCAACCTGGGGAAt                                                                                                                                        >  1:5644/1‑149 (MQ=255)           aTAAGAAAAAAACGGGAATCATGTTGCTTGAGCTTATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTAGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATc                                                                                                                                       <  3:322830/149‑1 (MQ=255)            tAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATc                                                                                                                                       >  7:252970/1‑148 (MQ=255)            tAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGAATGTCGTCGATTACTGTCGCGAGAACGATTCTGCGACCTggggctgc                                                                                                                                       >  2:404514/1‑144 (MQ=255)              aGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCaa                                                                                                                                     <  4:122675/148‑1 (MQ=255)              aGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAAc                                                                                                                                    >  7:221852/1‑149 (MQ=255)                aaaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACta                                                                                                                                  <  7:312879/149‑1 (MQ=255)                          gAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAata                                                                                                                        >  8:99118/1‑149 (MQ=255)                           caTCCTGTTTCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGACCAGTTATAACGCTGACTATTCTTTGTATGGCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAatat                                                                                                                       <  6:127070/148‑1 (MQ=255)                            aTCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCTGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATAtt                                                                                                                      <  8:221852/149‑1 (MQ=255)                            aTCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCTGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATAtt                                                                                                                      <  2:26917/149‑1 (MQ=255)                              ccTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTgg                                                                                                                    >  5:120913/1‑149 (MQ=255)                               cTGGTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGGTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGt                                                                                                                   <  4:347271/149‑1 (MQ=255)                                    gCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCATCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTc                                                                                                              >  4:306202/1‑149 (MQ=255)                                            tGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCt                                                                                                             <  4:5244/142‑1 (MQ=255)                                            tGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCt                                                                                                             >  3:5244/1‑142 (MQ=255)                                               tCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAttt                                                                                                   >  4:247471/1‑149 (MQ=255)                                                cccTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAtttt                                                                                                  >  5:105186/1‑149 (MQ=255)                                                cccTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAtttt                                                                                                  >  1:40644/1‑149 (MQ=255)                                                   ttttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                                                                                               >  4:361183/1‑149 (MQ=255)                                                   ttttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGAACTGGCGAATCAACTATACAAATATATTGGTTTGTCTGTGGGCATTTTTAc                                                                                               >  5:333835/1‑149 (MQ=255)                                                     ttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCAACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACt                                                                                              <  2:5644/148‑1 (MQ=255)                                                           tttttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCtt                                                                                        >  7:351298/1‑148 (MQ=255)                                                             tttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTtta                                                                                      <  3:306202/148‑1 (MQ=255)                                                             tttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTtta                                                                                      <  4:311316/148‑1 (MQ=255)                                                              ttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTttat                                                                                     <  3:361183/148‑1 (MQ=255)                                                               tCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAATCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTttatt                                                                                    <  3:106870/148‑1 (MQ=255)                                                                ccccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTtattat                                                                                  >  5:80881/1‑149 (MQ=255)                                                                   cTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCAt                                                                               >  3:306339/1‑149 (MQ=255)                                                                    ttCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTc                                                                                                              <  4:149509/117‑1 (MQ=255)                                                                    ttCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTc                                                                                                              >  3:149509/1‑117 (MQ=255)                                                                         gCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGca                                                                          <  8:252970/148‑1 (MQ=255)                                                                          cacacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGcaca                                                                        <  6:80881/149‑1 (MQ=255)                                                                            cacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCacac                                                                       >  6:311000/1‑148 (MQ=255)                                                                                          ttATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTg                                                                                                                     >  3:148034/1‑88 (MQ=255)                                                                                          ttATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTg                                                                                                                     <  4:148034/88‑1 (MQ=255)                                                                                             taACGCTGACTCTTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATTCTCACACGATTAGCCAGCGAACTAt                                                     <  2:40644/149‑1 (MQ=255)                                                                                               aCGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATg                                                    <  8:351298/148‑1 (MQ=255)                                                                                               aCGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATg                                                    <  3:247471/148‑1 (MQ=255)                                                                                                cGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGACCTATGGc                                                  >  8:84826/1‑149 (MQ=255)                                                                                                     aCTATTCTTTGTATGTCGTCGATTACTGTCCGGTCAACCATACTGGGATCTGTGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCggt                                              <  8:423656/148‑3 (MQ=255)                                                                                                     aCTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGtt                                              >  2:86492/1‑148 (MQ=255)                                                                                                          tCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAg                                         >  1:170871/1‑148 (MQ=255)                                                                                                            tttGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCat                                      >  7:215370/1‑149 (MQ=255)                                                                                                              tGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGATTCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTTTCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCa                                       <  1:404514/146‑1 (MQ=255)                                                                                                              tGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCa                                       <  5:408791/146‑1 (MQ=255)                                                                                                                tATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCattatt                                  <  7:99118/149‑1 (MQ=255)                                                                                                                    tcgtcgATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGt                              >  7:110082/1‑149 (MQ=255)                                                                                                                    tcgtcgATTACTGTCCGGACAACCATACTGCAACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATa                                          <  7:259593/137‑1 (MQ=255)                                                                                                                    tcgtcgATTACTGTCCGGACAACCATACTGCAACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATa                                          >  8:259593/1‑137 (MQ=255)                                                                                                                           ttACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCAt                       <  6:120913/149‑1 (MQ=255)                                                                                                                               tGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGc                    >  8:258211/1‑148 (MQ=255)                                                                                                                                   cGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTAtg                <  8:110082/148‑1 (MQ=255)                                                                                                                                    ggACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTAtgct              <  3:253814/149‑1 (MQ=255)                                                                                                                                        aaCCATACTGCTACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTCGGCATTTTTACTGTTCTTTATTATAATGGGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTAAGCTGCaa          >  4:373577/1‑149 (MQ=255)                                                                                                                                        aaCCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCaa          >  6:380118/1‑149 (MQ=255)                                                                                                                                           cATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTGTCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTaa                                            >  5:23065/1‑112 (MQ=255)                                                                                                                                           cATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTGTCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTaa                                            <  6:23065/112‑1 (MQ=255)                                                                                                                                               cTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAacttact    <  6:105186/148‑1 (MQ=255)                                                                                                                                               cTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGAATTATTATGTCTGGCATCGCTATGCTGCAacttact    >  1:173054/1‑148 (MQ=255)                                                                                                                                                  cGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTc  <  7:258211/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACTATGGCCGTTAATAGCATTATTATGTCTGGCATCGCTATGCTGCAACTTACTTC  >  NC_000913/1870943‑1871234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑