BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,865,007	A→G	T26A (ACC→GCC)	yeaG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,865,007	0	A	G	97.4%	360.2 / ‑0.1	117	T26A (ACC→GCC)	yeaG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (59/55); total (61/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 119 total aligned reads displayed.

AAAGGACGGCATACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTACCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATACAGCC  >  NC_000913/1864917‑1865096

aaaGGACGGCATACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGcc                                                                                         >  1:1658717/1‑93 (MQ=255)

aaaCGACGGCATACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGcc                                                                                         <  1:3281501/93‑1 (MQ=255)

  aGGACGGCATACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCAc                                                                                       <  1:4609580/93‑1 (MQ=255)

     aCGGCATACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTg                                                                                    <  2:3288234/93‑1 (MQ=255)

        ggaTACCATGAATATATTCGATCACTATCGCTAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCg                                                                                 <  1:732918/91‑1 (MQ=255)

        gCATACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCg                                                                                 >  1:2042642/1‑93 (MQ=255)

          aTACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGc                                                                               >  2:2057924/1‑93 (MQ=255)

          aTACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGc                                                                               >  1:3101451/1‑93 (MQ=255)

          aTACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGc                                                                               >  2:602416/1‑93 (MQ=255)

           tACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCa                                                                              >  1:4210986/1‑93 (MQ=255)

           tACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCa                                                                              >  2:1123976/1‑93 (MQ=255)

             ccATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAg                                                                            <  2:567453/93‑1 (MQ=255)

              caagAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGa                                                                           <  1:3657841/90‑1 (MQ=255)

              cATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGa                                                                           >  1:726723/1‑93 (MQ=255)

              cATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTACCACTTGTCGGCAAGa                                                                           >  2:4039188/1‑93 (MQ=255)

              cATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGGTTCTTGCCACTTGTCGGCAAGa                                                                           >  2:4142559/1‑93 (MQ=255)

              cATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGCTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGa                                                                           >  1:1396597/1‑93 (MQ=255)

               aTGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGAt                                                                          <  1:3298222/93‑1 (MQ=255)

               aTGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGAt                                                                          <  2:2430988/93‑1 (MQ=255)

                tGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATc                                                                         >  2:604498/1‑93 (MQ=255)

                 gAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCg                                                                        >  1:4078406/1‑93 (MQ=255)

                 gAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCg                                                                        >  2:777017/1‑93 (MQ=255)

                 gAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGCTCACACTGCAGGAGCTTCTTGCCACTTGTCGGCAAGATCg                                                                        >  1:70861/1‑93 (MQ=255)

                    tatatTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAg                                                                     <  1:3201107/93‑1 (MQ=255)

                     atatTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGt                                                                    <  1:1925815/93‑1 (MQ=255)

                      tatTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTg                                                                   <  2:1570268/93‑1 (MQ=255)

                      tatTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTg                                                                   <  2:648380/93‑1 (MQ=255)

                        ttCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCt                                                                 >  1:997013/1‑93 (MQ=255)

                         tCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCtt                                                                >  2:639169/1‑93 (MQ=255)

                          cGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTa                                                               <  1:2800330/93‑1 (MQ=255)

                          cGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTa                                                               <  2:2966160/93‑1 (MQ=255)

                          cGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTa                                                               >  1:1695170/1‑93 (MQ=255)

                           gATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTAt                                                              >  2:363426/1‑93 (MQ=255)

                            aTCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATg                                                             >  2:787941/1‑93 (MQ=255)

                             tCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGc                                                            >  2:331161/1‑93 (MQ=255)

                             tCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGc                                                            >  2:1311232/1‑93 (MQ=255)

                              cACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCCCACTGCAGGAGCTTCTTGCCACTTGTCGGCAAGATCGCAGCGCTTATGcc                                                           >  2:108740/1‑93 (MQ=255)

                              cACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGcc                                                           >  1:1717938/1‑93 (MQ=255)

                               aCTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCa                                                          >  1:2045593/1‑93 (MQ=255)

                                cTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCaa                                                         <  2:2144012/93‑1 (MQ=255)

                                cTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCaa                                                         <  1:995752/93‑1 (MQ=255)

                                  aTCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAAcg                                                       >  2:304189/1‑93 (MQ=255)

                                  aTCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAAcg                                                       >  1:3026314/1‑93 (MQ=255)

                                   tCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAAcgc                                                      >  2:4574563/1‑93 (MQ=255)

                                   tCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAAcgc                                                      >  2:4574564/1‑93 (MQ=255)

                                   tCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAAcgc                                                      <  1:374031/93‑1 (MQ=255)

                                    cGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAAcgcg                                                     <  1:2467371/93‑1 (MQ=255)

                                     gCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCgg                                                    >  2:2734454/1‑93 (MQ=255)

                                     gCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTACCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCgg                                                    >  2:1563664/1‑93 (MQ=255)

                                         gCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGa                                                <  1:1143162/93‑1 (MQ=255)

                                          cGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAg                                               >  2:2379987/1‑93 (MQ=255)

                                           gATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGc                                              <  1:3583939/93‑1 (MQ=255)

                                           gATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGc                                              >  2:4152589/1‑93 (MQ=255)

                                           gATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGc                                              >  1:1285742/1‑93 (MQ=255)

                                            atatGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCg                                             <  1:546021/93‑1 (MQ=255)

                                            atatGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCg                                             >  1:4596837/1‑93 (MQ=255)

                                            atatGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCg                                             >  1:3618932/1‑93 (MQ=255)

                                              atGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGc                                           <  2:4389556/93‑1 (MQ=255)

                                               tGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCt                                          <  2:3216931/93‑1 (MQ=255)

                                               tGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCt                                          >  1:927468/1‑93 (MQ=255)

                                               tGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCt                                          >  1:346423/1‑93 (MQ=255)

                                                 aaGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTAt                                        >  1:2799797/1‑93 (MQ=255)

                                                  accTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTAtt                                       <  1:4615394/91‑1 (MQ=255)

                                                  aGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTAtt                                       >  2:2973087/1‑93 (MQ=255)

                                                  aGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTAtt                                       <  2:1803811/93‑1 (MQ=255)

                                                  aGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTAtt                                       <  1:2692123/93‑1 (MQ=255)

                                                     tGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGAt                                    <  1:101342/93‑1 (MQ=255)

                                                        cAAGTACGTAGAGTACACACTTCAGGAGTTTCTTGCTTCCTGCCGGCAAGAACGTAGTGCCTTTGCCAACGCGGCTGAGCGGCTATTGATgtc                                 >  2:4332721/1‑91 (MQ=255)

                                                         aaGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCt                                <  2:2442366/93‑1 (MQ=255)

                                                         aaGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCt                                <  1:1440742/93‑1 (MQ=255)

                                                         aaGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCt                                >  1:283242/1‑93 (MQ=255)

                                                         aaGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCt                                <  1:305468/93‑1 (MQ=255)

                                                          aGGACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTa                               <  2:1343357/93‑1 (MQ=255)

                                                            gACGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATc                             <  2:2143521/93‑1 (MQ=255)

                                                              cGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCgg                           <  1:4146346/93‑1 (MQ=255)

                                                              cGAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCgg                           <  1:4462030/93‑1 (MQ=255)

                                                               gAAGAGTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGt                          >  1:1985164/1‑93 (MQ=255)

                                                                   agTTCACACTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGc                      <  1:322881/93‑1 (MQ=255)

                                                                       cacaCTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGt                  >  1:4059906/1‑93 (MQ=255)

                                                                       cacaCTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGt                  <  1:2479723/93‑1 (MQ=255)

                                                                        acacTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTc                 >  2:204692/1‑93 (MQ=255)

                                                                         cacTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCa                <  1:3042379/93‑1 (MQ=255)

                                                                          acTGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCAt               <  1:4203020/93‑1 (MQ=255)

                                                                          acTGCAGGAGTTTCTTACCGCTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCAt               <  1:3826792/93‑1 (MQ=255)

                                                                            tGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATgg             <  1:1398310/93‑1 (MQ=255)

                                                                            tGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATgg             <  1:1903727/93‑1 (MQ=255)

                                                                            tGCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATgg             <  1:1620896/93‑1 (MQ=255)

                                                                             gCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGt            <  2:3313815/93‑1 (MQ=255)

                                                                             gCAGGAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGt            <  1:1601310/93‑1 (MQ=255)

                                                                                ggAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGa         >  1:4144780/1‑93 (MQ=255)

                                                                                 gAGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGAt        >  1:3092432/1‑93 (MQ=255)

                                                                                  aGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATa       <  2:1368610/93‑1 (MQ=255)

                                                                                  aGTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATa       >  2:3710014/1‑93 (MQ=255)

                                                                                   gTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATAc      <  2:1385802/93‑1 (MQ=255)

                                                                                   gTTTCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATAc      >  2:675859/1‑93 (MQ=255)

                                                                                    tttCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATACa     <  2:2998114/93‑1 (MQ=255)

                                                                                    tttCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATACa     <  2:2998117/93‑1 (MQ=255)

                                                                                     ttCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATACAg    >  1:3287227/1‑93 (MQ=255)

                                                                                     ttCTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATACAg    <  2:2540705/93‑1 (MQ=255)

                                                                                       cTTGCCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATACAGcc  >  1:3884269/1‑93 (MQ=255)

AAAGGACGGCATACCATGAATATATTCGATCACTATCGCCAGCGATATGAAGCTGCCAAGGACGAAGAGTTCACACTGCAGGAGTTTCTTACCACTTGTCGGCAAGATCGCAGTGCTTATGCCAACGCGGCTGAGCGGCTATTGATGGCTATCGGTGAGCCTGTCATGGTCGATACAGCC  >  NC_000913/1864917‑1865096

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