BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,938,374	C→T	D145N (GAC→AAC)	yiePs ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,938,374	0	C	T	99.8%	2139.3 / ‑5.9	576	D145N (GAC→AAC)	yiePs	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (307/267); total (309/267)
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.38e-01

Only 100 of 585 total aligned reads displayed.

GCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATA  >  NC_000913/3938284‑3938465

  gcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTg                                                                                         >  2:4024460/1‑93 (MQ=255)

   ccTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTgg                                                                                        >  2:698349/1‑93 (MQ=255)

     tGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTa                                                                                      <  1:6271605/93‑1 (MQ=255)

      gACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTaa                                                                                     >  1:6795252/1‑93 (MQ=255)

      gACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTaa                                                                                     <  1:3395233/93‑1 (MQ=255)

       aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAt                                                                                    <  1:4475033/93‑1 (MQ=255)

       aCAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAt                                                                                    <  1:3074043/93‑1 (MQ=255)

           gCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGAc                                                                                <  1:5544930/93‑1 (MQ=255)

             tttAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGt                                                                              <  1:1000355/93‑1 (MQ=255)

             tttAAATGCCGCGTCGCCATCGCTTTGGAGAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGt                                                                              >  2:6222111/1‑93 (MQ=255)

              ttAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTg                                                                             >  1:539380/1‑93 (MQ=255)

               tAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGa                                                                            >  2:3669131/1‑93 (MQ=255)

               tAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGa                                                                            >  2:6626695/1‑93 (MQ=255)

                aaaTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGGCGTGaa                                                                           <  1:822948/93‑1 (MQ=255)

                aaaTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaa                                                                           <  2:6019811/93‑1 (MQ=255)

                     ccGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaag                                                                      <  2:1313397/93‑1 (MQ=255)

                      cgcgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaagt                                                                     <  2:7573393/93‑1 (MQ=255)

                      cgcgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaagt                                                                     <  1:9158036/93‑1 (MQ=255)

                      cgcgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaagt                                                                     >  1:943078/1‑93 (MQ=255)

                        cgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAgtgt                                                                   <  1:1816058/93‑1 (MQ=255)

                        cgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAgtgt                                                                   <  2:9769864/93‑1 (MQ=255)

                        cgTCCCCATCGCTTTGGATAATCGCATCGACCATCGCCTTATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAgtgt                                                                   <  2:885035/93‑1 (MQ=255)

                          tCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGa                                                                 <  1:4034902/93‑1 (MQ=255)

                           cGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGAt                                                                <  1:8665024/93‑1 (MQ=255)

                             ccATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaa                                                              >  1:5164412/1‑93 (MQ=255)

                             ccATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaa                                                              <  2:3788615/93‑1 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                             <  1:3232586/93‑1 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGGCGTGAAGTAAGTGTGATAAAcc                                                           >  1:8037600/1‑93 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                           >  2:9006211/1‑93 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                           >  1:4577881/1‑93 (MQ=255)

                                   cTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGaa                                                        <  2:708391/93‑1 (MQ=255)

                                    tttGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAAt                                                       <  2:3317802/93‑1 (MQ=255)

                                     ttGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATg                                                      >  2:587211/1‑93 (MQ=255)

                                      tGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATgg                                                     >  1:2662541/1‑93 (MQ=255)

                                       ggATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGa                                                    >  1:418932/1‑93 (MQ=255)

                                        gATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGaa                                                   <  2:5008215/93‑1 (MQ=255)

                                        gATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGaa                                                   >  1:1403555/1‑93 (MQ=255)

                                          tAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATa                                                 <  2:7845784/93‑1 (MQ=255)

                                           aaTCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAg                                                >  2:1556637/1‑93 (MQ=255)

                                            aTCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGc                                               >  1:3529363/1‑93 (MQ=255)

                                               gCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAg                                            <  1:1218094/93‑1 (MQ=255)

                                               gCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAg                                            <  2:3218888/93‑1 (MQ=255)

                                                 aTCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGc                                          >  1:5812596/1‑93 (MQ=255)

                                                 aTCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGc                                          <  1:9868045/93‑1 (MQ=255)

                                                 aTCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGc                                          <  1:9644864/93‑1 (MQ=255)

                                                  tCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCg                                         >  2:2365945/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGGCGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                       >  2:223814/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                       <  1:7202066/93‑1 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                       >  1:5335724/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                       >  2:3479079/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                       >  1:7386693/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                       <  1:8474721/93‑1 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                       >  1:7879757/1‑93 (MQ=255)

                                                     aCAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAAt                                      <  2:724088/93‑1 (MQ=255)

                                                       aaTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGa                                    <  2:4444932/93‑1 (MQ=255)

                                                           gCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGGCGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTc                                >  2:4736609/1‑93 (MQ=255)

                                                           gCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTc                                >  2:8051229/1‑93 (MQ=255)

                                                             cTGATCCAGATCCAGCTTTATCACTGTGTTGCTGGAAATTGACGTGTAGTAAGTGTGATAAACCGAATGGTATAGCGAGGCGTATGATGTCaa                              >  2:9205699/1‑93 (MQ=255)

                                                              tGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGGCGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaa                             >  2:5101427/1‑93 (MQ=255)

                                                              tGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaa                             >  2:1881732/1‑93 (MQ=255)

                                                              tGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaa                             >  1:2431434/1‑93 (MQ=255)

                                                               gATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaaa                            >  1:4711219/1‑93 (MQ=255)

                                                                  gCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGTGAATGAGGTCGAAAACg                         >  1:8963902/1‑93 (MQ=255)

                                                                  gCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACg                         <  2:7660682/93‑1 (MQ=255)

                                                                    aGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGa                       <  2:5393123/93‑1 (MQ=255)

                                                                      aTCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAtt                     >  1:3858515/1‑93 (MQ=255)

                                                                      aTCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAtt                     <  2:4891985/93‑1 (MQ=255)

                                                                       tCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAACACGGATTg                    >  1:5677250/1‑93 (MQ=255)

                                                                       tCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTg                    >  2:8693775/1‑93 (MQ=255)

                                                                       tCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTg                    >  1:2344495/1‑93 (MQ=255)

                                                                        ccAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTgg                   >  2:4279731/1‑93 (MQ=255)

                                                                         cAGCTGTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGgc                  <  1:7660959/93‑1 (MQ=255)

                                                                           gCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGgcgc                >  2:9441094/1‑93 (MQ=255)

                                                                           gCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGgcgc                >  2:2091783/1‑93 (MQ=255)

                                                                            cTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCt               >  2:2689512/1‑93 (MQ=255)

                                                                             tttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTc              <  2:6445317/93‑1 (MQ=255)

                                                                             tttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGGGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTc              >  2:5595152/1‑93 (MQ=255)

                                                                              ttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCa             <  2:2552420/93‑1 (MQ=255)

                                                                              ttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCa             >  2:8679726/1‑93 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt            <  1:1594781/93‑1 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt            <  1:8820847/93‑1 (MQ=255)

                                                                                aTCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAtt           <  1:7101686/93‑1 (MQ=255)

                                                                                 tCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAttt          >  1:5916885/1‑93 (MQ=255)

                                                                                 tCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAttt          >  1:7011137/1‑93 (MQ=255)

                                                                                 tCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAttt          >  2:2395313/1‑93 (MQ=255)

                                                                                   aCTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCa        >  2:1459840/1‑93 (MQ=255)

                                                                                   aCTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCa        >  1:5028883/1‑93 (MQ=255)

                                                                                   aCTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCa        >  2:2929335/1‑93 (MQ=255)

                                                                                    cTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGGTCATTTCAt       >  1:6117059/1‑93 (MQ=255)

                                                                                    cTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCAt       >  1:2874196/1‑93 (MQ=255)

                                                                                     tgtgTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATa      <  1:3659923/93‑1 (MQ=255)

                                                                                     tgtgTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATa      >  2:687925/1‑93 (MQ=255)

                                                                                     tgtgTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATa      >  1:208509/1‑93 (MQ=255)

                                                                                     tgtgTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATa      >  1:6430502/1‑93 (MQ=255)

                                                                                      gtgtTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAg     <  1:664075/93‑1 (MQ=255)

                                                                                       tgtTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGa    <  1:9315360/93‑1 (MQ=255)

                                                                                       tgtTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGa    <  2:8457142/93‑1 (MQ=255)

                                                                                       tgtTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGa    >  1:7631106/1‑93 (MQ=255)

                                                                                       tgtTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGAGAATGAGGTCAACAACGGATTGGCGCTCATTTCATAGa    >  1:4636144/1‑93 (MQ=255)

                                                                                         tCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGata  >  1:7520187/1‑93 (MQ=255)

GCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGATA  >  NC_000913/3938284‑3938465

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